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INTRODUCTION 
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more  advanced,  higher  performance  infrared  countermeasure  systems. 

This  edition  of  The  DTIC  Review  features  the  latest  requirements,  developments  and  technologies  in  advanced 
infrared  countermeasures  systems.  This  issue  goes  on  further  to  examine  current,  evolving  and  future  threats 
associated  with  infrared  countermeasures. 
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1  Introduction 


This  report  summarizes  the  results  obtained  in  a  three- year  investigation  into  the  application  of 
wavelet  methods  in  modulation  techniques  for  multi-user  spread  spectrum  communication  systems. 
The  investigation  was  undertaken  between  November  1992  and  January  1996. 

2  Objectives 

The  objectives  of  this  project  remained  essentially  as  outlined  in  the  project  proposal  and  pre- 
vious  status  reports: 

•  Task  1:  Development  of  and  evaluation  of  wavelet  symbols  for  multiple  access  spread  spectrum 
communications. 

•  Task  2;  Evaluation  of  communication  system  performance, 

•  Task  3:  Investigation  of  possible  countermeasures  to  wavelet-based  spread  spectrum. 

In  addition  to  addressing  these  objectives,  some  basic  research  on  wavelet  analysis  was  performed 
in  connection  with  the  project  which  led  to  a  class  of  multidimensional  generalizations  of  the  wavelet 
transform.  ^ 

3  Status  of  Effort 

The  project  was  completed  in  January,  1996  with  the  benefit  of  a  two-month  extension  of  the 
original  three-year  time  line. 

Task  1  led  to  development  and  analysis  of  a  class  of  bandlimited  orthonormal  wavelet  sets  having 
properties  well  suited  for  use  with  the  wavelet-based  scale-division  multiple  access  communication 
scheme  around  which  much  of  this  project  was  based. 

Crosstalk  effects  in  scale-division  multiple  access  under  asynchronous  and  distribute*^^ operation 
were  studied  under  Task  2.  Other  performance  questions  originally  planned  to  be  addressed  under 
this  Task  were  de-emphasized  in  later  stages  of  the  project  in  view  of  published  work  on  performance 
of  wavelet-based  communication  protocols  by  other  researchers. 

Task  3  led  to  introduction  and  evaluation  of  an  algorithm  for  detection  of  cyclostationary  signals. 
This  detector  was  used  to  evaluate  of  the  detectability  of  scale- division  multiple  access  signals  by 
cyclostationary  methods. 
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4  Accomplishments  /  New  Findings 

At  the  inception  of  this  project,  only  a  handful  of  researchers  had  looked  into  the  possibility  of 
using  wavelets  [8,  44,  49,  50]  or  fractal  signals  [57,  58]  in  the  modulation  of  communication  signals. 
The  research  group  at  Arizona  State  University  was  the  first  to  investigate  multiple  access  issues 
in  this  context  [8].  There  are  now  over  forty  published  papers  in  this  area  [2,  5,  9-13,  17-21,  24-30, 
33-36,  38-45,  47,  49-50,  53,  55-59]. 

Work  under  this  effort  focused  on  wavelet-based  multiple-access  techniques  in  which  digital 
communication  signals  are  encoded  on  orthogonal  w'avelet  sets,  as  was  introduced  in  [8]  and  described 
in  the  original  proposal.  Highlights  of  the  results  obtained  in  this  project  are  summarized  below. 

4.1  Orthonormal  Wavelet  Symbols 

Scale-division  multiple  access  (SDMA)  refers  to  a  multiple-access  communication  scheme  in 
which  the  orthogonal  symbols  are  obtained  by  dilation  and  time  shifting  of  a  single  prototype 
wavelet.  The  various  users’  messages  are  separated  into  channels  based  on  the  scale  of  dilation  of 
the  wavelet  symbol  on  which  their  messages  are  encoded  -  hence  the  term  SDM.4  [8,  10,  52].  As 
shown  in  these  references,  the  spectral  structure  of  a  transmitted  SDMA  signal  depends  on  the 
choice  of  the  prototype  wavelet  symbol  in  a  straightforward  way. 

The  goals  in  construction  of  wavelets  for  use  as  symbols  in  SDMA  are  different  from  those  en¬ 
countered  in  designing  wavelets  for  other  applications.  For  use  in  time- frequency  analysis  of  signals, 
for  example,  it  is  generally  desirable  for  wavelets  to  have  their  energy  simultaneously  concentrated 
in  both  the  time  and  frequency  domains  in  order  to  provide  signal  analysis  that  is  localized  in  the 
phase  plane.  In  spread  spectrum  applications,  however,  wavelets  should  be  compactly  supported  or 
highly  concentrated  in  the  time  domain  while  having  broadband  frequency  structure. 

In  general,  SDMA  relies  on  exact  synchronization  of  all  transmitters  in  the  network  to  maintain 
orthogonality  of  the  signals  and  avoid  channel  crosstalk.  Even  if  this  is  achieved,  time  delays 
introduced  by  propagation  through  the  transmission  medium  will  cause  the  received  signals  to  be 
non-orthogonal.  This  issue  may  be  addressed  by  requiring  that  the  channels  in  the  network  use 
symbols  with  Fourier  transforms  having  disjoint  supports.  In  this  project,  orthonormal  wavelet 
bases  generated  by  mother  wavelets  whose  frequency  spectra  are  supported  in  several  disjoint  bands 
were  constructed. 

The  simplest  example  of  a  dyadic  orthonormal  basis  of  bandlimited  wavelets  is  generated  by  the 
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mother  wavdet  defined  in  the  frequency  domain  by 


h((jj)  = 


1  TT  <  |u;|  <  2tk 
0  otherwise 


(1) 


with  ao  =  2  and  6o  =  1-  A  related,  but  more  complicated,  example  of  an  orthonormal  basis  of 
bandlimited  wavelets  was  described  by  Mallat  using  the  mother  wavelet  defined  by 


h{ijj) 


1  Iwj  €  (=7r,7r)U[4x,^7r) 
0  otherwise 


(2) 


Examination  of  the  relationship  between  these  two  wavelet  bases  gives  insight  about  how  to  construct 
other  bandlimited  wavelet  bases  in  which  the  dilations  of  the  mother  wavelet  h  do  not  overlap  in 
the  frequency  domain. 

In  what  follows,  the  connected  components  of  the  support  of  a  mother  w^avelet’s  Fourier  transform 
h  wiU  be  called  ‘‘slots.”  Throughout  this  section  h  will  be  a  real- valued  wavelet,  which  implies  that 
A(— u;)  =  and  hence  the  slots  are  symmetric  about  ij  =  0.  Accordingly,  when  the  values  or 

support  of  h  are  specified  on  the  positive  frequency  axis,  they  will  be  assumed  to  be  given  on  the 
negative  frequency  axis  by  symmetry.  The  class  of  wavelets  to  be  constructed  will  be  categorized 
according  to  the  number  of  slots  they  have  in  the  positive  frequency  axis.  The  wavelet  defined  by 
(1)  is  thus  a  ‘‘one-slot“  wavelet  while  the  one  defined  by  (2)  is  a  “two-slot”  wavelet.  The  following 
sections  describe  the  construction  of  “n-slot”  orthonormal  wavelet  bases. 

Suppose  the  support  So  of  h  is  the  union  of  n  disjoint  intervals 


So  =  [-So,  *51]  U  [S2-  '53]  U  •  •  •  U  [-Sn-l ,  Sn] 


Then  the  support  of  is 


5,n  =  [so2^^.  5i2-^]  U  (s22“^,  532“^]  U  •  •  *  U  ,  3^2^^] 


If  the  measure  of  5m  H  5^,1^  is  zero  for  all  non-equal  integers  m  and  m',  5o  will  be  called  an  orthogonal 
support.  In  particular,  an  orthogonal  support  consisting  of  n  disjoint  intervals  will  be  called  an  n-slot 
orthogonal  support. 

K  ||/i||  =  1  and  the  support  of  h  is  orthogonal, 

for  all  m  7^  m'.  In  this  case,  for  the  time-shifted  and  dilated  replicates  /i^  of  h  to  form  an  orthonor¬ 
mal  set  in  it  remains  only  to  ensure  that  h^  and  are  orthonormal  for  n  /  n'  and  all  m. 
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This  is  simplified  by  the  observation  that 

{h^ra^ht)  =  {hlX) 

Thus,  one  must  only  verify  that  all  of  the  time-shifted  replicates  of  h  are  orthonormal  at  a  single 
level  of  dilation. 

4.1.1  Construction  of  Orthogonal  Supports 
Consider  rfeaJ  numbers  0  <  cq  <  *  •  •  <  Cn  =  2co.  Then 

[co,2co]  =  [co,Ci]  U  [Ci,C2]  U  [C2,C3]U  -  U  [c„_i,2co] 

Let  be  distinct  integers  and  dilate  each  sub-interval  [c_,, Cj+i]  by  2^-'+'  to  form  a  set 

5o  =  [co2P^Cl2P>]u[cl2P^C22P^]U•••U[cn-l2P^2co2P"] 

=  [-30,51]  u  [52,53]  u  •• -U  [s2„,S2n-l] 

Then  So  is  clearly  an  n-slot  orthogonal  support.  Conversely,  suppose  So  is  an  n-slot  orthogonal 
support  of  the  form  (3).  Then  it  is  straightforward  to  show,  that  So  can  be  generated  by  the  method 
just  described  with  possibly  some  subintervals  deleted, 

4.1.2  Orthonormality  of  the  Wavelets 

Suppose  ||/i||  =  1  and  the  support  of  h  is  orthogonal.  Then  {h^}  is  an  orthonormal  wavelet  set 
if  and  only  if  {Hq,  Hq)  =  0  for  all  integers  n  and  n'  with  n  7^  n'.  This  is  equivalent  to 

0  =  f  h[t  —  n)h"{t  —  n')  dt 

2tV  Jf,  ZTT 

for  all  n  7^  n'.  With  fc  =  n  —  n'  7^  0,  this  expression  becomes 

/  du;  =  0  (4) 

Jn 

for  any  non-zero  integer  k.  The  following  section  presents  some  orthonormal  wavelet  sets  corre¬ 
sponding  to  certain  solutions  of  this  equation. 
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4.1.3  Flat  n-Siot  Orthonormal  Wavelets 


The  first  n-slot  orthonormal  wavelets  constructed  were  similar  to  the  wavelets  defined  by  (i) 
and  (2)  in  that  their  Fourier  transforms  are  constant  on  the  n  slots  and  zero  elsewhere.  These  will 
be  called  flat  n-slot  wavelets.  Construction  of  "non-flat”  n-slot  orthonormal  wavelets  is  undertaken 
in  the  next  section  of  this  report. 

Let 


V  =  [co2P^ ,  ci2P^]  U  [cl2P^  022^^]  U  •  •  •  U  Cn2^-] 

where  0  <  co  <  . . .  <  Cn  =  2co  and  •  •  • » Pn  are  integers  and  define 


a  = 


Then  if  h{u)  =  a  for  u;  €  V  and  h{u;)  =  0  otherwise,  ||ft||  =  1  and  equation  (4)  reduces  to 


7  e''^‘^\h{u 

/IS, 

=  «'(/ 

•-co2fl 

fCi2Pi 

+ 

/ 

J- 

-cr2'>i 

Jca^n 

=  2a^( 

J 

r-'i  2Pi 

rC3  2P2 

/  + 

/ 

'co2Pl 

/ci2P2 

2ry^  . 

—  ■■  ■"  1 

—  sin(co2^^A:)  -f 

Cn2P^ 


-  sin(Cn-i2^''fc)  -f  sin(Cn2^”A:)}  =  0 


for  all  A:  7^  0.  This  is  equivalent  to 

-  sin(co2^^  A;)  +  sin(ci2^^A:)  -  sin(ci2^-A:)  -f  sin(c22^-A:)  - - sin(cn-i2^^A:)  +  sin(cn2^"'A’)  =  0  (5) 

for  all  A:  0.  If  Co2^S  ,  Ci2p^,  022^^,  . . c^_i2P”,  and  satisfy  equation  (5),  then  so  will 

iV/co2^^,  Mci2P^,  iV/ci2P2^  Mc22^^^  . . and  Mcn2^^  for  any  natural  number  M.  Thus, 
there  exists  a  support  of  lowest  dilation  level  which  wiU  be  called  the  mother  support.  In  other 
words,  if  the  end  points  of  the  support 


V  =  [co2^\c^2^^]  U  [ci2PSc22P2]  U  •  -  •  U  [cn^^2^^ 0^2^-] 
satisfy  (5)  and  the  end  points  of  the  dilated  support 

[co2^^/M,Ci2^^/M]  U  [ci2PViV/,C22PV^-^]  U  •  -  U  [cn^i2P^ / M ,  €^2^- / M] 


do  not  satisfy  (5)  for  any  M  €  N,  then  V  is  the  mother  support.  • 
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Example  1  In  the  one-slot  case,  the  support  consists  of  only  one  interval  [co2^^ ,  2co2^^]  with  cq  >  0. 
Denoting  a  =  equation  (5)  becomes  -sin(afc)  +  sin(2aA;)  =  0,  or 

sin(Gfc)[2cos(aA:)  -  1]  =  0 

The  only  non- trivial  solution  of  this  equation  that  holds  for  any  non-zero  integer  k  is  a  =  m-n  with 
m  €  Hence,  the  support  is  [m7r,  2m7r],  the  mother  support  is  [^,27r],  and  an  orthonormal  wavelet 
basis  obtained  having  the  mother  wavelet  defined  in  equation  (1). 

Example  2  In  the  two-slot  case,  the  cut  points  cq  and  Ci  and  the  dilation  p  should  satisfy 

—  sin(coAr)  -I-  sin(cifc)  -  sin(ci2^A:)  -f  sin(2co2^fc)  =  0 


This  equation  can  be  decomposed  into  two  equations  (with  loss  of  some  solutions) 

!sin{2cQ2^k)  —  sin(coi')  =  0 

sm{ci2^k)  —  sin(ciA:)  =  0 

which  can  be  solved  simultaneously  under  the  additional  conditions  \\h\\  =  1  and  q  =  1  to  yield 

I  Co  =  25^^^ 

[  Cl  =  TT 

With  p  =  0  and  Co  =  tt,  the  wavelet  obtained  is  the  one  defined  in  equation  (1);  p  =  2  and  cq  =  =7r 
yields  the  wavelet  of  equation  (2).  For  p  =  1  and  cq  =  |7r,  a  new  bandlimited  orthonormal  mother 
wavelet  defined  by 

1  |u^*l  G  [^TT,  tt)  U  [27r,  ^tt) 

0  otherwise 


is  obtained. 


4.1.4  Non-Flat  n-Slot  Orthonormal  Wavelets 

Denoting  p  =  equation  (4)  becomes 

[  du=  [  du  =  2xt/;(fc)  =  0.  (6) 

it 

for  all  non-zero  integers  k.  This  formulation  leads  to  the  following: 

Theorem  1  Let  i;  be  a  real-valued  and  non-negative  function  with  orthogonal  support.  Further 
suppose  dcj  =  27r  and  i^{k)  =  0  for  all  non-zero  integers  k.  Then  with  \h\  =  and  the 

phase  of  h  chosen  arbitrarily  (except  for  conjugate  symmetry),  h  is  an  orthonormal  mother  wavelet. 
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In' view  of  this  result,  construction  of  the  desired  non-flat  wavelets  hinges  on  finding  functions 
^  having  the  properties  specified.  Such  functions  may  be  constructed  as  follows. 

Denote  sinc(0  =  and  observe  that  sinc(A:)  =  0  for  all  non-zero  integers  k.  Its  Fourier 

transform  has  values 

^  J  1  |w|  <  JT 

[  0  otherwise 

Suppose  5  =  U  •••U  [cn-.r2P”, Cn2P'*]  is  an  n-slot  orthogonal  support  in  which  each 

slot  has  measu;*e  larger  than  2x.  Let  /  be  any  non-trivial  function  supported  in  [co2^^  +7r,  Ci2^^  — 
tt]  U  •  •  -  U  [cn-i2P”  +  7r,Cn2P”  —  tt]  (figure  1)  and  w^hose  convolution  with  sine  is  non-negative.  Then 
=  f[t)smc{t)  is  zero  for  aU  non-zero  integer  values  of  =  (sine  ♦  /]  is  non-negative,  and  the 
support  of  ^  is  orthogonal.  Hence,  an  orthonormal  mother  wavelet  h  can  be  obtained  by  letting 


\h\^  = 


/jI  ^*(^)  duj 


and  setting  the  phase  of  h  arbitrarily. 


(a) 


k 


(b) 


c  2 
k-} 


(D 


Figure  1:  (a)  Each  interval  in  the  original  n-band  orthogonal  support  5  has  length  larger  than  27r. 
(b)  The  support  of  /  is  formed  by  shrinking  each  interval  in  5  by  tt  on  its  left  and  right  ends. 


Note  that  h  may  be  a  dilated  version  of  another  orthonormal  mother  wavelet.  If  after  rescaling, 
there  exists  some  N  €  N,  such  that  h{2'“^t)  satisfies  equation  (4),  i.e., 


for  all  nonzero  integers  A:,  and 


Jr 


dcj  =  0 


for  some  nonzero  integer  k,  then  h'  defined  by  h'{t)  =  ^t)  is  also  an  orthonormal  mother 

wavelet. 


i 


The  first  part  of  the  following  result  is  a  consequence  of  the  above  construction;  the  converse 
portion  is  proven  in  [10]. 

Theorem  2  For  any  n-slot  orthogonal  support  Soy  there  exist  orthonormal  mother  wavelets  having 
support  So-  Conversely f  any  n-slot  orthonormal  mother  wavelet  can  be  constructed  in  the  way 
discussed  above. 


Example  3  Let  cq  =  7r,pi  =  1,  ci  =  77r/4,  and  C2  =  2.  Then  the  support  of  ft  is  [7r,77r/4]u[7;r,87r]. 
Expand  this  support  by  a  factor  of  4  to  [47r,77r]u  [28T,327r]  to  make  the  length  of  each  interval 
larger  than  27r.  Choose  the  support  of  /  to  be  [57r,67r]U  [297r,31?r]  and  generate  /  subject  to  the 
condition  that  /  *  sine  must  be  real  and  nonnegative;  here  define  /  by 


Then 


m  = 


1  5?:  <  \u\  <  Qn  and  297r  <  |u;|  <  ZIk 
0  elsewhere 


=  [/  *  slncjp) 


and  lb  is  as  shown  in  figure  2.  Taking  the  square  root  of  y;  as  the  absolute  value  of  ft  and  imposing 
zero  phase,  the  orthonormal  mother  wavelet  with  Fourier  transform  depicted  in  figure  3  is  obtained. 


2n 


Figure  2:  Convolution  of  /  with  the  Fourier  transform  of  a  sine  function  ensures  =  0  for  all 
non-zero  integers  k. 


4.2  Orthonormal  Wavelet  Symbols  with  Frequency  Overlap 

If  bandlimited  wavelets  are  to  be  used  as  orthogonal  symbols  in  a  low-probability  of  exploitation 
(LPE)  communication  system,  it  may  be  desirable  to  have  the  bands  overlap  to  avoid  the  possibility 
of  an  unintended  receiver  separating  the  channels  using  a  relatively  simple  bank  of  bandpass  filters. 
Following  a  technique  developed  by  Suter  and  Oxley  [46],  bandlimited  wavelet  symbols  have  been 
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Figure  3:  An  orthogonal  mother  wavelet  h  is  obtained  as  the  inverse  Fourier  transform  of  the  square 
root  of 

developed  which  (i)  generate  orthonormal  sets  and  (ii)  whose  replicates  at  different  scales  overlap 
in  the  frequency  domain.  The  approach  is  summarized  below. 

The  construction  given  in  [46]  assumes  that  the  real  Line  has  been  partitioned  into  disjoint 
intervals  Ij,  j  G  Z.  If  €  Z}  is  an  orthonormal  basis  for  then  {/j.itjj, A:  €  Z}  is  an 

orthonormal  basis  for  X^(iR).  A  new  orthonormal  basis  is  constructed  by  extension  of  the  functions 
fj^k  off  the  interval  /j  in  a  specific  way,  resulting  in  a  basis  in  which  the  elements  •‘overlap.*’ 

This  construction  is  modified  to  the  wavelet  case  to  construct  ‘'overlapped**  bandlimited  or- 
thoriormal  wavelet  symbols  as  follows.  Let  h  be  a  bandlimited  w'avelet  having  the  properties: 

1.  The  Fourier  transform  h  of  h  \s  real-valued  and  even.  The  steps  in  the  construction  described 
below*  apply  to  the  portion  of  h  on  the  positive  frequency  axis  and  are  assumed  to  also  be 
applied  to  the  portion  of  h  on  the  negative  frequency  axis  to  maintain  even  symmetry. 

2.  The  support  of  h  is  orthogonal  in  the  sense  introduced  above;  i.e.,  the  portion  on  the  positive 
frequency  axis  is  of  the  form 

(2‘^‘7r,2^‘ai]  U  •  •  -U  [2‘'''a„_i,2^"27r] 

where  n  >  1,  tt  <  aj  <  •  •  •  <  <  2",  and  d\,...,dn  are  integers.  The  construction  given 

here  also  assumes  di  ^  2rfn- 

For  integers  j  and  k,  define 

hj,k{t)  =  2-U  (7) 

and  denote  by  hj^k  the  Fourier  transform  of 

Suppose  {hj^k\jjk  G  Z}  is  an  orthonormal  basis  of  Z^(IR)  (some  examples  of  such  bases  are 
constructed  above  and  others  are  described  in  [6]).  With  e  <  min{(ai  -  7r)/2,(27r  - 
define  an  extension  h  of  ft  by 


1.  Constructing  the  odd  extension  of  /i  about  the  point  2^‘7r  into  the  interval  [2^‘(7r  -  €),2^^n]; 

2.  Constructing  the  even  extension  of  A  about  the  point  2^^27r  into  the  interval  [2^'*27r,2‘^”(27r-{-c)]; 
and 

3.  Repeating  these  steps  on  the  negative  frequency  axis  to  preserve  even  symmetry. 

An  example  is  shown  in  figure  4  using  Mallat’s  wavelet  defined  in  equation  (1). 

Denote  the  extended  wavelet  just  constructed  by  h  and  let  w  denote  an  frequency  domain 
window  function  with  support  identical  to  that  of  h  and  with  the  amplitude  normalization  properties 
described  in  [46].  An  example  of  w  is  depicted  in  figure  4.  With  u  defined  by 

u(u;)  =  ii;(u;)^(u;)  (8) 

the  set  2}  is  shown  to  be  an  orthonormal  wavelet  basis  of  Z'(IR).  Also,  by  construe* 

tion,  the  dilated  replicates  of  u  have  non-trivial  overlap.  Figure  5  shows  a  particular  overlapped 
bandlimited  orthonormal  wavelet  symbol  u{^)  constructed  in  the  above  example. 


•>  0  i  10-10  -S  0  S  to 


Figure  4:  Initial  frequency  extension  h{^')  of  the  wavelet  h  (left)  and  the  window  function  w{^') 
applied  to  it  (right). 

4.3  Countermeasures 

A  detector  for  cyclostationary  signals  was  developed  during  this  effort.  This  detector  (depicted 
in  figure  6)  uses  magnitude-squared  coherence  (MSC)  estimation  as  a  measure  of  the  spectral  cor¬ 
relation  in  a  signal.  The  use  of  MSC  estimation  allows  the  statistical  behavior  of  the  detector 
in  a  noise-only  environment  to  determined  analytically  and  threshold  values  corresponding  to  de¬ 
sired  false  alarm  probabilities  to  be  computed  [15,  16],  These  capabilities  represent  a  substantial 
improvement  over  similar  cyclostationary  detectors  described  in  the  research  literature. 
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Figure  5:  Example  of  an  overlapped  bandlimited  orthonormal  wavelet  symbol  u(a;). 


Figure  6:  Discrete-time  version  of  cyclostationary  feature  detector. 


Receiver  operating  characteristic  (ROC)  curves  for  this  detector  against  various  types  of  cy¬ 
clostationary  communication  signals  have  been  determined  by  (i)  using  the  analytical  results  to 
set  detection  a  threshold  for  a  particular  false  alarm  probability,  and  then  (ii)  running  computer 
simulations  with  signal  present  to  estimate  the  corresponding  detection  probability.  Performance 
results  were  obtained  by  applying  the  detector  against  standard  communication  signals  (e.g.,  polar 
BPSK)  and  against  some  LPE-type  communication  signals,  including  direct  sequence  spread  spec¬ 
trum,  transmitted  reference  spread  spectrum,  and  SDMA  signals.  The  case  of  SDMA  signals  was 
of  particular  interest  in  this  project,  and  is  discussed  further  below. 

Consider  an  M-channel  SDMA  signal  r  described  by 

A/— 1  CO 

KO  =  E  E  ^rnMt)  (9) 

m=0  k—~-oo 

where  bm,k  represents  the  bit  on  the  channel  and  is  time-shifted  replicate  of  the 
wavelet  symbol  ip  at  dilation  level  m.  The  spectral  correlation  density  (SCD)  of  such  a  SDMA  signal 
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was  found  to  be 


where  i;  is  the  Fourier  transform  of  i;,  f  is  the  spectral  frequency,  a  is  the  cyclic  frequency,  and  To 
is  the  fundamentaJ  bit  period  (i.e.,  the  time  shift  required  for  orthonormality  of  the  wavelet  symbols 
at  the  lowest  scale  of  dilation).  This  expression  shows  that  the  SCD  of  a  SDMA  signals  will  be 
non-zero  only  at  cyclic  frequencies  that  are  integer  multiples  of  I/Tq.  This  behavior,  which  has  also 
been  verified  empirically,  indicates  that  SDMA  signals  exhibit  cyclostationary  characteristics  that 
can  possibly  be  exploited  for  the  purpose  of  signal  detection. 

The  cyclostationary  feature  detector  was  applied  to  detect  the  presence  of  a  SDMA  signal  in 
noise.  A  six-channel  SDMA  signal  was  used  as  the  input  signal  for  the  single-cycle  detector.  The 
wavelet  symbol  used  in  the  SDMA  signal  was  the  Daubechies-4  wavelet.  An  initial  analysis  was 
performed  to  determine  at  which  values  of  a  the  spectral  correlation  density  would  be  nonzero.  It 
was  found  that  the  SCD  had  a  strong  component  at  a  =  2^/To  and  /  =  0.  The  single-cycle  detector 
was  set  to  operate  at  these  values  of  /  and  a.  The  number  N  of  independent  samples  used  in  the 
correlations  performed  by  the  single-cycle  detector  w'as  set  to  32.  The  bandwddth  of  the  single-sided 
bandpass  filters  was  set  to  1,28/To.  Simulations  were  run  at  three  different  SNR's;  OdB.  3dB  and 
9dB.  One  thousand  outputs  were  used  in  obtaining  the  empirical  cumulative  distribution  for  each 
SNR  setting.  The  ROC  curves  obtained  by  this  procedure  are  showm  in  figure  7.  These  results 
suggest  that  SDM.A.  signals  are  vulnerable,  at  least  to  some  extent,  to  detection  by  cyclostationary 
methods. 

Another  detector  for  polycyclic  signals  based  on  generalized  coherence  estimation  [7]  has  been 
defined  and  its  application  against  various  types  of  communications  waveforms,  including  SDMA, 
is  currently  under  investigation. 

4.4  Generalized  Wavelet  Transforms 

Work  in  connection  with  this  project  has  led  to  consideration  of  how  the  continuous  wavelet 
transform  can  be  generalized  to  higher  dimensions.  Some  mathematical  results  arising  from  this 
work  and  a  connection  with  image  analysis  are  summarized  in  this  section. 

Several  multidimensional  generalizations  of  the  one-dimensional  wavelet  transform  are  currently. 
In  the  most  widely  known  (and  obvious)  generalization,  one-dimensional  transforms  are  applied 
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Figure  7:  ROC  curves  for  the  single-cycle  detector  with  a  six-channel  SDMA  signal  at  its  input, 

separately  in  each  orthogonal  direction.  In  what  follows,  this  will  be  referred  to  as  the  rectilinear 
wavelet  transform;  it  is  the  method  used  by  most  wavelet-based  image  compression  techniques. 
Another  wavelet  generalization  is  the  two  dimensional  transform  by  Murenzi[37],  literature  on  which 
can  be  more  easily  found  in  the  paper  by  Antoine  et  aL  [1].  Here,  it  will  be  called  the  circular 
wavelet  transform,  since  it  generalizes  the  dilations  of  the  wavelet  transform  to  the  set  of  dilations 
and  rotations  on  15.^.  Neither  of  these  generalizations  encompasses  the  other.  Thus,  a  question  of 
interest  is  w'hether  there  is  a  stronger  generalization  that  encompasses  both.  The  remainder  of  this 
section  describes  a  generalized  wavelet  transform  that  includes  as  special  cases  the  one-dimensional, 
rectilinear,  and  circular  wavelet  transforms.  Among  the  multidimensional  wavelet  transforms  that 
arise  from  this  approach  are  several  that  appear  to  be  unknown  in  engineering  applications.  The 
circular  wavelet  transform  is  shown  to  arise  as  a  result  of  the  relationship  between  and  C  in  the 
context  of  this  generalization  and  is  further  shown  to  be  strongly  related  to  the  cortex  transform 
of  Watson  [54].  An  implementation  of  the  cortex  transform  from  the  circular  wavelet  transform  is 
developed  as  an  example  of  this  relationship. 

4.4.1  Theoretical  development 

Recall  first  the  essentials  of  one-dimensional  wavelet  analysis,  details  of  which  can  be  found  in 
[14].  With  a,6,x  €  E,  /i  €  i^(E),  and  h  the  Fourier  transform  of  h,  define 
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If  Ch  is  finite,  h  is  admissible.  In  this  case,  the  reconstruction  theorem  states  that  any  /  G  i^(IR) 
can  be  reconstructed  from  its  wavelet  coefficients  ^/,  by  the  formula 


The  rectilinear  wavelet  transform  generalizes  the  one-dimensional  transform  to  a  separable  mul¬ 
tidimensional  transform  by  letting  a,6,x  G  i.e,, 


with  the  following  alterations  in  the  definitions  and  theorem: 


(.t)  = 


/ai  •  •  * 


C,  =  / 

pri  •  -Xn\ 

f=-^[  [  - 

Ch  af \  f 

This  generalization  is  equivalent  to  taking  the  one-dimensional  wavelet  transform  successively  in  each 
orthogonal  direction  of  the  multidimensional  space.  For  n  =  1,  the  rectilinear  wavelet  transform 
reduces  to  the  one- dimensional  wavelet  transform. 

An  alternative  generalization  of  the  one-dimensional  wavelet  transform  to  is  the  circular 
wavelet  transform  [1,  37],  for  which  a  G  E  is  replaced  with  the  pair  (a,  5)  G  E  x  [0,2?:)  and  6  G  E  is 
replaced  by  6  =  [6i,  62]"^  €  Let  tq  denote  the  matrix  which  performs  rotation  in  by  Q  degrees. 
Then  the  definitions  and  theorem  for  circular  wavelet  analysis  are 

^^(a, 6)^1  =  i  /i  -  bfj 

C,  =  / 

Jn-^-  llilr 

f=^l  a(/6d0c^a. 

Ch  Jo  Jw.'^  '  '  ' 
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Clearly,  the  circular  wav^et  transform  has  a  different  character  than  the  rectilinear  wavelet 
transform  on  In  particular,  it  allows  non-separable  wavelets  in  two  dimensions  and  replaces 
independent  dilation  in  each  dimension  by  rotation  and  dilation  simultaneously  in  both  dimensions. 
Despite  their  distinct  characters,  both  the  rectilinear  and  circular  wavelet  transforms  arise  as  special 
cases  of  a  more  general  multidimensional  wavelet  transform,  the  formulation  of  w'hich  is  based  on 
ring  multiplications  on  E^. 

Consider  E^  as  a  vector  space  over  the  field  E,  If  a  ring  multiplication  is  introduced  to 
so  that  E”  becomes  an  algebra  under  the  multiplication  then  there  is,  for  each  a  6  E"^,  a  left  regular 
representation  Xa  :  E^  E”  defined  by  LaX  =  n  o  or.  There  is  also  an  associated  linear  operator 
5r  *  defined  by  5xU  =  La^x.  Suppose  the  ring  multiplication  has  the  property  that  the  set 

of  elements  with  ring  inverses  are  dense  in  E^.  Then,  with  h  €  and  a,  6,  a:  6  E^,  the  following 

definitions  and  theorem  form  the  mathematical  foundation  of  a  generalized  multidimensional  wavelet 
analysis: 


\hix)r- 

I  det  5x1 


As  examples  of  this,  suppose  the  ring  multiplication  on  is  defined  by 


ai 

■  bi  ' 

i 

0 

: 

= 

: 

-  . 

.  b.  . 

then  La  and  Sa  are  the  diagonal  matrices  defined  by  La  =  Sa  —  diag(ai,  02, . . . ,  an).  This  makes  the 
generalized  definitions  and  theorem  equivalent  to  the  definitions  and  theorem  of  rectilinear  wavelet 
analysis.  Since  rectilinear  wavelet  analysis  with  n  =  1  is  the  same  as  one-dimensional  wavelet 
analysis,  both  one- dimensional  wavelet  analysis  and  rectilinear  wavelet  analysis  arise  as  special 
cases  of  this  formulation. 

■*> 

Similarly,  circular  wavelet  analysis  comes  from  the  ring  multiplication  ""0”  on  E^  defined  by 


ax 

0 

’  bi  ' 

aibi  —  a2i>2 

a2 

axb2  ”1“ 

This  ring  multiplication  is  same  as  the  multiplication  of  the  complex  numbers.  In  this  case 


La 


ai  —ci'i 

(1‘2  CLl 


H’-Za’ 
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and  the  generalized  definitions  and  reconstruction  theorem  become  equivalent  to  the  definitions  and 
reconstruction  theorem  of  circular  wavelet  analysis. 

4.4.2  Connection  to  the  cortex  transform 

Circular  wavelet  analysis  has  fundamental  variables  of  scale  (which  corresponds  to  logarithmic 
frequency),  rotation,  and  two-dimensional  positional  shift.  These  correspond  to  the  basic  organiza¬ 
tion  of  Watson’s  cortex  transform[54],  which  was  devised  as  a  means  to  process  images  in  a  similar 
way  to  the  processing  in  the  human  visual  cortex.  Hence  there  is  a  relationship  between  the  circular 
wavelet  transform,  the  cortex  transform,  and  image  processing  in  the  human  visual  cortex.  It  is 
shown  here  that  this  relationship  can  be  strengthened  by  actually  deriving  the  cortex  transform 
from  the  circular  wavelet  transform. 

The  cortex  transform  was  implemented  by  Watson  by  partitioning  the  frequency  domain  into 
separate  angularly  oriented  regions  at  logarithmically  spaced  frequency  steps,  bandlimiting  the 
image  to  each  region,  and  sampling  the  bandlimited  pieces  as  efficiently  as  possible  using  available 
techniques  from  sampling  theory. 

Following  a  similar  procedure,  the  frequency  domain  can  be  divided  according  to  the  figure 
below,  where  zero  frequency  is  in  the  center: 


This  division  corresponds  to  a  discrete  sampling  of  the  frequency  and  rotation  variables  of  the 
circular  wavelet  transform.  In  order  to  bandlimit  the  image  to  a  given  region,  the  ‘‘mother  wavelet” 
h  is  chosen  so  that  its  Fourier  transform  h  is  the  characteristic  function  a  6m,  which  consists  of  two 
opposing  frequency  regions  as  shown  below: 


Once  the  image  is  bandlimited  to  a  bin,  positional  samples  are  used  to  represent  the  bandlimited 
function.  These  positional  samples  produce  frequency  replications.  Aliasing  is  avoided  by  preventing 
these  frequency  replications  from  overlapping.  An  example  of  an  optimal  sampling  for  a  bin  is  shown 
below. 


This  frequency  replication  pattern  is  according  to  theory  the  optimally  efficient  replication  that 
avoids  aliasing,  but  in  practice  aliasing  is  not  avoided  due  to  filter  imperfections.  Thus,  it  is  often 
desirable  to  introduce  a  controlled  amount  of  inefficiency  in  order  to  reduce  aliasing.  In  this  case, 
the  frequency  replication  pattern  looks  more  like  this: 

The  efficiency  of  this  sampling  is  controlled,  so  it  is  easy  to  make  it  more  efficient  than  the  fixed 
efficiency  sampling  proposed  by  Watson. 

Preliminary  results  show  that  the  transform  adequately  represents  an  image  and  reconstructs 
it  from  its  representation  (figure  8).  In  order  to  obtain  these  results,  special  consideration  must 
be  paid  to  preserving  frequencies  near  zero,  since  the  transform  has  a  singularity  there.  It  is  also 
crucial,  that  the  image  be  of  relatively  large  size  (at  least  256  x  256  pixels),  since  the  processing 
assumes  an  image  is  continuous,  and  the  violation  of  this  assumption  by  a  coarsely  sampled  image 
introduces  significant  artifacts. 
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Personnel  contributing  to  this  research  effort  were: 


•  D.  Cochran,  Principal  Investigator 

•  S.  Enserink,  Research  Assistant 


Figure  8:  The  Leena  image  (left)  and  the  result  after  it  has  been  decomposed  and  reconstructed  by 
the  transform  (right). 

•  S.  Han,  research  Assistant 
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•  C.  Wei,  Research  Assistant 
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Multidimensional  Wavelet  Transforms  and  an  Application  to  Image  Processing. 
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[6a-8]  S.  Han  and  D.  Cochran,  “Orthonormal  Bases  of  Bandlimited  Wavelets  with  Frequency 
Overlap,”  Proceedings  of  the  IEEE  International  Conference  on  Acoustics,  Speech,  and  Signal  Pro¬ 
cessing,  vol.  3,  pp.  1541-1543,  Detroit,  May  1995. 

[6a-9]  S.  Enserink  and  D.  Cochran,  “On  Detection  of  Cyclostationary  Signals,”  Proceedings  of  the 
IEEE  International  Conference  on  Acoustics,  Speech,  and  Signal  Processing,  vol,  3,  pp.  2004-2007, 
Detroit,  May  1995. 

[6a-10]  D.  Cochran,  H.  Gish,  and  D.  Sinno,  “A  Geometric  Approach  to  Multiple-Channel  Signal 
Detection,”  IEEE  Transactions  on  Signal  Processing,  vol.  SP-43(9),  September  1995. 

6b.  Papers  accepted  for  publication 

[6b-l]  D.  Cochran  and  J.J.  Clark,  “Time- Warped  Bandlimited  Functions:  Sampling,  Bandlimit- 
edness,  and  Uniqueness  of  Representation,”  To  appear  in:  SIAM  Journal  on  Applied  Mathematics. 

[6b- 2]  A.  Clausen  and  D.  Cochran,  “An  Invariance  Property  of  the  Generalized  Coherence  Esti¬ 
mate,”  To  appear  in:  IEEE  Transactions  on  Signal  Processing. 


19 


6c.  Recent  submissions 

[6c- 1]  A.  Clausen  and  D.  Cochran,  “Generalized  Coherence:  Invariance  of  Estimate  to  Distri¬ 
bution  of  One  Channel,”  Submitted  to  the  30th  Asilomar  Conference  on  Signals,  Systems,  and 
Computers,  November  1996. 

[6c-2]  L.L.  Scharf  and  D.  Cochran,  “On  the  Distribution  of  Linear  Prediction  Error,”  Submitted 
to  the  30th  Asilomar  Conference  on  Signals,  Systems,  and  Computers,  November  1996. 

[6c-3]  J.Q.  Trelewicz  and  D.  Cochran,  “A  Sampling  Approach  to  Bandlimited  Instantaneously 
Companded  Signals,”  Submitted  to  the  30th  Asilomar  Conference  on  Signals,  Systems,  and  Com¬ 
puters,  November  1996. 

6d.  Other  related  publications 

[6d-l]  D.  Cochran,  “Review  of  Cyclostationarity  in  Communications  and  Signal  Processing,'" 
(Invited  review)  in  IEEE  Communications  Magazine^  vol.  32,  no.  4,  April  1994.  Also  in  Proceedings 
of  the  IEEE,  October  1994. 

7  Interactions/Transitions 

7a,  Meetings,  Conferences,  and  Seminars 

Numerous  university  and  industrial  colloquium  presentations  have  been  presented  on  the  research 
performed  in  this  project.  These  include: 

•  D.  Cochran  presented  a  talk  entitled  “Wavelet-Based  Multiple- Access  Spread  Spectrum”  in  the 
Motorola  Government  Systems  and  Technologies  Group  Signal  Processing  Series  [16  attendees] 
in  March  1993. 

•  D.  Cochran  presented  a  talk  entitled  “Detection  of  Cyclostationary  and  Polycyclostationary 
Signals”  at  the  DSP  group  meeting  at  MIT  [17  attendees]  in  October  1994.  During  the  same 
visit,  he  met  with  A.  VVillsky,  C.  Karl,  and  R.  Learned  to  exchange  ideas  on  wavelet  and 
wavelet-packet  methods  in  spread  spectrum. 

•  D.  Cochran  presented  a  talk  entitled  “Multidimensional  Wavelet  Generalizations  and  Watson’s 
Cortex  Transform”  in  the  Vision  Seminar  at  Harvard  [25  attendees]  in  September  1994. 

•  D.  Cochran  presented  a  talk  entitled  “W^avelet  Applications  in  Communications”  in  a  seminar 
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sponsored  jointly  by  the  DSP  and  Communications  groups  at  Georgia  Tech  [45  attendees]  in 
October  1994. 

In  addition,  D.  Cochran  participated  in  the  following  conference  activities  related  to  the  topic 
of  this  project: 

•  He  organized  and  chaired  an  invited  session  entitled  “Applications  of  Wavelets  and  Chaotic 
Signals  in  Communications”  at  the  Asilomar  Conference  on  Signals,  Systems,  and  Comput¬ 
ers  in  November  1993.  Participants  included  R.  Orr  (Atlantic  Aerospace),  S.  Isabelle  (MIT), 
M.  Motamed  and  A.  Zakhor  (UC  Berkeley),  C.  Wei  (Arizona  State  University),  and  J.S.  Gold¬ 
stein  (USAF  and  Georgia  Tech). 

•  He  was  a  keynote  speaker  (with  J.J.  Benedetto,  M.V.  Wickerhauser,  and  W.  Sweldens)  at  a 
wavelet  applications  workshop  in  Melbourne,  Australia  in  February  1995.  The  topic  of  this 
presentation  was  wavelet  applications  in  communication  systems. 

•  He  is  scheduled  to  be  plenary  speaker  (with  L.  Cohen)  at  the  ANZIIS-96  Conference  in  Novem* 
ber  1996. 

7b.  Consultive  and  Advisory  Functions 

During  this  project,  D.  Cochran  consulted  informally  with  J.  Stephens  at  Wright  Laboratory 
and  with  Capt.  J.S.  Goldstein  and  Dr.  A.  Lindsey  of  Rome  Laboratory  on  applications  of  wavelets 
and  cyclostationary  signal  processing  in  covert  communications  and  countermeasures.  He  has  also 
been  collaborating  with  the  Centre  for  Signal  and  Information  Processing  Research  (CSSIP)  and 
the  Defence  Science  and  Technology  Organization  (DSTO)  in  Australia  on  wavelet  applications  in 
communications  and  other  topics. 

7c.  Transitions 

None. 


8  New  Discoveries 

No  inventions  or  patent  disclosures  have  resulted  from  this  research  program. 
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9  Honors  /  Awards 

None. 
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1.0  INTRODUCTION 


This  report  describes  efforts  to  develop  a  closed-loop  Direct  Write  Scene  Generation 
(DWSG)  Focal  Plane  Array  Test  Capability  (FPATC)  at  Arnold  Engineering  Development  Cen¬ 
ter  (AEDC)  (Refs,  1-8).  The  Closed-Loop  DWSG  (CLDWSG)  development  efforts  support 
current  and  future  space-sensor  test  efforts.  Recent  strategic  and  tactical  programs  such  as  Alert 
and  Report  Missiles  (ALARM),  Brilliant  Eyes  (BE),  Exo-atmospheric  Kill  Vehicle  (EKV),  Fol¬ 
low-on  Early  Warning  System  (FEWS),  Space-Based  Infrared  System  (SBIRS),  and  Theater 
High- Altitude  Area  Defense  (THAAD)  help  establish  the  need  for  test  capabilities  that  range 
from  chip-level  Focal  Plane  Array  (FPA)  characterization  to  complete  sensor  subsystem  calibra¬ 
tion  and  mission  simulation. 

DWSG  high-speed  photonics  processing  is  required  for  closed-loop  operation  and  simula¬ 
tion  of  dynamic  and  interactive  sensor  commands  that  reposition  the  sensor’s  field  of  view 
within  the  field  of  regard,  and  for  high-fidelity  simulation  of  optical  blurring  and  temporal  effects 
such  as  jitter.  The  real-time  CLDWSG  method  requires  performance  of  a  number  of  operations 
within  the  framing  rate  of  the  focal  plane  array  of  the  sensor.  These  operations  include: 

•  Selection  and  transfer  of  the  scene  background  field  of  regard, 

•  Calculation  and  inclusion  of  target  intensities  and  positions, 

•  Rotation  and  translation  of  the  composite  background  and  target  scenes  to  simulate  interac¬ 
tive  sensor  pointing  commands, 

•  Simulation  of  sensor  and  satellite  jitter, 

•  Simulation  of  optical  blurring  resulting  from  aberrations  and  diffraction  for  broadband 
spectral  radiation, 

•  Integration  of  the  photon  flux  over  each  FPA  pixel  element, 

•  FPA  pixel  responsivity  calibration, 

•  and  compensation  for  acousto-optic  (AO)  cell  optical  modulation  and  other  system  (RF 
generation,  power  amplifiers,  etc.)  inefficiencies. 

AEDC  is  developing  laser-based  CLDWSG  methods  to  simulate  dynamic  sensor  operations 
and  complex  infrared  scenes.  New  photonic  image-synthesis  methods  are  being  developed  to 
employ  image-to-object  Whittaker-Shannon  sampling,  anisoplanatic  optical  convolution  by 
quasi-isoplanatic  spatial  decomposition,  and  high-speed  digital  electronics  for  acousto-optic  mod¬ 
ulation.  Because  of  the  large  terabyte  volume  of  data  to  be  processed,  the  increased  bandwidth 
requirements,  and  the  increased  simulation  fidelity  required  for  DT&E  and  OT&E  of  FPA  sen¬ 
sors,  the  laser-based  DWSG  methodology  is  being  extended  to  accommodate  optical  and 
computational  decomposition  methods  to  better  exploit  highly  and  massively  parallel  real-time 
image  processing  schemes. 
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Section  2.0  describes  the  process  used  to  derive  and  establish  testing  methodologies  and 
requirements  from  sensor  system  missions,  background/target  radiation  and  corresponding  scien¬ 
tific  phenomena,  sensor  parameters,  and  evaluation  and  testing  objectives.  Section  3.0  provides  a 
brief  explanation  and  overview  of  the  DWSG  methodology,  and  Section  4.0  explains  the  method¬ 
ology  used  for  deriving  and  simulating  sensor  optical  properties  and  image  synthesis  methods  to 
include  the  effects  of  aberrations  and  diffraction.  The  CLDWSG  configuration  selection  method¬ 
ology  is  described  in  Section  5.0,  along  with  a  program  developed  to  help  support  test 
configuration  decisions.  The  initial  efforts  to  implement  a  CLDWSG  configuration  using  exist¬ 
ing  FPATC  hardware  is  described  in  Section  6.0.  A  Proof-of-PrincipIe  demonstration  using  the 
existing  FPATC  hardware  is  then  described  in  Section  7.0. 

2.0  SENSOR  MISSION  AND  SOURCE  RADIATION  ISSUES 

Testing  and  evaluation  (T&E)  and  modeling  and  simulation  (M&S)  of  electro-optical  sensor 
systems  to  assess  their  operational  characteristics,  performance  limits,  and  to  support  engage¬ 
ment  scenario  simulations  can  be  accomplished  and  understood  using  a  top-down  approach  for 
establishing  comparative  models  and  corresponding  testing  criteria.  Figure  1  illustrates  a  mis¬ 
sion-driven  approach  to  establish  CLDWSG  requirements  using  a  sensor’s  mission  and 
characteristics  and  related  sensor  T(&E  objectives. 

As  illustrated  in  Fig.  1,  identification  of  the  mission  objectives  and  profiles  and  the  expected 
functional  requirements  are  primary  steps  in  designing,  developing,  and  testing  electro-optical 
systems.  This  is  true  whether  an  optical  system  is  designed  and  developed  for  detection,  discrimi¬ 
nation,  search  and  track,  ranging,  warning,  or  enhanced  imaging,  as  is  done  with  traditional 
FLIR-type  electro-optical  systems.  Depending  on  its  function  and  mission,  a  sensor  will  gener¬ 
ally  have  a  particular  set  of  important  design  and  test  parameters  that  characterize  its  operation, 
function,  and  performance.  Not  all  parameters  will  be  used  universally  for  all  optical  systems. 
Also,  some  parameters  will  vary  in  importance,  depending  on  the  mission  and  the  function  of  the 
electro-optic  sensor  system.  For  example,  the  signal-to-noise  ratio  in  a  detection  or  warning  sys¬ 
tem  is  important  to  increase  the  probability  of  detection  and  decrease  the  probability  of  false 
alarms  with  little  concern  devoted  to  reproducing  the  exact  target-signal  waveforms.  In  contrast, 
imaging  systems  emphasize  the  reproduction  of  scene  imagery  with  little  distortion  of  the  origi¬ 
nal  object  scene.  The  signal-to-noise  ratio  sometimes  receives  very  little  emphasis  in  some 
imaging  systems,  as  compared  to  detection-  or  warning-type  optical  systems. 

Another  important  step  is  identification  of  the  radiation  sources  and  related  phenomenolo- 
gies  that  will  be  exploited  to  accomplish  the  sensor’s  mission,  i.e..  the  target  and  background 
spectral  radiances,  shapes,  sizes,  spatial  exients/distributions,  Weiner  spectrums  (spatial  power 
spectral  densities),  temporal  variations,  etc.  A  target’s  velocity  and  trajectory  are  also  critical  in 
establishing  a  sensor’s  operation  (scan,  step-stare,  staring),  field  of  view,  and  data  processing 
schemes  employed,  as  well  as  the  sensor  testing  methods  chosen  for  accurate  and  adequate  mea¬ 
surement  programs  and  mission  simulations. 
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3.0  DWSG  OPERATIONAL  DESCRIPTION 


The  DWSG  methodology  illustrated  in  Fig.  2  provides  high-speed  photonics  for  visible  and 
infrared  FPA  diagnostics  using  lasers  for  real-time  evaluation  and  scene  simulation.  The  DWSG 
Focal  Plane  Array  Test  Capability  (FPATC)  includes  lasers  operating  at  0.514,  1.06,  5.4,  and/or 
10.6  pm,  and  acousto-optic  (AO)  modulators  for  laser  beam  control  (Ref.  2).  High-speed  RF  elec¬ 
tronics  use  Direct  Digital  Synthesis  (DDS)  to  drive  the  AO  modulators  with  multifrequency 
input  for  muitibeam  output.  DDS  frequency  modulation  provides  extremely  accurate  and  inde¬ 
pendent  deflection  of  the  multibeam  laser  ‘Take”;  DDS  amplitude  modulation  similarly  provides 
accurate  and  independent  multibeam  intensity  control.  Laser  beam  scan  optics  focus  the  laser 
beams  onto  subpixel  regions  for  pixel-to-pixel  photoflux  deposition  with  minimal  stray  radiation 
from  diffraction.  The  laser  beam  rake  is  then  step-scanned  across  the  FPA  during  each  integra¬ 
tion  period,  as  illustrated  in  Figs.  2  and  3. 

Subpixel  irradiation  (illustrated  in  Fig.  4)  is  also  used  to  minimize  DWSG  optically  induced 
cross-talk  signals  and  related  signal  perturbations.  Modular  and  multiple  laser  beam  systems  and 
acousto-optic  photonics  also  help  to  ensure  that  most  of  the  photons  designated  for  specific  sub¬ 
pixel  regions  during  a  given  FPA  frame  are  delivered  to  the  desired  pixel  for  high-fidelity  scene 
simulation  having  similar  signal  and  noise  properties  as  those  with  normal  thermal-source  FPA 
modes  of  operation.  The  DWSG  resolution  (laser  beam  spot  size  and  separation)  is  designed  to 
be  approximately  twice  the  normal  Rayleigh  resolution  criteria  of  most  acousto-optic  modulators. 

DWSG  acousto-optic  noise  sources  (Fig.  5)  are  also  minimized  to  mitigate  undesirable  pho¬ 
tonic  fluctuations.  One  should  note  that  the  equivalence  of  using  lasers  instead  of  thermal  sources 
for  photodetection  evaluation  depends  not  only  on  generating  the  same  mean  number  of  detect¬ 
able  photoevents,  but  also  on  generating  the  same  photonic  noise.  Effective  signal-to-noise  ratios 
and  detection  equivalence  also  depend  on  the  absence  of  signal  artifacts  or  photonic  perturba¬ 
tions  resulting  from  using  lasers  instead  of  thermal  sources.  The  useful  domain  and  relevance  for 
using  lasers  instead  of  thermal  sources  for  photodetecior  evaluation  have  been  investigated  to 
ascertain  differences  and  similarities  regarding  spectral  distributions,  coherence  domains  and 
integration  volumes,  degree/states  of  polarization,  and  aspects  of  photon  packets  (Refs.  3-6). 
Lasers  and  thermal  sources  are  known  to  have  inherently  different  optical  properties  and  differ¬ 
ent  photon  statistics  that  can  be  observed  with  special  photon  counting  and  coherence 
experiments.  However,  in  theoretical  and  experimental  DWSG  domains  of  visible  and  infrared 
wavelengths  and  thermal  sources  with  moderate  temperatures,  lasers  and  thermal  sources  are 
practically  indistinguishable,  exhibiting  near-Poisson  photon  statistics  with  similar  FPA  signals 
and  photonic  fluctuations  (Refs.  4-6)  Experimental  investigations  indicate  that  DWSG  methods 
are  suitable  for  optical  diagnostics  and  evaluation  of  a  variety  of  FPAs  (Ref.  3). 

The  DWSG  laser  beam  is  directed  through  a  2-D  acousto-optic  scanner  to  direct  the  beam 
through  the  scan  optics,  and  focused  onto  the  FPA  subpixel  region.  The  acousto-optic  modulator 
actually  clips  or  truncates  the  laser  beam,  producing  an  intensity  pattern  that  is  neither  Gaussian 
nor  Bessel-function  distributed.  For  a  laser  beam  waist  co  and  circular  aperture  of  radius  a,  a  trun¬ 
cation  ratio  is  defined  as  7  =  a/o)  (Ref.  7).  For  uniform  irradiance  onto  the  AO-cell  aperture,  the 
intensity  pattern  at  the  focal  plane  is  represented  by  the  traditional  Bessel-function  distributed 
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Airy  irradiance  pattern  as  illustrated  in  Fig.  6.  As  the  truncation  ratio  increases,  the  irradiance  pat¬ 
tern  transitions  from  an  Airy  pattern  to  a  Gaussian  beam  pattern,  as  illustrated  in  Figs.  7  and  8. 

There  are  numerous  definitions  for  the  spot  size  of  an  irradiance  pattern,  and  for  a  l/e** 
radius,  the  encircled  energy  for  increasing  truncation  ratios  exhibits  a  maximum  near  7  =  1.4 
(illustrated  in  Fig.  9).  One  could  also  define  and  use  an  84-percent  encircled-energy  spot  size  or 
some  other  value  such  as  90-  or  95-percent  encircled  energy.  Figure  10  illustrates  the  radius  of 
the  laser-beam  spot  as  a  function  of  the  truncation  ratio.  The  current  design  criterion  for  the 
DWSG  is  to  provide  approximately  90  percent  of  the  energy  or  photons  into  the  desired  pixel 
region.  For  T>  1.0,  the  irradiance  patterns  are  very  similar  for  either  a  circular  or  a  square  aper¬ 
ture  which  may  be  desirable  to  further  increase  the  relative  flux  onto  the  desired  pixel  (Ref.  6). 

4.0  BACKGROUNDS,  TARGETS,  AND  SENSOR  SIMULATION  CONCEPTS 

New  CLDWSG  photonic  image  synthesis  methods  are  being  incorporated  to  provide  more 
realistic  optical  simulations  with  either  real  or  synthetic  scenes.  These  photonic  image  synthesis 
methods  can  more  accurately  generate  and  represent  temporal  power  spectral  densities  and  spa¬ 
tial  Wiener  spectrums  for  complex  background  and  target  scene  simulations  that  are  more 
representative  of  scenes  and  convolved  FPA  images  anticipated  in  and  typical  of  real-world  FPA 
sensor  operation. 

4.1  RADIATION  SOURCES 

FPA  sensor  evaluation,  modeling,  and  simulation  ultimately  depend  on  the  mission  and 
objectives  of  the  electro-optical  FPA  sensor  system  and  a  sensor’s  modes  of  operation.  Some 
emphasis  has  been  devoted  to  understanding  the  mission,  background,  and  target  phenomenolo- 
gies;  object-to-image  mappings;  scene  distortions  resulting  from  3-D  to  2-D  radiance  mappings 
from  space  to  sensor  FPA  coordinates;  and  temporal  variations.  The  “reality”  or  “truth”  of  any 
validation  and  verification  effort  has  to  be  cast  in  context  of  the  intended  purpose  and  use  of  the 
sensor  being  evaluated  and  tested. 

Furthermore,  no  absolute  truth  table  or  matrix  exists  to  determine  or  provide  “absolute” 
scenes  for  testing  and  simulation.  In  source  generation,  detection,  and  simulation,  there  are  and 
will  be  natural  photonic  fluctuations  and  spatial/temporal  variations  that  are  necessary  to  provide 
the  inherent  fluctuations  anticipated  in  real-world  engagements.  There  are  many  stochastic  and 
quasi-periodic  fluctuations  that  cannot  be  simulated  absolutely/exactly.  Scene  simulation  and 
detection  validation  should  be  based  on  a  statistical  ensemble  instead  of  any  concept  of  absolute 
scene  “truth”  data.  Any  potential  object  scene  is  then  considered  as  one  sample  from  an  ensem¬ 
ble  of  possible  scenes  that  represent  the  stochastic  statistical  population.  Scene  "truth"  is  then  a 
relative  measure  of  the  ensemble  mean,  variance,  and  higher-order  moments  that  ultimately 
describe  the  population’s  probability  distribution.  Accurate  signal-to-noise  simulations  are  often 
needed  to  judge  the  success  or  failure,  not  of  any  singe  event  or  single  test  case,  but  of  an  ensem¬ 
ble  of  detection  events  that  can  be  described  by  the  normal  laws  of  stochastic  processes  and 
related  detection  criteria  such  as  probability  of  detection,  probability  of  false  alarms,  etc. 
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4.2  SENSOR  VIEWING  CONSIDERATIONS 

To  help  better  understand  potential  methodologies  being  considered  for  closed-loop  DWSG 
operation,  a  number  of  concepts  are  briefly  discussed  to  provide  an  overview  of  some  basic  prin¬ 
ciples  being  considered  and  evaluated.  A  sensor's  geometry  for  above-the-horizon  or  below-the- 
horizon  viewing  is  illustrated  in  Figs.  1 1  and  12.  The  background  and  targets  detected  by  a  sen¬ 
sor  will  generally  depend  on  the  sensor’s  field  of  regard  and  field  of  view,  both  of  which  can  be 
time-dependent,  semi-deterministic  processes  or  stochastic  processes.  Typical  orbital  positions 
will  vary  the  field  of  regard  that  is  viewed  from  space-based  sensor  systems,  and  corresponding 
background-to-FPA  and  target-to-FPA  mappings  have  to  be  considered.  In  addition,  differences 
resulting  from  using  different  wavebands  for  normal  sensor  operations  in  atmospheric  absorption 
bands  or  transmission  windows  of  interest  must  also  be  accommodated.  The  background  and  tar¬ 
get  radiance  levels  reaching  the  sensor  depend  not  only  on  the  emission  or  scattering  properties 
of  the  sources,  but  also  on  the  transmission  of  the  intervening  medium  such  as  the  atmosphere 
and  clouds. 

A  sensor’s  spatial  location  and  relative  orientation  help  to  determine  the  field  of  regard  as 
illustrated  in  Figs.  13  and  14.  A  3-D  geometry  of  space  can  be  viewed  as  being  either  function¬ 
ally  or  optically  mapped  onto  a  sensor’s  FPA,  as  illustrated  in  Figs.  14  and  15.  One  can  either 
think  of  the  background  and  targets  within  the  field  of  regard  as  being  projected  onto  the  FPA,  or 
as  the  FPA  being  projected  into  the  3-D  space,  as  illustrated  in  Fig.  15.  A  set  of  radiance  field-of- 
regard  mappings  for  selected  orbital  parameters  can  then  be  viewed  as  two-dimensional  map¬ 
pings  onto  an  extended  FPA  and  could,  in  principle,  be  precomputed  or  preselected  for  typical 
orbital  parameters.  The  individual  sensor  FPA  image  frames  could  then  be  processed  to  include 
optical  effects  for  time-dependent  pointing  vectors  within  the  field  of  regard. 

4.3  OPTICAL  EFFECTS  AND  IMAGING 

Sensor  imaging  of  quasi-monochromatic,  thermal,  or  blackbody  sources  generally  depends 
on  the  theory  of  partial  coherence.  In  many  cases  this  generalization  can  be  simplified  to  one  of 
the  two  extreme  cases  of  partial  coherence,  either  totally  coherent  or  totally  incoherent,  with  the 
understanding  that  it  is  an  approximation  with  known  errors  acceptable  for  the  intended  purpose 
of  modeling  and  simulation.  Neither  of  these  two  extremes  (coherent  or  incoherent  radiation) 
exists  completely  for  real  sources  which  are  always  partially  coherent  to  some  degree  and  par¬ 
tially  polarized,  also.  Thermal  or  blackbody  sources  are  normally  modeled  using  the  simple, 
incoherent  superposition  assumptions  for  adding  or  integrating  the  statistically  independent 
object  sources  to  obtain  the  integrated  images  with  very  good  approximations  to  the  “real  world” 
for  many  cases.  Each  case  is  normally  evaluated  on  an  individual  basis  to  ensure  that  the  assump¬ 
tions  and  simplifications  are  indeed  representative  of  real-world  engagement  “reality.” 

Imaging  and  photodetection,  as  illustrated  in  Fig.  16,  follow  the  laws  of  optical  diffraction, 
and  optical  components  are  not  perfect.  The  images  of  “point”  sources  are  not  “point”  sources 
because  of  the  inherent  wave-nature  of  photons  and  resulting  diffraction  and  aberration  effects. 
As  a  “point”  source  is  moved  within  the  field  of  view  of  a  sensor,  the  corresponding  image  of  the 
“point”  source  may  vary  noticeably,  depending  on  the  wavelength  and  the  degree  of  aberrations 
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present.  For  practical  purposes,  many  optical  systems  can  be  considered  spatially  invariant  and 
the  image  blur  or  point-spread  function  does  not  vary  as  measured,  and,  for  a  given  wavelength 
range,  the  optical  system  can  be  considered  diffraction  limited.  In  other  cases,  there  may  be  con¬ 
siderable  distortions  and  variations  in  the  point-spread  function  resulting  from  optical  aberrations 
which  need  to  be  taken  into  consideration.  The  optical  components  effectively  map  the  object 
plane  onto  the  image  plane  with  potential  aberrations,  blurring,  distortions,  magnification,  or  min- 
ification  that  transform  the  object  into  the  image.  If  one  can  compute  this  transformation 
function,  then  the  image  can  be  accurately  determined  and  properly  simulated.  Simulations 
should  also  consider  the  nature  of  the  sensor  system  mission,  any  integrated  detector  effects,  and 
the  corresponding  data  processing  algorithms  used.  Verification  and  validation  methods  should 
be  applied  to  demonstrate  and  further  ensure  that  the  evaluation  methods  and  corresponding  sim¬ 
ulations  are  relevant  for  the  intended  purpose. 

There  are  a  number  of  concepts  and  methods  for  determining  a  sensor’s  point-spread  func¬ 
tion  (PSF)  or  its  optical  transfer  function  (OTF)  as  illustrated  in  Fig.  17.  A  sensor’s  image  and 
the  FPA’s  integrated  photoflux  can  then  be  determined  by  basic  integration  methods  or,  for  near 
spatially  invariant  FOV  regions,  by  convolution  methods  or  by  way  of  the  convolution  theorem 
using  OTF  frequency-domain  methods,  as  illustrated  in  Fig.  18.  The  method  chosen  depends  on 
the  desired  fidelity  as  well  as  the  computational  performance  desired.  As  with  any  computational/ 
scientific  model,  these  methods  have  inherent  assumptions  and  simplifications.  These  have  to  be 
understood  to  provide  a  quantitative  method  to  account  for  and  potentially  accommodate  any 
resulting  errors. 

A  geometrical  point  source  is  imaged  into  a  diffraction  pattern.  Only  as  the  geometrical 
image  size  increases  relative  to  the  PSF  dimension  do  we  see  any  structure  or  effects  of  the  shape/ 
size  of  the  “point”  source.  Only  after  the  geometrical  image  dimension  increases  to  on  the  order 
of  magnitude  of  the  Airy  radius  do  we  even  see  the  effects  of  its  shape  and  size.  Even  square 
sources  or  arbitrary  source  shapes  appear  to  be  “point”  sources  when  their  maximum  geometrical 
image  dimension  is  small  compared  to  the  PSF  dimension.  The  spatial  and  structural  features  are 
not  observed  until  the  geometrical  image  dimension  is  on  the  order  of  magnitude  of  the  optical 
system’s  PSF  dimension.  In  the  frequency  domain,  one  says  the  high-frequency  content  is 
stripped  off  and  the  image  does  not  have  sharp  edges  or  discontinuous  spatial  features.  If  a  small 
square’s  diffracted  image  has  the  same  PSF  as  that  produced  by  a  small  circular  source,  they  can 
be  indistinguishable.  Under  some  circumstances  we  can  then  conceivably  simulate  images  with 
sources  that  are  distinctly  different,  yet  yield  the  same  effective  image  as  detected  by  the  FPA. 
Using  the  Whittaker-Shannon  sampling  theory,  we  can  even  use  an  array  of  sources  that  generate 
the  same  irradiance  pattern  if  intensities  are  properly  selected  for  OTF-filiered,  bandwidth-lim¬ 
ited  optical  images. 

One  distinct  feature  of  the  process  of  convolution  is  that  when  the  point-spread  function  is 
“small,”  the  image  and  object  can  be  very  similar.  The  image  is  also  said  to  be  of  high  fidelity 
when  there  are  few  or  no  aberrations  present.  In  the  limit  of  aberrationless,  linear,  shift-invariant, 
delta-function  PSFs,  the  image  will  be  an  “exact”  duplicate  of  the  object  scene  (an  idealism).  As 
the  point-spread  function  increases  in  relative  dimension,  the  image  will  lose  much  of  its  clarity 
and  fidelity;  considerable  differences  between  the  object  and  the  image  can  result.  In  the  limit  of 
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very  large  PSFs,  the  image  will  blur  into  a  uniform  irradiance  pattern  with  near-zero  information 
content  or  modulation  (maximum  entropy),  especially  when  the  PSF’s  dimension  is  large  com¬ 
pared  with  the  largest  geometrical  image  features  present. 

In  the  realm  of  Fourier  optics,  the  object  and  PSF  convolution  can  be  viewed  as  a  filtering 
process.  With  this  perspective,  the  object  is  viewed  in  the  frequency  domain  and  the  OTF  filters 
or  attenuates  the  high-frequency  components  (amplitude  and  phase)  of  the  object  scene  resulting 
in  a  low-frequency  image.  For  incoherent  imaging  and  when  the  geometrical  image’s  spatial  fre¬ 
quency  is  on  the  order  of  magnitude  of  the  reciprocal  of  half  the  PSF  dimension,  the  contrast  or 
modulation  reduces  to  zero  and  the  image  details  or  information  is  lost.  Any  spatial  frequency 
beyond  the  cut-off  frequency  of  the  optical  system  will  be  attenuated  in  the  image  plane.  The 
high-frequency  details  are  essentially  filtered  out  by  the  optical  sensor’s  low-pass,  spatial  fre¬ 
quency  OTF. 

4.4  ANISOPLANATIC  OPTICAL  SYSTEMS 

For  our  applications,  the  scene  and  FPA  are  quantized  into  small  grid-sampled  regions  (as 
illustrated  in  Fig.  19)  to  perform  the  integrated  image  and  detector  photoflux  computations.  The 
degree  or  level  of  sampling  needed  for  a  scene  and  FPA  combination  depends  on  the  desired 
DWSG  simulation  fidelity. 

For  optical  scene  simulations  requiring  real-time  closed-loop  operation,  the  satellite/sensor 
interactively  updates  the  line-of-sight  pointing  vector  relative  to  the  sensor’s  position  in  space 
and  to  the  background  field  of  regard.  To  incorporate  and  support  interactive  CLDWSG  opera¬ 
tion  and  to  accommodate  simulation  of  nonuniform  anisoplanatic  optical  systems,  the  scene  is 
decomposed  into  quasi-isoplanatic  regions  for  narrow  spectral  bands,  as  illustrated  in  Fig.  20. 
For  this  method,  the  scene  is  decomposed  into  segments  that  allow  for  an  application  of  direct 
convolutions  or  OTF  methods  via  the  convolution  theorem  for  regions  of  quasi-isoplanatic 
patches.  The  spatial  decomposition  not  only  provides  for  improved  optical  fidelity  in  simulating 
anisoplanatic  sensor  optics,  but  also  provides  a  method  to  decompose  and  compute  the  scene  seg¬ 
ments  in  parallel.  When  these  methods  are  used,  it  is  important  to  examine  the  effects  of  spatial 
decomposition  and  to  mitigate  edge  effects,  and  diffraction  losses,  while  at  the  same  time  provid¬ 
ing  anisoplanatic  optical  simulation. 

With  the  scene  decomposed  into  quasi-isoplanatic  regions,  one  can  use  either  direct  meth¬ 
ods  or  the  convolution  theorem  for  image  computation  and  synthesis.  For  broadband  radiation, 
effective  PSFs  can  be  used  when  sample  points  have  the  same  relative  spectral  distribution. 
When  there  are  regions  in  the  scene  that  have  considerable  variation  in  spectral  content  (from 
region  to  region),  then  one  could  employ  complete  spectral  decomposition  and  image  synthesis 
for  effective  photon  flux. 

Spatial  FOV  decomposition,  computation,  and  synthesis  of  the  object/image  scene  provides 
a  useful  method  that  facilitates  multiple  PSFs  to  be  used  to  determine  and  directly  simulate  an 
anisoplanatic  optical  system.  Figures  21  and  22  illustrate  that  the  spatial  decomposition  and  com¬ 
putational  methods  may  lead  to  anomalous  computational  artifacts  resulting  from  simulated 
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diffraction/PSFs,  sampling,  and  other -windowing  effects.  Various  methods  to  improve  the  optical 
simulation  and  to  mitigate  computational  artifacts  have  been  investigated  and  illustrated  in  Fig. 
23.  In  these  examples,  extended  domains  help  to  improve  the  simulation  fidelity  and  subse¬ 
quently  improve  experimental  testing  results.  Figure  24  illustrates  the  improvement  obtained 
from  using  extended  computational  domains  to  simulate  anisoplanatic  optical  aberrations  with  a 
reduction  in  computational  artifacts  and  related  eirors. 

For  direct  computation  methods,  the  image  is  integrated  over  each  FPA  pixel  element  to 
determine  the  number  of  photoevents  for  each  frame.  However,  the  integrated  detector  photonux 
can  be  viewed  more  generally  as  a  convolution  of  the  image  irradiance  with  a  detector,  as  illus¬ 
trated  in  Fig.  25.  The  detector  convolution  can  then  be  computed  directly,  and  integrated 
photoflux  values  can  be  obtained  from  the  respective  grid  point  locations  as  desired,  or  arbitrary 
detector  repositioning  can  be  introduced  to  simulate  perturbations  in  the  detector's  position,  row 
offsets,  and  alignment  errors. 

For  many  applications  and  missions,  a  sensor  will  be  detecting  broadband  radiation,  and  one 
needs  to  account  for  the  spectral  distribution  of  radiation  as  well  as  the  spectral  respon.se  of  the 
sensor  system.  The  PSF  and  OTF  can  vary  considerably  due  to  broad  spectral  distributions  (see 
Figs.  26  and  27).  There  are  also  other  blurring  effects  such  as  jitter  which  vary  the  image  consid¬ 
erably  in  real  systems.  When  we  introduce  aberrations  and  polychromatic  or  broadband 
radiation,  the  total  detectable  irradiance  on  the  FPA  will  deviate  from  the  simple  deterministic 
one-case  PSF/OTF  normally  used  for  simulation.  However,  in  practice,  a  composite  PSF  from  a 
polychromatic  point  source  imaged  with  optical  aberrations  and  vibrational  blur  can  be  approxi¬ 
mated  with  a  two-dimensional  Gaussian-distributed  PSF,  as  illustrated  in  Fig.  27. 

To  help  diagno.se,  investigate,  and  simulate  the  etfects  of  potential  real-time  closed-loop 
DWSG  operation  that  includes  optical  diffraction  and  aberration  effects,  a  PC-based  program  has 
been  developed.  An  example  of  the  program's  output  using  the  OTF  method  for  simulation  is 
illustrated  in  Fig.  28. 


5-0  CLDWSG  CONFIGURATION  SELECTION  METHODOLOGY 

The  determination  of  the  optimum  DWSG  configuration  for  testing  a  FPA  is  not  necessarily 
simple  and  straightforward-  FPAs  are  developed  with  various  sizes,  aspect  ratios,  pitches,  speeds, 
and  spectral  bandpa.sses,  and  typically  the  SGTC  AODs  must  be  specifically  arranged  to  meet  the 
needs  of  the  test  article.  In  light  of  frequent  inquiries  from  FPA  vendors  as  to  the  lest  potential  of 
the  SGTC,  an  Excel®  spreadsheet  has  been  created  to  enable  a  DWSG  engineer  to  develcM^  the 
lest  configuration  in  a  systematic  and  timely  manner,  so  that  even  while  the  prospective  is 
watching,  his  questions  can  be  answered.  An  executable  program  written  in  Visual  Basic  has  .use 
been  created  to  provide  a  visual  perspective  of  the  projection  footprints  of  each  individual  AOD 
onto  the  surface  of  the  FPA  under  test. 

Two  operational  constraints  of  the  DWSG  hardware  play  a  role  in  determining  the  lest 
configuration: 
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•  Acousto-optic  deflector  access  time 

•  Acousto-optic  deflector  aperture  (this  with  (1)  and  the  RF  frequency  bandwidth  determine 

the  time-bandwidth  product,  or  resolution  of  the  deflector) 

There  are  several  key  parameters  associated  with  the  test  article  which  also  influence  the 
DWSG  test  configuration: 

•  Number  of  pixels  in  each  axis  of  the  FPA 

•  FPA  pixel  pitch  (this  with  (1)  determines  total  spatial  extent  of  FPA) 

•  FTA  integration  time 

•  Spatial  offsets,  if  any,  between  groups  of  rows  or  columns 

The  interaction  of  these  parameters  is  shown  in  Fig.  29.  The  highly  flexible  modular  nature 
of  the  DWSG  layout  can  accommodate  a  wide  variety  of  test  configurations.  Three  basic 
approaches  to  these  configurations  are  shown  in  Fig.  30. 

The  main  configuration  options  for  the  FPATC  are  shown  in  Fig.  31.  The  TDWSG  configu¬ 
ration  is  similar  to  option  B.  except  that  it  has  only  12  modules  which  cover  75  percent  of  the 
512  X  512  area.  Option  C  illustrates  the  case  when  there  is  a  half-pixel  offset  between  regions  of 
an  FPA.  There  is  also  a  limited  capability  to  perform  continuous  projection  as  in  option  D,  where 
each  AOD  limits  the  scan  to  one  column. 

The  object  in  using  the  Excel  spreadsheet  is  to  be  able  to  respond  very  quickly  to  a  potential 
user's  request  to  test  their  FPA,  and  to  diagnose  difficult  testing  situations  so  that  the  possible 
options  can  be  determined.  It  considers  a  number  of  items: 

•  Projection  extent  of  the  test  article  (pixel  extent  in  each  axis,  presence  of  individual  sub- 
module  FPAs,  if  the  projection  is  in  sections,  or  to  the  complete  FPA,  if  spatial  averaging 
is  used), 

•  Metrics  and  operational  parameters  of  the  FPA  (pixel  pitch,  active  area,  spectral  response, 
mission  bandpass,  fill  factor,  integration  time,  and  type  of  reset), 

•  Projection  hardware  (wavelength,  number  of  AODs  at  operational  wavelengths,  optics, 
aperture,  and  facility),  and 

•  Drive  electronics  (multiple  DRES  components,  multiplexing,  and  data  throughput). 

The  number  of  AODs  and  the  manner  of  their  projection  (coincident  or  adjacent)  through 
each  scan  lens  and  the  number  of  scan  lenses  needed  are  also  considered.  As  the  data  are  input 
and  decisions  are  made  concerning  the  best  means  of  projection,  various  messages  (special  con¬ 
figuration  information  in  blue,  warnings  in  red,  and  suggestions  in  green)  are  produced.  The  user 
can  also  add  notes.  The  Excel  spreadsheet  can  be  printed  out  in  a  convenient  two-page  format 
that  describes  the  optical  and  electronic  configurations  (see  examples  in  Figs.  32  and  33). 
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Once  the  configuration  is  determined  the  derived  optical  dwell  time  must  be  consistent  with 
test  article  compatibility  measurements.  The  optical  power  needed  to  fully  illuminate  the  projec¬ 
tion  must  be  available  with  the  facility  hardware.  Another  Excel  spreadsheet  (illustrated  in  Fig. 
34)  has  been  developed  to  aid  in  this  determination. 

The  configuration  conceptualization  program  illustrates  the  footprints  of  the  AOD  projec¬ 
tions  and  the  relation  to  the  surface  of  the  FPA.  An  example  of  its  output  is  shown  in  Fig.  35. 
Some  additional  concept  layout  may  need  to  be  performed  to  ensure  that  the  derived  configura¬ 
tion  is  physically  possible. 

6.0  INITIAL  CLDVVSG  IMPLEMENTATION  EFFORTS 

This  section  describes  the  cuirent  status  of  recent  AEDC  technology  efforts  tasked  with  the 
implementation  of  real-time,  closed-loop  scene  manipulation  methods  using  AEDC’s  FPATC 
mission  simulation  test  capability.  Figure  36  shows  the  current  FPATC  open-loop  capability/con¬ 
cept  using  the  Direct  Write  Scene  Generation  method  to  test  sensors  at  the  component  or 
subsystem  level.  In  this  configuration,  test  scenarios  are  played  into  the  sensor  under  test  and 
then  captured  at  the  sensor's  output.  Also,  there  is  no  feedback  between  the  radio  frequency  (RF) 
electronics  that  control  the  projection  and  data  acquisition  electronics  that  capture  the  sensor 
output. 

Figure  37  illustrates  a  concept  for  a  closed-loop  FPATC  capability.  The  test  article  is  a  sen¬ 
sor  subsystem  including  FPA  operational  hardware  and  Oighl  navigation  hardware.  Under  closed- 
loop  conditions,  Hight  hardware  tracking  algorithms  output  new  line-of-sight  coordinates  that  are 
accepted  by  FPATC  hardware  and  used  to  create  the  next  sequential  frame  to  be  projected. 
Because  the  DWSG  does  not  project  through  a  sensor's  optics,  real-time  digital  simulation  of  sen¬ 
sor  optical  effects  is  a  critical  part  of  realistic  simulations. 

The  initial  closed-loop  proof-of-principle  (PoPl)  demonstration  focuses  on  using  the 
FPATC  first  to  interactively  manipulate  or  construct  scenario  frames  by  responding  to  changes  in 
line-of-sight  or  pointing-vector  coordinates,  and  second,  to  demonstrate  computational  capabili¬ 
ties  for  simulating  optical  effects  in  the  constructed  image  scene.  The  pointing-vector 
information  can  be  computer  generated,  input  by  a  user  who  is  interactively  viewing  the  sce¬ 
nario,  or  both.  Before  a  PoPl  demonstration  capability  could  be  designed  and  integrated  into  the 
FPATC,  it  was  necessary  to  understand  some  of  the  computational  requirements  for  closed-loop 
operation  in  general,  and  the  PoPl  demonstration  specifically. 

Figure  38  outlines  the  steps  required  to  interactively  compute  a  2-D  infrared- rendered  sce¬ 
nario  frame  based  on  a  pointing-vector  input.  These  steps  include  rotation,  translation,  and 
extraction  of  the  scene  image;  application  of  the  optical  transfer  function  using  a  2-D  FFT  then 
inverse  FFT  approach;  detector  integration  for  greater  than  1:1  image  sampling:  radiometric, 
laser,  and  RF  calibration;  and  projection  buffer  loading  before  outputting  the  frame.  An  equation 
was  then  derived  by  predicting  the  number  of  compute  cycles  required  for  the  whole  process  and 
summing  the  total.  The  equation  is  written  as  a  function  of  lesi  parameters  such  as  oversampling 
ratio,  number  of  FPA  rows  and  columns,  and  point-spread  function  size.  To  be  conservative,  no 
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optimization  was  done  on  the  calculation  of  compute  cycles  per  step.  A  compute  cycle  is 
assumed  to  be  the  time  required  to  perform  an  add.  subtract,  multiply,  divide,  or  memory  move. 

There  are  some  optimizations  that  could  be  done  to  combine  steps  in  the  process  and  speed 
up  the  calculations.  Figure  39  shows  the  equation  that  would  result  if  the  FFT  approach  were 
substituted  with  direct  convolution  using  a  restricted  kernel.  Direct  convolution  is  not  always, 
however,  the  most  optimal  method  to  choose.  Issues  concerning  kernel  size  relative  to  sampling 
ratio  can  cause  fidelity  problems  in  the  simulation  and  create  excessive  computational  loading. 
However,  for  applications  that  permit  the  fidelity  levels  associated  with  a  smaller  point-spread 
function  kernel,  direct  convolution  is  simpler  and  efficient  to  implement.  Figure  40  tabulates 
direct  convolution  versus  FFT  methods  and  indicates  a  crossover  point  between  FFT  and  direct 
convolution  performance  near  a  9  x  9  kernel  size. 

After  considering  worst-case  computational  requirements  for  closed-loop,  a  three-fold 
implementation  strategy  was  adopted:  (1)  capabilities  of  existing  FPATC  hardware  designs  were 
examined  for  PoPl  demonstrations:  (2)  off-the-shelf  hardware  and  software  subsystems  that 
could  be  integrated  into  the  FPATC  were  investigated;  and  (3)  integration  of  off-the-shelf  sub¬ 
systems  and  components  with  in-house  developments  that  could  meet  closed-loop  requirements 
were  investigated. 

Figure  41  illustrates  an  in-house-designed  DWSG  RF  control  module.  The  module  uses  a 
TMS320C30  (C30)  Digital  Signal  Processor- (DSP)  to  control  16  RF  channels  (or  FPA  scan 
lines)  during  a  mission  simulation.  In  open-loop  simulation,  the  system  uses  the  DSP  as  a  data 
pump  to  produce  new  frames  rapidly  during  simulation.  However,  in  closed-loop  simulation,  the 
DSP’s  computational  capabilities  could  be  tapped  to  execute  the  required  ciosed-loop  processing 
steps  on  scene  data  stored  in  each  processor’s  memory.  Thirty-two  C30’s  work  together  in  a  Sin¬ 
gle  Instruction  Multi-Data  (SIMD)  architecture  in  the  FPATC  DWSG  electronics.  Because  the 
SIMD  architecture  is  very  compatible  with  image  processing  and  the  decomposition  methods 
required  for  anisoplanatic  optical  simulation,  the  PoPl  demo  was  implemented  using  existing 
FPATC  designs  and  hardware.  Reference  8  documents  the  hardware  architectures  developed  for 
the  DWSG  control  electronics. 

Figure  42  illustrates  the  anticipated  computational  capabilities  using  all  of  the  DWSG  RF 
control  electronics.  For  example,  the  chart  predicts  that  a  128  x  128  FPA  could  be  operated  at  13 
frames  per  sec  with  full  5x  oversampled  fidelity.  This  frame  rate  is  considered  more  than  ade¬ 
quate  for  proof-of-principle  demonstration.  However,  many  128  x  128  FPAs  operate  at  frame 
rates  in  excess  of  1,000  frames  per  sec.  Such  a  frame  rate  would  require  a  computational  capabi! 
ity  of  better  than  50  GFLOPs  for  high-fidelity  simulation.  These  computational  loads  reinfon^^  ^ 
the  need  for  advanced  hardware  development  beyond  the  PoPl  demonstration  capability  and  con¬ 
tinued  investigation  of  smart  fidelity  trade-offs  and  advanced  algorithm  developments  that  will 
reduce  the  hardware’s  computational  requirement. 

Using  the  FPATC  hardware  and  software  for  the  PoPl  demonstration  is  cost  effective 
because  it  delays,  by  at  least  one  year,  large  material  procurements  for  more  advanced  sub¬ 
systems.  This  approach  recognizes  that  using  an  opportunity  to  wait  for  new  electronic 
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capabilities  and  cost  decreases  while  demonstrating  the  principles  of  closed-loop  operation  is  an 
optimal  path  to  take  in  a  time  of  rapidly  changing  electronics  developments. 

Based  on  results  from  the  PoPl  demonstration,  a  more  detailed  development  plan  will  be 
completed  for  advanced  hardware  and  software  designs.  During  and  after  developing  the  PoPl 
demonstration  with  existing  hardware  designs,  it  is  also  important  to  investigate  other  potential 
vendors  for  signal  processing  hardware  and  software. 

7.0  PROOF-OF-PRINCIPLECLDWSG  DEMONSTRATION  EFFORTS 

The  first  closed-loop  Direct  Write  Scene  Generation  Proof  of  Principle  (PoPl)  Demonstra¬ 
tion  helped  to  evaluate  the  capacity  and  functionality  of  the  current  FPATC  configuration  to 
interactively  manipulate  a  dynamic  scene  projection  by  responding  to  changing  pointing  vectors. 
For  this  PoPl  demonstration,  the  dynamic  scene  consists  of  a  static  background  scene  and  a  mul¬ 
tiple  frame  target  database  overlaid  onto  the  background  scene.  The  PoPl  interactive  interface  is 
mouse-driven,  which  allows  the  user  to  select  the  pointing  vector  coordinates  via  the  mouse  with 
a  real-time  display  depicting  the  current  frame  being  projected. 

The  CLDWSG  PoPl  Demo  consisted  of  three  major  software  components  operating  in  con¬ 
junction  with  the  other.  The  first  software  component  is  the  CLDWSG  Scene  Projection  software 
executing  in  the  DRES.  This  software  component  computes  the  new  frame  information  based 
upon  the  target  file  database  and  the  new  pointing  vector.  This  software  component  is  distributed 
over  multiple  C30  processors  operating  in  parallel.  Each  processor  performs  the  computations  for 
a  portion  of  the  scene  to  be  projected.  A  number  of  steps  are  required  in  the  computation: 

•  Clear  targets  out  of  background  scene  from  previous  frame 

•  Overlay  new  targets  onto  background  scene 

•  Retrieve  new  pointing  vector  from  SCRAMNET  shared  memory 

•  Using  the  new  pointing  vector,  translate,  rotate,  and  extract  the  new  scene  image  using 
nearest-neighbor  method 

•  Perform  Optical  Sensor  PSF  diffraction  effects  using  Direct  Convolution  with  a  3  x  3 
kernel 

•  Integrate  photoHux  over  FPA  pixel  elements 

•  Convert  from  FPA. output  levels  to  RF  dB  levels  using  a  5th-order  polynomial 

•  Convert  from  RF  dB  levels  to  attenuator  commands  using  a  lookup  table 

•  Load  RF  electronics  with  attenuator  commands 

The  second  software  component  in  the  CLDWSG  PoPl  demonstration  is  the  closed-loop 
control  program  which  was  executed  on  a  66- MHz  486  PC.  This  program  interactively  creates 
new  pointing  vectors  consisting  of  2D  axial  and  angular  parameters.  The  angular  values  are  cre¬ 
ated  from  interactively  clicking  the  mouse  button.  A  left-mouse  button  click  increments  delta 
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theta  in  a  positive  angular  direction,  whereas  a  right-mouse  button  click  will  increment  delta 
theta  in  a  negative  angular  direction.  Once  the  pointing  vector  has  been  created,  it  is  transferred 
to  the  CLDWSG  Scene  Projection  software  operating  on  the  ORES  via  SCRAMNET'Shared 
memory.  The  pointing  vector  is  transferred  to  the  CLDWSG  real-time  display  software  operating 
on  a  DIPIX  frame  grabber  interlace. 

The  third  CLDWSG  PoPl  demonstration  software  component  is  the  CLDWSG  real-time 
display.  This  software  component  operates  on  a  C30-based  DIPIX  frame  grabber.  This  display 
assists  the  user  in  creating  the  desired  pointing  vector.  The  entire  background  scene  with  overlaid 
targets  is  displayed  in  real  time  synchronously  with  the  scene  projection.  For  visual  comparison, 
the  translated,  rotated,  and  extracted  image  to  be  projected  is  also  displayed  to  the  side  of  the 
background  scene.  Using  these  two  displays  in  conjunction  with  the  mouse,  the  user  can  select 
the  desired  portion  of  the  background  scene  to  be  projected. 

8.0  CONCLUDING  REMARKS 

In  order  to  provide  more  optimized  optical  simulation  fidelity  and  to  reduce  computational 
burdens,  closed-loop  DWSG  image  synthesis  methods  are  being  developed  which  employ  image- 
to-object  Whittaker-Shannon  sampling,  anisoplanatic  optical  convolution  by  quasi-isoplanatic 
spatial  decomposition,  and  high-speed  digital  electronics  for  acousto-optic  modulation.  Optical 
and  computational  decomposition  will  not  only  provide  high-fidelity  optical  simulation  for 
anisoplanatic  optical  sensors  and  complex  infrared  scenes,  but  will  also  facilitate  high-speed  par¬ 
allel-processing  schemes  for  real-time  CLDWSG  and  sensor  operations. 

The  current  trend  in  FPA  development  is  to  develop  devices  that  operate  at  very  short  inte¬ 
gration  times.  Thus,  to  keep  up  with  FPA  technology,  there  must  be  corresponding  increases  in 
acousto-optic  cell  and  other  scanner  technologies  (Ref  9).  A  Phase  I  SBIR  program  (AF  95-0(35) 
was  funded  by  AEDC  to  develop  new  techniques  for  acousto-optic  deflectors  with  higher  access 
times  (5  p.sec  or  less  for  256  x  256  resolvable  spots;  50-  to  100-MHz  driving  frequencies)  than 
are  currently  possible  with  standard  designs.  One  contract  was  awarded  to  Physical  Optics  Corpo¬ 
ration  based  on  a  multisectional  concept  which  uses  available  acousto-optic  materials.  Prototype 
hardware  for  this  effort  will  be  delivered  to  AEDC.  An  additional  contract  was  awarded  to 
Aurora  Associates,  who  is  currently  procuring  acousto-optic  materials  for  further  study. 

A  computer  testbed  has  been  constructed  and  is  being  used  to  develop  and  evaluate  the  high¬ 
speed  computational  algorithms  required  for  scene  extraction  and  convolutions.  The  proof-of- 
principle  demonstration  provides  assurance  that  the  computation  and  photonics  methods  used  are 
not  only  feasible,  but  the  concept  is  scaleable.  This  technology  effort  helps  to  provide  new  opti¬ 
cal  diagnostics  for  cost-effective  and  systematic  DT&E  and  OT&E  of  large  FPA  sensors  using 
parametric  and  statistical  methods  not  amenable  to  costly  field  or  flight  testing  methods. 
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Driving  Factors  for  Test  and  Evaluation 


Mission  Objectives 


Backgrounds 

Targets 

Atmospheric  Effects 
Sensor  Environment 
Flight/Orbital  Parameters 
Engagement  Scenarios 
Stochastics  &  Ensembles 
Imaging  Systems 
Detection  and  Tracking 
Warning  Systems 
Autonomy/Communications 
Human  in  the  Loop 


Sensor  Properties 


Operational  Modes 
Field  of  View 
Field  of  Regard 
Frame  Rate 
Dynamic  Range 
Resolution 
Spectral  Band(s) 

Focal  Length/Aperture  (F/#) 
Detector  Parameters 
Electronic  Parameters 
Algorithms  &  Processing 
Communications  Links 


Testing*  Requirements 


Functional  Operation 
Performance  Charac. 
Imaging  Properties 
Ranging  Accuracy 
Probability  of  Detection 
Probability  of  False  Alarm 
Tracking  Accuracy 
System  Integration 
Data  Fusion 
Hardware  in  the  Loop 
Interactive  Closed-Loop 
Human  in  the  Loop 


Sensor  Testing  Methodologies  and  Requirements 


Simulation/Testing  Technology  Developments 

Testing  Methodologies 

Models  and  Databases 

Scene  Generation 

•  Performance  &  Efficiency 

•  Component/Subsys./Systems 

•  Ground,  Field,  Flight 

•  Interactive  Suite 

Simulations 

•  Backgrounds 

•  Targets 

•  Atmospheric  Effects 

•  Dynamics  &  Stochastics 

•  Digital  Simulation 

•  Signal  Injection 

•  Direct  Write 

•  Image  Projection 

Closed-Loop  DWSG  Technology  Developments 


•  Sensor/DWSG  HW/SW  Interfaces  &  Environments 

•  Sensor/DWSG  Operation  &  Synchronization 

•  Calibration,  Diagnostics,  Verification,  &  Validation 

•  Bandwidth  Limitations  &  AO  Cell  Performance 

•  Background/Target  Scene  Synthesis  Methodologies 

•  Data  Rates  &  Processor  Throughput 

•  Decomposition  Schemes  &  Parallel  Processing 

•  Closed-Loop  Methodologies 


Figure  1.  Sensor  testing  and  evaluation  diagram  illustrating  closed- 
loop  DWSG  development. 
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Figure  2.  Illustration  of  the  laser-based  Direct-Write  Scene  Generation 
methodology  for  focal-plane-array  diagnostics  and  scene- 
generation  schemes. 
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Tested  Wavelengths: 
Ar  Ion  0.514  pm 
HeNe  0.6328  pm 
Nd:YAG  1.06  pm 
NdiYAG  1.32  pm 
HeNe  3.39  pm 
CO  5.4  pm 

^  COglO.Spm  ) 


Blanks  the 
Beam  During 
Pixel-to-PIxel 
Transitions  ^ 

Acousto-Optic 

Modulator 

Blanker 


Directs  the 
Laser  Beam 
into  the  AOD 
at  the  Correct 
Angle  and 
,  Position  ^ 


Beam 

Steering 

Mirrors 


Converts  RF  ' 
Frequency  into 
Angular  Position 
and 
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Output  Beam  onto 
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Electro-Optic 

Stabilizer 
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Laser  Beam 
Stabilization 


RF 

Electronics 
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Frequencies  Corresponding 
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Directions  with  Beam 
Intensity  Control 


Provides  Cryogenic  , 
Environment  and  ] 
Electrical  Connections 
for  Proper  Operation 
and  Converting 
Photoelectron  Events 
into  Electronic  Signals 
for  Evaluation 


Figure  3.  Laser-based  Direct-Write  Scene  Generation  methods  for 
photo-  characterization  and  evaluation  of  focal-plane- 
array  sensor  systems. 
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J  =  sin  =  (twice  the  Rayleigh  Resolution) 
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total  scan  angle  (p 

smallest  resolvable  angle  Aip 
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2a 


Position  (Airy  Radii) 


Position  (Airy  Radii) 


Figure  4.  DWSG  operation  for  laser-based  photonics  illustrating  the 
resolvable  spots  for  sub-pixel  irradiation  at  twice  the  Rayleigh 
resolution. 
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High-Order  Diffraction  ^ 
Optical  Jitter 
RF  Electrical  Purity 
BRDF/BTDF 
Scene  "Truth" 
Rotational  Grid  Sampling 
Translation  Grid  Sampling 
Digital  A/D  Quantization 
Roundoff/T  runcation 
Computational  . 
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Figure  5.  Potential  noise  sources  for  DWSG  photonics  and  FPA 
operation. 
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Radial  Distance  (Airy  Radii) 

Figure  6.  Normal  Airy  pattern  for 

r=o.o. 


Radial  Distance  (Airy  Radii) 

Figure  7.  FPA  DWSG  laser  beam 
irradiance  pattern  for 
small  truncation  ratios 

r  =  0,  1.0,  and  1.5. 


Figure  8.  FPA  DWSG  laser  beam 
irradiance  pattern  for 
small  truncation  ratios 

T  =  0,  2.0,  2.5,  and  3.5. 


Figure  9.  Percent  Encircled  energy 
for  a  radius  set  at  the 
l/e^  irradiance  level. 


Figure  10.  FPA  laser-beam  irradiance  spot  size  for  various  encircled  energies. 
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Figure  11.  Sensor  view  of  earth  for  above-  and  below-the-horizon 
perspectives  require  quite  different  data  bases  for 
simulation  of  real-world  background  and  target 
engagements  and  related  phenomena  exploited  for  efficient 
detection,  discrimination,  and  tracking. 
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Figure  12.  Multisensor  perspective  and  simulation  will  depend  on 
mission  as  well  as  on  environmental  and  operational 
characteristics. 


AEDC-TR-95 


Figure  13.  Illustration  of  methodology  to  view  backgrounds  and 
targets  for  relative  coordinate  systems  from  a  space-based 
sensor  in  an  earth  orbit. 
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Figure  14.  Illustration  of  a  3D'to-2D  spatial  mapping  of  space-based 
backgrounds  and  targets  onto  the  FPA  for  DWSG  optical 
simulation  efforts. 
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Figure  15.  Concept  illustration  of  the  scene  being  mapped  onto  the  FPA, 
which  could  also  be  usefully  viewed  as  the  FPA  field  of  view 
being  mapped  onto  the  2-D  projection  of  the  background  and 
target  scene. 
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Figure  16.  Illustration  of  the  object  plane  being  transformed  to  the 
image  plane  for  FPA  irradiation  via  the  optical  system  for 
a  simple  background  and  target  map. 
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Figure  17.  Different  concepts  and  methods  for  computing  the  PSF  or 
the  OTF  for  either  direct  convolution  methods  or  OTF 
methods  of  determining  and  simulating  a  sensor’s  imaging 
properties. 


Figure  18.  Optical  convolution  can  be  accomplished  by  either  direct 
methods  or  by  indirect  FFT  methods  using  the  convolution 
theorem  and  the  OTF  of  the  sensor  system  (the  method 
chosen  depends  on  computational/photonic  speed  as  well  as 
fidelity  issues/trades). 
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Figure  19.  Background,  target,  and  FPA  sampling  provide  efficient 
and  flexible  computational  schemes  for  scene  extraction 
and  image  convolution  by  employing  Whittaker-Shannon 
sampling  ideas. 


Figure  20.  Illustration  of  a  spatial  decomposition  method  to  support 
anisoplanatic  optical  simulation  that  could  also  support 
parallel  and  massively  parallel  processing  schemes  for  the 
laser-based  DWSG  photonics. 
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Object  Isoplanatic  PSF  Image  =  ObJect*PSF 


Figure  21.  Illustration  of  a  4x4  optical  decomposition  to  support  parallel 
processing  schemes  and  anisoplanatic  optical  simulation 
without  any  regard  for  diffraction  effects/loss  or  spectral 
leakage. 


4x4  Segmented  Image 
Object  Space  Scene  (although  isoplanatic) 


Figure  22.  Illustration  of  potential  simulation  artifacts  that  can  result 
from  decomposition  methods  that  fail  to  account  for 
diffraction  effects  and  spectral  leakage. 
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Scene  Segment  Hanning  Operator  Convolution  with  quasi-isoplanatic  PSF  Convolved  Image 


Figure  23.  Methodology  for  using  extended  computational  regions  to 
not  only  support  anisoplanatic  decomposition  for  better 
fidelity  and  simulation  but  to  also  support  real-time 
closed-loop  parallel  processing  schemes. 


Object  Quasi-isoplanatic  PSF  Regions  image=Object*PSF 


Figure  24.  Illustration  of  a  decomposition  method  to  arrive  at  the 
composite  field  of  view  for  an  anisoplanatic  optical  system 
being  simulated  in  tbe  laser-based  DWSG  scene  generation 
methodology. 
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Figure  25.  Illustration  of  methodology  to  derive  the  integrated 
photoflux  from  a  double  convolution  process  that  can  be 
accomplished  from  a  more  efficient  system  transfer 
function  operation  and  application  of  the  convolution 
theorem. 
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Figure  26.  Illustration  of  the  eff'ects  of  wavelength  dependent  variation 
of  the  PSF  and  optical  cutoff  frequency  as  seen  in  the 
magnitude  of  the  optical  transfer  functions  (modulation 
transfer  function,  MTF). 
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Figure  27.  Illustration  of  the  concept  of  an  effective  PSF  resulting  from  a 
number  of  contributing  sources  such  as  broadband  radiation, 
optical  diffraction,  deterministic  and  stochastic  jitter  and 
resulting  blur,  or  other  sources. 


Figure  28.  Example  of  the  screen  of  the  closed-loop  simulator  for  algorithm 
development,  diagnostic  efforts,  and  demonstrations. 
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Figure  29.  Dependence  of  DWSG  parameters  determining  jj  test  configuration. 
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Figure  30.  Modular  flexibility  of  the  DWSG  test  configuration. 
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A 
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4 
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B 
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16 
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0* 
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6 
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I^^SUnkOMiirJI 


MUX«1  MUX  #2 
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Figure  31.  FPATC  configuration  options  for  various  FPA  testing  methods. 
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Laser  Irradiance  from  TOWSG  SPA  for  128  x  128  Module  Nd:Yag  diode<pumped 


Amoco  1064EH  350 


100  micron  pixel  size 

1 .06  microns 

80  micron  90  %  encircled  energy  spot 

0.45  mm  beam  diameter 

128  step  by  128  rake  spots 

0.35  W  laser  output 

1  AOD(s)  per  scan  lens 

9.3  mm  AOD  aperture 

1  AOD(s)  per  laser 

25  %  total  diffraction  eff 

100.00  ms  integration  time  (estimated) 

3  dB  bandwidth  rolloff 

98.08  dwell/step  percentage 

1 5  ps  access  time 

1.00E+11  radiometric  background  warm 

std  bandpass  filter 

Conditions  at  exit  of: 

condition 

%  capacity 

transmission 

photon  rate 
(ph/sec) 

beam  dia 
(mm) 

bii 

Laser 

0.35  W 

100.00 

1 

1.87E+18 

0.45 

1.17E+21 

stabilization 

100.00 

1.00 

1.87E+18 

0.45 

1.17E+21 

modulator/blanker 

67.00 

0.9 

1.13E+18 

0.45 

7.08E+20 

correction  wedge 

0.95 

1.07E+18 

0.45 

6.72E+20 

aperture 

0.95 

1.02E+18 

0.45 

6.39E+20 

collimator  +  (diverger) 

0.95 

9.65E+17 

9.3 

1.42E+18 

neutral  density  filter 

1 

9.65E+17 

9.3 

1.42E+18 

1.00 

mi^sm 

9.3 

1.42E+18 

0.96 

9.27E+17 

9.3 

1.36E+18 

AO  deflector  X 

100.00 

0.50 

4.63E+17 

9.3 

6.82E+17 

scan  dwell  time  fraction 

0.70 

1 

3.23E+15 

9.3 

4.75E+15 

bandwidth  rolloff 

0.50 

1.62E+15 

9.3 

2.38E+15 

AO  deflector  Y 

1 

10.00 

0.50 

8.09E+13 

9.3 

1.19E+14 

number  of  beams  in  rake 

128 

0.71 

1 

5.75E+11 

9.3 

8.46E+11 

bandwidth  rolloff 

0.50 

9.3 

4.24E+11 

AO  turning  mirrors 

I  0.97 

2.79E+11 

9,3 

4.11E+11 

Main  beamsplitter/combiner 

1 

2.79E+11 

9.3 

4.11E+11 

Scan  lens  system 

0.9 

2.51  E+11 

9,3 

3.70E+1 1 

Dewar  window 

0.95 

2.39E+11 

1 

3.04E+13 

Laser  line  filter 

0.7 

1.67E+11 

0.5 

8.52E+13 

FPA  (single  pixel) 

96.11 

1 

1.67E+11 

0.1 

2.05E+15 

Throughput  8.96E-08 

Dynamic  Range  for  one  pixel  (above  background  level:  1  E+04  is  essential)  2.05E+04 

Assumes  all  pixels  full  on 

Consider  Additional  Safety  Factor:  -  2  if  two  spectral  bandwidths  are  to  be  illuminated 

Figure  34.  Example  DWSG  laser  irradiance  calculation  spreadsheet. 
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Figure  35.  DWSG  configuration  conceptualization  program  sample  output. 
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Figure  36.  Open-loop  FPATC  concept. 
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Figure  37.  Closed-loop  FPATC  concept. 
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Figure  38.  Computational  steps  required  for  FPATC  closed-loop  simulation. 
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Figure  39,  Closed-loop  computational  steps  using  direct  convolution. 
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Figure  40.  Comparison  of  direct  convolution  and  FFT  computational  methods. 


Figure  41.  DWSG  RF  control  electronics  module. 
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Figure  42.  Theoretical  determination  and  comparison  of  the  FPATC 
RF  control-electronics  capabilities. 
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Platform  Control  /  Status 


The  subject  of  this  report  is  the  detection  range  performance  of  the  HISS  Phase  2  system  against  the  TLX- 1  towed 
height-keeping  target  configured  with  an  APC-6  plumer  and  hot  metal  emitter  (HME).  The  APC-6  has  an  HME  to  produce  IR 
radiation  with  spectral  content  more  characteristic  of  the  body  radiation  from  a  supersonic  (Mach  2+)  cruise  missile 

Originally  the  HME  had  a  diameter  of  6.54  in.  and  a  measured  IR  signature  of  6  W/sr  in-band.  For  the  purpose  of  this 
test,  it  was  replaced  with  a  new  configuration  that  has  a  lO-in.  HME.  The  IR  signature  for  the  lO-in.  HME  was  calculated  to 
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HISS  Phase  2  Recorded  Image 
TLX,  altitude:  30  feet,  Range:  15.7  Nmiles 
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Range  (Nautical  miles) 
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The  following  data  results  are  based  on  an  analysis  of  25  tests  events  where  the  target  was  an  inbound  TLX  towed 
target  and  the  HISS  was  operated  in  a  normal  scan  mode.  Test  events  where  there  was  evidence  that  the  TLX  plumer  did  not 
bum  properly  or  the  HME  did  not  get  up  to  the  required  temperature  were  eliminated  from  this  data  analysis. 
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Using  the  definitions  detailed  previously,  this  figure  compares  the  first  detection,  transient  firm  track,  and  firm  track 
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predicted  results 


The  performance  of  the  system  can  be  predicted  for  other  operating  areas  through  the  use  of  statistical  weather 
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1  Executive  Summary 

The  goal  of  this  study  was  to  investigate  the  feasibility  of  fabricating  a  new  type  of  monolithic 
focal  plane  multipixel  staring  array  for  infrared  (IR)  radiation  detection.  This  sensor  works  by 
the  recently  postulated  process  of  photofluxonic  conversion  in  which  an  infrared  photon  creates  a 
vortex-antivortex  pair  in  a  superconducting  material.  The  novelty  of  the  new  sensor  lies  in  the 
detection  mechanism  itself,  and  in  the  way  the  vortex  (or  antivortex)  is  trapped  in  the  device  al¬ 
lowing  individual  sensing  pixels  to  be  readout  without  direct  electrical  connections.  The  new 
system  offers  tremendous  potential  for  improvement  over  current  infrared  detection  devices  via 
both  the  enhanced  detection  sensitivity  and  the  greatly  reduced  system  complexity. 

Technical  Objectives.  The  overall  objective  of  Phase  I  was  to  determine  the  feasibility  of 
building  an  inductively  coupled,  non-contact  staring  array  based  on  the  new  photofluxonic  de¬ 
tection  elecments.  To  this  end,  we  have  performed  a  series  of  experiments  on  a  single  array  ele¬ 
ment,  and  considered  pixel  requirements  for  constructing  an  array,  as  summarized  below. 

1.  Investigate  the  photofluxonic  effect.  We  have  experimentally  investigated  reproducibility 
of  the  photofluxonic  effect  in  a  2  mm  diameter  loop  of  YBa2Cu307.  We  have  investigated 
the  effect  as  a  function  of  light  intensity,  temperature,  and  frequency  of  chopping  the  IR  light 
source. 

2.  Consider  the  physics  of  the  pixel.  We  consider  the  optimum  coupling  of  the  photodetection 

effect  to  the  flux  state  of  a  pixel  element,  and  consider  the  effect  of  a  vortex  pair  on  the  su¬ 
percurrent  producing  the  Lorentz  force.  We  compare  various  geometries  for  the  supercon¬ 
ducting  array.  ' 

3.  Consider  array  readout  methods.  We  compare  various  methods  for  reading  the  magnetic 
flux  signal  out  using  SQUIDs  or  magneto-optic  effects.  We  consider  cooling  requiremetns, 
and  ways  of  refreshing  the  array. 

4.  Compare  recent  literature.  We  correlate  experimental  results  and  design  considerations 
with  recent  published  scientific  results  for  both  high  and  low  Tc  superconductors,  and  esti¬ 
mate  the  sensitivity  of  the  prototype  to  be  built  in  Phase  II. 

Work  Performed.  We  have  observed  the  conversion  of  quanta  of  light  into  quanta  of  magnetic 
flux  via  vortex-antivortex  pair  creation.  In  addition  we  have  performed  ancillary  experiments 
which  rule  out  thermally  activated  flux  creep  as  a  possible  explanation  of  the  phenomenon  o  :- 
served.  The  principal  technical  successes  resulting  from  this  study  are  as  follows. 

1.  We  have  observed  a  reproducible  change  in  the  flux  trapped  in  a  high  Tc  superconducting 
loop  on  illuminating  a  small  part  of  that  loop  with  light. 

2.  We  have  observed  distinct  steps  in  the  photoresponse  of  high  Tc  thin  film  loops  of 
YBa2Cu307.  The  magnitude  of  these  steps  is  quantized  in  units  which  closely  correspond 
with  the  predicted  value  for  a  single  flux  quantum. 

3.  The  time  dependent  behavior  of  the  response  rules  out  flux  creep  as  a  possible  explanation. 


Potential  Applications.  The  new  sensor  has  applications  in  all  situations  in  which  existing 
cooler  IR  detectors  are  currently  used.  In  addition,  it  has  the  potential  to  greatly  exceed  the  sen¬ 
sitivity,  bandwidth  and  spectral  range  capabilities  of  existing  devices.  Given  the  basic  simplicity 
of  the  individual  pixels  and  the  promise  of  reading  them  out  inductively  (i.e.,  without  making 
electrical  connections)  there  is  tremendous  promise  for  constructing  large  area,  low  cost,  arrays. 
Current  applications  for  infrared  imaging  include:  thermography,  infrared  spectroscopy,  astron¬ 
omy,  and  night  vision  devices  for  the  police  and  military  forces. 

2  Introduction 

Infrared  detectors  have  traditionally  been  semiconductors  devices,  with  both  extrinsic  and  intrin¬ 
sic  detectors  being  used.  In  intrinsic  detectors,  the  incident  photons  cause  transitions  between  the 
valence  and  conduction  bands.  Intrinsic  materials  are  more  suited  to  production  line  processing 
using  standard  photolithographic  techniques.  In.  extrinsic  detectors,  the  incident  photons  cause 
transitions  between  impurity  levels  in  the  bandgap  and  an  adjacent  band.  Extrinsic  device  mate¬ 
rials  have  a  smaller  absorption  coefficient  and  hence  need  detector  thicknesses  of  the  order  of  a 
millimeter  at  meaningful  impurity  concentrations.  In  addition  extrinsic  materials  are  more  noisy 
and  require  more  cooling. 

Superconductors  have  a  lower  energy  gap  compared  to  semiconductors  ( 1  meV  as  compared  to  1 
eV).  This  makes  the  superconductors  more  sensitive  to  radiation,  especially  at  relatively  longer 
wavelengths  (at  which  the  energy  per  photon  is  relatively  small).  Additionally,  superconducting 
devices  have  low  power  dissipation;  when  one  is  dealing  with  large  array  sizes  (which  would 
otherwise  require  a  lot  of  cooling  power)  this  makes  superconducting  arrays  very  attractive.  A 
variety  of  detection  mechanisms  have  been  explored  ranging  from  devices  based  on  the  tempera¬ 
ture  dependance  of  the  superconducting  energy  gap  [1]  to  bolometric  resistive  transition  edge 
thermometers  [2]. 

The  new  sensor  investigated  in  this  program  represents  a  revolution  in  photodetection  technol¬ 
ogy;  both  in  the  physics  of  the  detection  itself  and  in  the  simplicity  of  the  pixel  design,  particu¬ 
larly  for  large  arrays. 

•  The  detection  mechanism  is  b2ised  on  a  process  in  which  the  photon  is  transformed  into  one 
(or  possibly  several)  flux  quanta.  These  quanta  are  stored  and  read  magnetically.  The  de¬ 
tection  speed  is  of  order  100  ps  (10'’°  s).  Once  trapped,  the  quanta  may  be  stored  indefi¬ 
nitely. 

•  The  pixel  design  is  such  that  the  array  may  be  read  out  and  refreshed  inductively  via  mag¬ 
netic  fields;  No  Electrical  Connections  are  required  to  the  array.  In  addition,  the  pixels 
have  No  Power  Dissipation.  This  reduces  the  heat  load  onto  the  array  simply  to  that  of  the 
radiation  field  being  surveyed. 

In  the  recently  postulated  process  of  photofluxonic  conversion,  an  incident  photon  creates  a  vor- 
tex-antivortex  pair  by  local  depression  of  the  superconducting  energy  gap  in  a  current  carrying 
superconductor.  The  supercurrent  is  diverted  around  a  spot  of  scale  of  the  coherence  length, 
which  is  topologically  equivalent  to  a  closely  spaced  vortex-antivortex  pair  superimposed  on  an 
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otherwise  uniform  current  flow  as  shown  in  Figure  1.  Previously,  photofluxonic  processes  have 
only  been  inferred  from  anomalous  resistance  measurements  on  very  long  detector  strips.  In  this 
report  we  present  direct  observation  of  the  creation  of  magnetic  flux  via  the  absorption  of  inci¬ 
dent  optical  photons. 

Each  pixel  of  our  device  will  be  simply  a  loop  of  superconducting  material  of  order  10  -  50  nm  in 
diameter.  A  large  circulating  supercurrent  is  trapped  in  the  superconducting  ring  by  applying  and 
then  removing  a  small  dc  magnetic  field.  The  vortex  pair  is  pulled  apart  by  the  Lorentz  force 
created  by  the  circulating  current  as  shown  in  Figure  2.  One  vortex  (or  antivortex  -  no  absolute 
sign  convention  is  possible)  experiences  force  to  the  right  and  the  other  a  force  to  the  left.  The 
vortices  migrate  under  these  forces  with  the  net  result  being  a  transfer  of  one  flux  quantum  (or 
"anti-quantum")  from  the  left  to  the  right  of  the  current  flow;  i.e.,  into  the  superconducting  ring. 
The  change  in  flux  is  maintained  by  a  small  change  in  the  circulating  supercurrent  in  the  ring. 
The  arrival  of  vortices  (or  antivortices)  at  the  inner  edge  of  the  ring  reduces  the  flux  stored  in  the 
ring  by  an  easily  measurable,  quantized  amount.  In  short,  the  device  is  a  photon  to  fluxon  con¬ 
verter. 

Arrays  based  on  our  device  will  be  much  simpler  to  fabricate  than  any  existing  photodetector. 
The  important  difference  is  that  the  pixel  elements  (the  superconducting  loops)  are  electrically 
isolated  from  each  other,  i.e.,  there  is  no  electrical  interconnects  between  the  elements  or  the  rest 
of  the  system.  This  considerably  simplifies  the  array  fabrication  process.  In  addition,  even 
though  pixel  uniformity  is  desirable,  it  is  not  absolutely  essential  for  the  array  to  function  well. 
As  long  as  the  loops  are  reasonably  similar,  the  critical  currents  will  be  nearly  the  same.  Since 
vortex-antivortex  creation  is  not  a  strong  function  of  the  current  within  the  loop,  there  is  room  for 
some  variation  in  the  loop  parameters  from  pixel  to  pixel.  This  is  a  further  help  to  the  fabrication 
process. 

The  proposed  device  is  actually  the  magnetic  analogue  of  a  photodiode,  and  the  proposed  device 
array  is  the  counterpart  of  a  CCD  used  in  optical  imaging.  The  photofluxonic  model  is  analogous 
to  photoconductivity  through  the  formal  duality  between  electric  charge  and  magnetic  flux  in  two 
dimensions  [3].  A  photon  incident  on  the  junction  region  of  the  sensor  creates  an  excited  vortex- 
antivortex  pair  (c.f.  electron-hole  pair)  which  is  pulled  apart  by  a  current  flowing  through  the 
junction  (c.f.  voltage  for  photoconductive  sensors).  The  quantized  magnetic  flux  created  (c.f. 
electric  charge)  is  stored  in  the  superconducting  ring  (c.f.  the  capacitors  in  a  CCD  array)  before 
being  read.  The  fundamental  advantage  over  CCD  photosensors  is  that  the  greatly  reduced  en¬ 
ergy  gap  for  vortex  pair  creation  in  the  superconductor  based  sensor  opens  up  the  opportunity  for 
new  regimes  of  IR  detection. 
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Figure  1.  Creation  of  a  Vortex-Antivortex  Pair  From  a  Non  Superconducting 
Region. 
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Figure  2.  Separation  and  Trapping  of  Vortex-Antivortex  Pair. 
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3  Task  1:  Determine  Pixel  Requirements. 

The  major  part  of  the  project  was  to  expefimentally  observe  vortex-antivortex  pair  creation  in  a 
superconductor  to  determine  requirements  for  an  array  of  pixels.  By  a  new  method  of  detecting 
the  created  vortex  (or  antivortex)  we  hoped  to  circumvent  many  of  the  problems  previously  en¬ 
countered  by  other  researchers.  The  optical  characteristics  of  potential  HTS  materials  will  dictate 
if  arrays  can  be  fabricated  with  dimensions  and  pixel  sizes  within  usable  limits. 

3.1  Superconducting  System 

Although  the  photofluxonic  model  was  first  proposed  to  explain  the  infrared  response  of  granular 
low  Tc  superconducting  materials  [4,5]  it  has  also  been  applied  to  the  so  called  nonbolometric 
response  of  high  Tc  superconductors.  Since  cooling  to  80  K  is  considerably  more  convenient 
than  cooling  to  5  K,  there  is  a  significant  advantage  in  fabricating  the  new  sensor  from  high  Tc 
materials.  In  order  to  better  characterize  the  interaction  of  the  pixels  with  the  incident  light,  we 
made  measurements  at  different  temperatures. 

Two  high  Tc  films  of  passivated  YBa2Cu307  (YBCO)  were  obtained  from  IBM,  one  consisting 
of  loops  of  3  mm  diameter  and  the  other  2  mm  diameter  loops.  The  HTS  chips  are  composed  of 
a  YSZ  substrate  with  the  HTS  superconductor  thin  film  loops,  40  microns  wide,  and  0.2  microns 
thick.  The  substrate  is  nearly  transparent  to  visible  and  near  IR;  therefore,  most  of  the  light  that 
doesn’t  hit  the  HTS  film  will  pass  through  instead  of  being  thermally  absorbed  by  the  substrate. 

3.2  Experimental  Apparatus 

During  the  course  of  the  experiment,  a  probe  was  built  to  evaluate  experimentally  the  effects  dis¬ 
cussed  in  the  Phase  I  work  plan.  The  instrument  is  illustrated  in  Figure  3.  The  instrument  is  in¬ 
stalled  in  an  aluminum  liquid  helium  dewar.  The  upper  portion  of  the  instrument  is  composed 
mostly  of  “nonmagnetic”  type  304  stainless  steel.  The  lower  assembly.  Figure  3b,  is  the  crux  of 
the  experiment.  The  bottom  of  the  probe  are  free  of  all  magnetic  alloys,  being  composed  of  brass 
and  GIO,  and  a  few  other  materials  (lead,  niobium,  etc.)  named  below. 

The  entire  experiment  is  enclosed  in  a  superconducting  lead  shield.  This  shield  is  completely 
closed  at  the  bottom,  and  mostly  closed  at  the  top;  it  has  some  holes  to  accommodate  screws  and 
structure.  Inside  of  the  shield  is  a  GIO  vacuum  cell.  The  cell  is  mounted  to  a  brass  flange  with 
brass  screws  and  utilizes  an  indium  vacuum  seal.  The  GIO  cell  has  6  mm  thick  side  walls  and  a 
1.5  mm  thick  “window”  at  the  bottom.  The  machined  surfaces  of  the  GIO  are  further  protected 
against  vacuum  leaks  and  permeability  by  a  thin  coat  of  epoxy. 

Wound  around  the  outer  diameter  of  cell  are  100  turns  of  multifiliment  Cu/NbTi  wire.  This  sin¬ 
gle  layer  bias  solenoid  is  centered  with  respect  to  the  gradiometer.  This  solenoid  supplies  about 
45  Oe/amp  at  the  center  and  is  used  to  induce  a  current  in  the  HTS  superconducting  film. 

Inserted  into  bottom  of  the  GIO  cell  on  the  helium  bath  side  is  a  first  order  gradiometer  with  6 
turn,  9  mm  diam  windings  on  a  baseline  of  6  mm.  This  gradiometer  is  made  of  #40  NbTi  wire. 
The  top  of  the  gradiometer  is  butted  up  against  the  GIO  cell’s  “window”.  This  gradiometer  is 
attached  to  a  Quantum  Design  dc  SQUID.  The  SQUID  is.  controlled  and  data  is  acquired  by  a 
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Quantum  Design  dcSOOO  digital  controller  linked  to  a  486  computer  via  GPIB  interface.  An 
analog  output  on  the  dcSOOO  was  utilized  by  a  spectrum  analyzer  and  an  oscilloscope. 

Inside  of  the  cell  (vacuum  space)  is  the  HTS  chip  platform.  The  platform  is  mounted  such  that 
the  HTS  chip  will  be  3  mm  from  the  upper  pickup  loop  of  the  gradiometer,  centered  in  the  lead 
shield.  The  platform  is  thermally  isolated  from  the  rest  of  the  instrument  by  three  small  thin- 
walled  GIO  tubes.  The  platform  is  also  shielded  from  thermal  radiation  by  two  layers  of  superin¬ 
sulation.  With  this  isolation  about  8  mW  of  power  is  required  to  raise  and  sustain  the  platform 
temperature  from  the  4.2  K  ambient  temperature  to  80  K. 

Mounted  on  the  platform  are  a  germanium  resistance  thermometer,  a  heater  winding,  and  the  end 
of  a  fiber  optic.  The  fiber  optic  is  0.4  mm  in  diameter  and  is  held  in  a  Teflon  tube.  The  fiber  is 
locked  into  position  by  a  #0-80  brass  set  screw  that  gently  squeezes  the  Teflon  tube.  The  fiber  is 
positioned  such  that  the  end  is  almost  flush  with  the  bottom  of  the  platform.  The  chip  is 
mounted  centered  on  the  bottom  with  the  fiber  emitting  light  on  a  portion  of  the  superconducting 
loop. 

The  heater  winding  is  composed  of  #40  twisted  pair  phosphor-bronze  wire.  The  twisted  pair 
was  soldered  at  one  end  and  wound  around  (50  turns,  two  layers)  the  platform.  This  corresponds 
to  a  heater  resistance  of  212  Q.  Under  the  heater  wire  is  a  aluminum  split  ring  to  ensure  uniform 
heating  of  the  platform.  This  split  ring  is  heat  sunk  to  the  platform  using  GE703 1  heat  conduct¬ 
ing  varnish,  as  is  the  germanium  thermometer  leads.  The  thermometer  itself  is  inserted  into  a 
slip  fit  hole  with  vacuum  grease  filling  any  gap. 

The  top  of  the  dewar  is  composed  of  two  sets  of  feedthrus:  one  for  the  vacuum  space  and  another 
for  the  helium  vapor  space.  A  NW  vacuum  flange  system  was  employed  for  the  vacuum 
feedthrus,  allowing  versatility  in  changing  the  setup.  The  system  is  composed  simply  of  a  NW 
flange  pipe  cross  (4  ports)  bulkhead  mounted  to  the  top  plate.  One  side  port  consists  of  a  flange 
mounted  circular  multipin  Deteronics  vacuum  feedthru  that  carries  the  4  thermometer  leads  (#36 
phosphor-bronze)  and  the  two  heater  leads  (#36  phosphor-bronze).  The  other  side  port  has  a 
vacuum  valve  to  facilitate  pump  dow'n  and  venting.  The  top  port  has  an  epoxy  vacuum  feedthru 
flange  for  the  fiber  optic.  All  the  vacuum  space  wire  pairs  and  the  fiber  optic  were  run  down  in 
#24  Teflon  tubing.  The  vapor  feedthrus  consist  of  the  following:  the  dc  SQUID  flex  probe  mul¬ 
tipin  connector,  a  4  pin  connector  for  the  helium  level  meter,  and  a  BNC  feedthru  for  the  bias 
solenoid. 

Custom  electronics  were  manufactured  to  ensure  minimal  noise  and  ground  loop  problems.  Each 
circuit  had  a  unique  battery  power  source.  The  following  circuits  were  made:  a  selectable  low 
noise  current  source  (10  pA,  100  pA,  or  1  mA)  for  the  thermometer;  a  low  noise  current  source 
for  the  bias  coil;  and  a  voltage  source  for  the  heater  windings.  The  thermometer  was  read  out 
using  a  battery  powered  digital  multimeter.  All  external  cables  were  coaxial,  with  the  exception 
of  the  He  level  meter  and  the  SQUID  cables. 
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The  pickup  coil  detects  quantized  changes  in  flux  from  photons  striking  the  loop  surface  and 
transmits  this  information  to  the  SQUID  coupled  with  the  SQUID  electronics,  a  Quantum  Design 
Model  5000  dcSQUID  Controller.  The  final  output  was  plotted  on  an  XY  chart  recorder,  and/or 
collected  digitally  by  a  simple  computer  program  written  in  C  which  reads  the  output  of  the 
dc5000  via  a  IEEE  488  card. 

An  optical  isolation  bench  was  configured  with  a  laser  source  directed  through  a  variable  density 
filter  and  a  neutral  density  filter  of  one  percent  transmission  for  intensity  control.  A  fiber  coupler 
with  a  lOX  microscope  objective  was  used  to  align  this  reduced  intensity  beam  with  the  single 
mode  fiber  optic.  The  fiber  feeds  through  a  potted  epoxy  joint  at  the  top  of  the  probe  with  the 
lowest  3-5  cm  terminating  in  the  fixture  at  the  bottom  of  the  probe.  This  end  of  the  fiber  was 
positioned  approximately  0.5  mm  above  the  high  Tc  superconducting  loop  to  be  illuminated. 

3.3  Biasing  the  Superconducting  System 

The  Current  Required  to  Bias  the  System.  A  large  electrical  current  flowing  normal  to  the 
vortex-antivortex  pair  is  required  to  overcome  the  attractive  force  between  the  vortices.  Kadin 
has  postulated  that  the  electrical  current  density,  J,  required  to  pull  apart  the  vortex-antivortex 
pair  should  be  of  order  the  critical  current  density,  Jc  [6].  The  maximum  deviation  below  the 
critical  current  for  which  a  pair  will  be  separated  depends  on  the  frequency /of  the  radiation  in¬ 
volved,  and  the  characteristic  size  of  the  vortices.  Kadin  gives  the  following  expression,  in 
which  Jo  is  the  mean  field  critical  current  density  and  d  the  film  thickness: 


The  characteristic  size  a  depends  strongly  on  the  nature  of  the  superconducting  state;  for  a  ho¬ 
mogeneous  film  it  is  of  order  the  Ginzberg  Landau  coherence  length  for  an  intergranular  array 
of  Josephson  junctions  it  is  of  order  the  grain  size.  Unfortunately,  high  Tc  films  exhibit  varying 
degrees  of  granularity  in  their  conduction  properties  depending  on  their  quality.  For  photons  of 
wavelength  600  nm,  equation  1  gives  an  acceptable  deviation  from  the  critical  current  density  of 
order  lO’*  A/m^  for  a  homogeneous  film  and  10*  A/m^  for  a  granular  film.  Since  typical  values 
of  Jo  are  of  order  lO'*  A/m^  we  see  that,  provided  that  the  bias  current  density  is  close  to  the 
transport  critical  current  density  of  lO"  A/m^  ,  there  should  be  sufficient  force  to  separate  the 
pairs. 

Establishing  the  Current.  The  simplest  way  to  establish  current  densities  of  this  magnitude  in  a 
closed  superconducting  loop  is  to  trap  a  magnetic  field  within  the  loop.  There  are  two  ways  to  do 
this: 

a)  Cool  the  loop  through  its  superconducting  transition  in  the  presence  of  a  magnetic  field. 
Then  remove  the  field.  The  density  of  magnetic  flux,  B,  produced  by  this  procedure  is  shown 
schematically  in  Figure  4a. 

b)  Apply  a  magnetic  field  of  magnitude  twice  that  required  to  just  penetrate  the  film.  Then  re¬ 
move  the  field.  As  shown  in  Figure  4b  this  results  in  the  same  distribution  of  magnetic  flux. 
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In  the  above  analysis  wc  have  used  the  critical  state  model  in  which  pinning  sites  within  the  su¬ 
perconducting  material  create  a  gradient  in  the  magnetic  flux  density  within  the  sample  [7],  This 
gradient  is  formally  equivalent  to  a  constant  current  density  equal  to  Jc  flowing  throughout  the 
sample  cross-section.  The  magnitude  of  the  applied  field  required  to  trap  the  current  is  a  factor 
of  two  greater  when  the  flux  is  pushed  into  the  film  over  when  the  film  is  cooled  in  the  applied 
field. 

There  are  several  potential  complications  possible  with  this  simple  model.  Firstly  it  assumes  a 
homogeneous  material.  Secondly,  it  requires  the  flux  to  uniformly  flow  out  of  the  sample  on  re¬ 
moval  of  the  field  and  not  leave  by  catastrophic  flux  jumps.  In  addition,  for  approach  a),  it  is 
necessary  for  the  sample  to  uniformly  expel  the  flux  within  it  on  cooling.  The  major  experimen¬ 
tal  complication  arises  if  an  insufficient  magnetic  field  is  applied  in  the  flux  trapping  step.  If  the 
applied  field  is  less  than  that  at  which  flux  just  penetrates  across  the  body  of  the  loop  to  reach  the 
center,  then  upon  removal  one  will  trap  a  complex  flux  and  current  distribution.  A  vortex- 
antivortex  pair  created  in  a  region  of  uniform  current  will  be  pulled  apart  as  required.  However, 
on  reaching  the  boundary  between  two  regions  of  oppositely  circulating  current,  a  vortex  will  be¬ 
come  trapped.  The  magnitude  of  the  flux  change  in  the  ring  will  not  be  (po,  but  some  small  frac¬ 
tion  of  it. 

Experimental  Verification.  To  verify  that  the  high  Tc  loops  could  indeed  support  a  bulk  super- 
current,  and  to  investigate  the  reproducibility  of  trapping  such  a  current,  a  series  of  measurements 
was  performed  in  a  commercial  SQUID  magnetometer. 

Initially,  a  simple  sequence  of  measurements  was  performed  at  0  gauss,  with  measurements  taken 
at  5  Kelvin  increments  starting  at  100  K  down  to  50  K.  This  verified  a  transition  temperature  at 
just  above  80  K  for  the  YBCO.  Subsequent  sequences  were  created  to  study  various  effects  in 
the  ring,  including  persistence  and  flux  creep.  The  5  mm  square  wafers  were  mounted  such  that 
the  magnetic  field  was  perpendicular  to  the  plane  of  the  loop. 

We  investigated  both  methods  of  trapping  a  circulating  current.  In  the  first  the  sample  was 
w'armed  to  110  K,  i.e.,  above  the  superconducting  transition  temperature,  the  desired  field  ap¬ 
plied  and  the  sample  cooled  to  5  K.  The  sample's  magnetic  moment  was  measured  and  then  the 
applied  field  reduced  to  zero  while  maintaining  a  temperature  of  5  K  and  a  second  measurement 
made.  In  the  second  test  run  the  sample  was  held  at  a  constant  5  K  and  measurements  were  taken 
as  follows:  apply  a  field,  measure  the  susceptibility,  apply  the  next  higher  field,  etc.,  until  a  field 
of  300  gauss  was  reached.  The  field  is  then  reduced  in  steps  to  0  G,  to  look  for  a  permanent  dia¬ 
magnetic  moment. 

The  results  of  the  second  type  of  measurement  on  the  YBCO  loop  are  shown  in  Figure  5.  We  see 
that,  as  expected,  the  diamagnetic  moment  of  the  loop  increases  as  the  applied  field  is  increased. 
We  see  that  upon  removal  of  the  applied  field  the  sample  retains  a  permanent  magnetic  moment. 
This  is  due  to  a  trapped  supercurrent  flowing  around  the  loop.  As  the  field  was  lowered,  at  an 
applied  field  of  approximately  100  G  the  high  Tc  loop  is  no  longer  stable  in  the  applied  field. 
The  pinning  potential  is  unable  to  support  the  flux  gradient,  flux  "jumps"  into  the  center  of  the 
loop,  and  a  supercurrent  flows  in  the  opposite  direction  around  the  inside  circumference  of  the 
loop  to  maintain  the  flux  within  it.  The  magnetic  moment  is  thereby  reduced.  As  the  external 
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a)  Temperature  Cycling 


Apply  Field 
—Sample  Normal 


Cool  Somple  in  Field 
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Remove  Field 
-Flux  Trapped 


b)  Constant  Sample  Temperature 


Apply  Field,  Sample  Superconducting 
-Internal  Shielding 


Figure  4.  Two  Methods  of  establishing  uniform  circulating  current  densities. 
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field  is  increased  further  the  sample  undergoes  additional  flux  jumps  but  its  maximum  magnetic 
moment  always  remains  of  the  same  order. 

It  is  clear  that  the  preferred  method  of  trapping  flux  is  the  method  of  applying  a  field  to  the  al¬ 
ready  cold  superconductor.  This  is  because  the  method  of  cooling  in  an  applied  field  involves 
waiting  for  thermal  equilibrium  in  the  pixel  array,  which  in  general  would  result  in  a  much  longer 
dead  time  for  the  detector.  One  scenario  under  which  one  could  imagine  a  heating  method  would 
be  feasible  is  in  the  case  of  pixels  of  small  mass  that  are  well-coupled  to  a  relatively  massive 
substrate.  If  the  pixels  were  small  enough,  one  could  heat  them  with  a  short  laser  pulse  to  insti¬ 
gate  the  normal  state  at  the  same  time  the  applied  field  is  raised  to  a  value  where  the  trapped  cur¬ 
rent  is  a  maximum.  This  method  of  “charging”  the  detector  will  be  considered  further  in  Phase  n 
work. 


Calculation  of  the  Circulating  Current  Trapped  in  the  HTS  Loop.  One  may  calculate  the 
induced  circulating  current  from  the  measured  magnetic  moment.  In  Figure  5  the  maximum 
screening  current  in  the  sample  occurs  at  an  applied  field  of  150  G.  The  measured  moment  at 
this  field  is  equal  to  8  x  10'^  emu  =  8  x  10'®  Am^.  The  area  of  the  loop  is  3. 1  x  10’^  m^  implying 
a  circulating  current  I  of  mm/A  =  0.025  A.  By  using  the  relation  =  t  b  Jc,  with  the  thickness 
of  the  film  t  =  0.2  pm  and  the  width  b  =  0.004  cm,  the  critical  current  of  the  film  L  =  3  x  10^ 
A/cm^.  From  Figure  5  we  see  that  the  residual  moment  of  7  x  10'®  Am^  upon  removing  the  ap¬ 
plied  field  implies  a  current  flowing  of  0.038  A,  or  a  current  density  of  4.5  x  10*  A/cm“.  This 
compares  favorably  with  the  vendor-supplied  value. 

It  is  instructive  to  calculate  the  magnetic  field  trapped  in  the  loop  on  removal  of  the  applied  field. 
The  inductance  of  the  loop  is  given  by  [8]  L  =  0.0047ta  [loge  (8a^)  -  1/2].  This  equation  was 
used  to  describe  the  YBCO  loop  where  the  radius  a  =  0.1  cm  and  the  width  b  =  0.004  cm.  L  = 
0.012  pH  was  the  resulting  value  of  the  inductance.  We  can  estimate  the  magnetic  field  created 
by  a  current  flowing  in  the  loop  from  B  =  L  Imax  /  A,  where  B  is  the  magnetic  field  in  gauss,  L 
and  Imax  are  the  inductance  and  maximum  current  of  the  loop  as  described  above,  and  A  is  the 
area  of  the  loop.  Using  these  expressions  we  find  that  a  trapped  current  of  0.02  A  creates  a  mag¬ 
netic  field  of  about  1  G  at  the  center  of  the  loop.  For  reference  it  should  be  noted  that  a  maxi¬ 
mum  of  approximately  2x10*  flux  quanta  may  be  trapped  in  the  high  Tc  loops  used  in  this  study. 


3.4  Calculation  of  the  coupling  factors 


It  is  important  to  calculate  the  flux  coupled  in  to  the  SQUID  for  a  1  tpo  change  in  the  flux  state  of 
the  high  Tc  loop  in  order  to  compare  the  experimental  results  with  theory.  Simple  circuit  analysis 
leads  to  the  following  expression  for  the  flux  coupled  into  the  SQUID: 


(pSQ  - 


Msq  Mp 
L  (  Lsq  +  Lp) 


(po 


The  self-inductance  and  mutual  inductance  of  the  SQUID  were  Lsq  =  1.85  pH  and  Msq  =  0.01 
pH,  respectively.  The  self  inductance  of  the  2  mm  diam  YBCO  loop  is  0.006  pH  and  the  self 
inductance  of  the  NbTi  pickup  coil  was  2  pH.  ( (po=  2x10''*  Wb,  the  numerical  value  of  one 
flux  quanta.) 
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Moment  vs.  Field  for  HTC  Loop 
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Figure  5.  Trapping  a  magnetic  moment  within  the  sample 
by  applying  a  fleld  and  then  removing  it. 
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The  model  used  to  calculate  the  mutual  inductance  Mp  between  the  HTS  loop  and  SQUID  is  that 
of  a  multiturn  loop  and  a  coaxial  circular  filament,  where  the  filament  is  in  a  plane  perpendicular 
to  the  axis  of  the  solenoid  [8]: 

M  =n- f  --jA-a 

where  the  radius  a  of  the  high  Tc  loop  is  0.1  cm,  the  radius  A  of  the  NbTi  coil  is  0.45  cm,  n  is  the 
number  of  turns  on  the  coil,  /is  a  function  of  a,  A,  and  the  separation  between  the  pickup  coil 
and  the  pixel.  The  resulting  mutual  inductance  of  the  NbTi  coil  and  the  high  Tc  loop  was 
0.002  pH. 

Using  a  value  of  0.002  pH  for  M  leads  to  a  coupling  ratio  of  8.9  x  10  '^  <l)o  for  a  one  flux  quantum 
change  in  the  sample  loop.  The  flux  to  voltage  conversion  factor  for  the  SQUID  is  1.26  <])o  =  IV 
giving  a  theoretical  voltage  step  size  of  0.7  mV. 


4  Task  2:  Analyze  Array  Detection  Issues 
4.1  Experimental  Parameters 

To  our  knowledge  there  was  little  experimental  work  performed  by  others  in  the  period  between 
writing  our  Phase  I  proposal  and  beginning  our  experimental  investigation.  Thus,  we  had  little 
information  to  reduce  the  large  range  of  experimental  parameters.  The  issues  of  photon  coupling 
into  an  array  and  the  heat  load  of  the  IR  signal  were  measured  during  this  Phase  I  work  on  a  sin¬ 
gle  pixel  element.  Measurements  on  arrays  of  elements  will  be  performed  in  Phase  II  work. 

Measurement  Bandwidth.  One  of  the  major  concerns  was  how  much  filtering  it  would  be  pru¬ 
dent  to  use.  The  vortex  crossing  occurs  very  rapidly,  on  the  time  scale  of  nanoseconds,  so  a  di¬ 
rect  observation  of  the  event  in  progress  using  a  SQUID  readout  was  impractical.  However,  it  is 
important  to  observe  distinct  individual  crossings  in  order  to  verify  their  amplitude.  To  do  so 
accurately  requires  as  high  a  measurement  bandwidth  as  possible,  which  was  48  kHz  for  our 
SQUID  electronics. 

In  our  probe,  the  flux  coupled  into  the  SQUID  for  one  90  change  in  the  high  Tc  loop  was  about 
10'^  9o-  Looking  for  a  10'^  tpo  change  in  a  48  kHz  measurement  bandwidth  requires  a  probe  noise 
of  about  5x10'^  <Po/  VHz.  This  is  about  the  white  noise  levels  of  the  dc  SQUID  used  when  cou¬ 
pled  to  a  coil  of  similar  inductance.  However,  this  is  the  white  noise  level  for  a  fully  shielded 
SQUID,  for  which  there  is  a  negligible  contribution  from  the  environment.  When  connected  to  a 
coil  there  is  always  a  contribution  from  ambient  noise  and  motion  in  the  earth's  magnetic  field. 
For  the  maximum  flexibility  it  was  decided  not  to  roll  off  the  pickup  coil  circuit  at  all  but  attempt 
to  do  all  filtering  in  the  SQUID  electronics. 

4.2  Experimental  Results 

4.2.1  Observation  of  the  Photofluxonic  Effect 

As  soon  as  the  SQUID  was  operational  we  obtained  a  clear  change  in  the  flux  state  of  the  YBCO 
loops  on  applying  laser  light.  Figure  6  shows  the  data  from  two  runs  taken  with  the  YBCO  film 
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at  5  K.  The  field  used  to  trap  the  circulating  current  was  15  mT  corresponding  to  a  circulating 
current  of  order  5  x  10*  A/cm^.  The  light  from  the  1.3  |i,m  semiconductor  diode  laser  was  turned 
on  for  30  seconds,  then  off  for  30  seconds.  The  data  was  collected  in  the  computer  through  an 
A/D  board  at  a  rate  of  6  kHz.  The  data  was  reduced  by  averaging  in  50-point  blocks.  Figure  6 
clearly  indicates  the  presence  of  steps  in  the  SQUID  response.  However,  in  addition  to  quantized 
steps  there  is  a  clear  drift  which  occurs  only  while  the  light  is  on  and  for  a  short  period  (1-5  sec¬ 
onds  depending  on  the  circulating  current)  immediately  after  it  is  blocked.  This  same  behavior 
was  seen  in  many  other  runs  at  5  K  similar  to  those  shown  here.  Several  runs  with  light  levels 
about  1%  of  those  shown  in  Figure  6  showed  similar  background  noise  levels,  but  no  drift  was 
observable  within  the  resolution  of  our  measurements.  Test  runs  without  the  HTS  film  present 
produced  no  increase  in  the  drift  signal  when  the  light  wa.s  present.  Furthermore,  tests  with  the 
film  present  but  with  no  circulating  current  within  it  also  gave  zero  response.  Given  its  clear 
magnetic  nature,  the  simplest  interpretation  of  the  drift  is  that  it  is  due  to  flux  vortices  which  be¬ 
come  trapped  in  the  film  and  move  in  small  steps  rather  than  pass  across  the  body  of  the  loop  is  a 
single  sweep.  These  steps  are  too  small  to  be  resolved  by  our  apparatus. 

The  main  question  is  whether  the  distinct  steps  correspond  to  one  flux  quantum  in  the  regions  in 
which  the  flux  state  is  changing.  We  expect  the  trapping  of  a  single  vortex-antivortex  pair  to 
produce  a  dc  voltage  step  of  order  0.7  mV.  Figure  6a  shows  a  flux  jump  of  about  0.4  mV,  and 
Figure  6b  a  jump  of  about  0.7  mV.  These  numbers  are  consistent  with  an  interpretation  in  which 
the  former  is  a  single  vortex-antivortex  pair  being  created,  and  the  latter  a  pair  of  photofluxonic 
transitions  occurring  too  close  together  in  time  for  us  to  resolve  with  our  data  collection  limita¬ 
tions.  We  see  a  remarkable  quantized  nature  in  the  amplitude  of  the  steps  demonstrating  that  we 
are  indeed  measuring  the  arrival  of  quantized  units  of  magnetic  flux  into  the  center  of  our  high  Tc 
loop. 

Another  question  which  arises  in  interpreting  the  data  is  why  the  incidence  of  apparent  photo¬ 
fluxonic  vortices  is  so  low.  The  drift  rate  of  the  signal  in  the  presence  of  the  laser  light  indicates 
that  a  large  number  of  photons  are  being  absorbed.  Studying  the  efficiency  of  different  materials 
for  producing  the  photofluxonic  effect  will  be  a  major  Task  in  Phase  II  work. 

In  another  set  of  measurements,  we  used  an  electronic  circuit  to  chop  the  laser  light  at  various 
frequencies.  Typical  results  are  shown  in  Figure  7.  For  this  data,  the  light  was  chopped  at  a  fre¬ 
quency  of  10  Hz.  In  this  Figure,  we  see  that  the  general  drift  that  was  present  in  the  unchopped 
light  data  is  still  present,  which  indicates  that  the  film  is  responding  to  an  average  power  input. 
The  large  jump  occurred  when  the  light  was  on,  and  has  the  same  magnitude  as  the  other  jumps 
shown  in  Figure  6;  thus  we  interpret  this  as  also  occurring  through  the  photofluxonic  effect. 

In  yet  another  set  of  measurements,  we  increased  the  temperature  of  the  pixel  to  100  K,  applied  a 
field  of  12  mT,  and  cooled  the  HTS  loop  in  the  field  to  60  K.  The  field  was  then  removed  and 
measurements  taken  with  the  laser  light  on  for  60  seconds,  and  off  for  60  seconds.  The  results 
are  shown  in  Figure  8.  During  these  measurements,  the  temperature  was  allowed  to  drift.  We 
see  that  there  is  a  large  drift  in  the  signal  corresponding  to  the  changing  temperature.  Any  jumps 
due  to  the  photofluxonic  effect  we  unfortunately  not  distinguishable  due  to  this  large  drift.  In 
Phase  n  we  will  improve  the  temperature  control  of  our  test  apparatus  to  allow  for  measurements 
at  constant  elevated  temperature. 
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Figure  6.  Photoresponse  Data  for  YBCO  Loop  at  5  K,  trapping  Field  at  12  mT.  The  laser 
was  on  for  the  first  30  seconds  of  data  collection,  then  turned  off  for  30  seconds. 
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Figure  7.  Photoresponse  of  YBCO  Loop  at  5  K,  with  the  laser 
light  chopped  at  a  frequency  of  10  Hz. 


Figure  8.  Photoresponse  of  YBCO  Loop  at  60  K,  with  laser  light  on  for  60  seconds,  then  off 
for  60  seconds.  The  large  drift  is  correlated  with  the  changing  temperature  of  the  pixel. 
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4.2.2  Flux  Creep 

It  is  important  to  consider  alternative  explanations  for  the  behavior  obserx'ed.  The  most  likely 
alternative  phenomenon,  and  perhaps  the  only  other  possible  explanation  for  the  observed  pho¬ 
toresponse,  is  thermally  activated  flux  creep.  Essentially,  bolometric  heating  of  the  film  raises  its 
temperature  thereby  activating  fluxons,  trapped  when  the  circulating  current  is  set,  to  jump  past 
their  pinning  sites.  In  the  flux  creep  model  there  is  no  creation  of  vortex-antivortex  pairs;  rather 
the  photons  free  flux  quanta  already  trapped  in  the  film. 

With  the  laser  blocked  a  large  level  of  flux  creep  clearly  took  place  each  time  a  new  current  was 
set  in  the  HTS  loop.  Indeed  immediately  after  trapping  a  higher  circulating  current  the  sample 
response  is  so  dominated  by  flux  creep  that  the  photoresponse  is  almost  unmeasurable;  depend¬ 
ing  on  the  current  level  it  was  necessary  to  wait  between  approximately  one  hour  for  the  non- 
illuminated  jump  and  drift  signal  to  fall  to  acceptable  levels.  The  flux  creep  rate  has  been  quan¬ 
tified  in  earlier  measurements  by  us  by  placing  a  similar  sample  in  a  field  of  100  G  at  5  K  in  a 
SQUID  magnetometer  and  measuring  its  moment  every  minute  for  60  minutes.  The  field  was 
then  dropped  to  0  G  and  again,  60  measurements  were  made  a  minute  apart.  After  an  initial 
rapid  decrease  the  creep  rate  was  approximately  lO  ’  emu/s  =  I0‘‘^  Am^/s  after  an  hour.  This 
corresponds  to  approximately  5000  tpo  which  is  at  least  a  factor  of  500  too  large  to  allow  accurate 
measurements  of  photofluxonic  processes.  At  70  K  the  creep  rate  was  a  fact  of  3  greater. 

The  crucial  question  as  far  as  photofluxonics  is  concerned  is  the  following:  is  the  flux  change 
signal  when  the  laser  is  on  due  to  the  creation  of  vortex-antivortex  pairs  or  simply  the  thermally 
activated  escape  of  trapped  flux?  The  following  general  observation  can  be  made:  The  general 
drift  rate  builds  up  in  the  first  minute  after  beginning.  The  photoresponse  then  proceeds  by  peri¬ 
ods  of  decreasing  drift  followed  by  "flux  avalanches"  in  which  groups  of  trapped  vortices  appear 
to  move  at  random  intervals  in  bundles  over  a  period  of  1  to  10  seconds. 

In  discussing  flux  creep  it  is  important  to  note  that  we  really  do  not  know  the  thermal  time  con¬ 
stant  Td  for  this  system.  There  are,  however,  two  possibilities:  that  the  sample  reaches  thermal 
equilibrium  on  a  time  scale  much  shorter  than  our  experiments  (i.e.,  Td  is  of  order  10'^  to  lO' s), 
or  that  the  sample  temperature  generally  changes  during  the  period  of  illunaination.  Even  given 
the  uncertainty  between  these  two  cases  the  following  arguments  may  be  given  against  a  flux 
creep  interpretation. 

Al:  if  Xd  is  short  then  one  expects  the  film  to  be  raised  to  a  higher  temperature  with  higher  light 
intensity  and  thus  take  longer  to  return  to  normal. 

A2:  if  Td  is  long  then  this  clearly  contradicts  a  thermal  model  in  which  the  film  temperature  in¬ 
creases  over  the  time  scale  of  our  experiments  (approximately  60  s). 

A3:  The  observation  that  the  drift  rate  builds  up  after  the  light  is  turned  on  is  more  consistent 
with  a  phenomenon  in  which  the  number  of  flux  quanta  increases  with  time  than  with  one  based 
on  the  freeing  of  trapped  quanta  since  in  the  latter,  the  number  of  trapped  vortices  is  actually  a 
maximum  at  the  beginning  of  the  measurement  period  and  should  not  build  up. 
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In  summary,  we  believe  that  the  jumps  in  our  data  are  inconsistent  with  thermally  activated  mo¬ 
tion  of  flux  already  trapped  in  the  sample.  However  this  is  a  question  of  some  subtlety  and 
should  be  studied  further  in  Phase  n  (note  that  an  IR  sensor  based  on  thermally  activated  flux 
creep  may  be  of  great  practical  interest  anyway,  though  not  as  exciting  scientifically). 

Effect  of  Varying  the  Light  Intensity 

For  scientific  reasons  it  is  obvious  that  the  sample  photoresponse  should  in  some  way  depend  on 
the  intensity  of  the  illumination.  Frome  Figure  6,  it  is  clear  we  see  a  clear  increase  in  the  pho¬ 
toresponse  with  increasing  illumination.  Further  measurements  demonstrated  that  the  increase  is 
monotonic.  The  data  at  low  intensities  is  unreliable  owing  to  variations  in  the  background  drift, 
and  the  fact  that  the  total  flux  change  is  only  of  order  1-5  flux  quanta  which  means  that  multi¬ 
quantum  steps  have  a  relatively  large  effect. 

Effect  of  Varying  the  Circulating  Current 

A  second  major  question  is  how  the  photoresponse  is  affected  by  the  circulating  current  in  the 
film.  The  variation  of  the  photoresponse  with  intensity  was  measured  for  several  values  of  the 
trapping  field.  We  see  a  clear  increase  in  the  response  with  increasing  circulating  current.  There 
appears  to  be  some  saturation  of  the  photoresponse  with  trapping  field  at  the  higher  light  intensi¬ 
ties. 

The  value  of  the  circulating  current  can  only  be  inferred  from  the  trapping  field  and  is  not  pre¬ 
cisely  known.  We  see  from  Figure  5  that  the  flux  trapped  in  the  sample,  and  hence  the  circulat¬ 
ing  current,  is  not  a  linear  function  of  the  trapping  field  but  tends  to  a  maximum.  Thus  the  lack 
of  increase  in  the  photoresponse  with  trapping  field  could  simply  be  due  to  the  circulating  current 
maximizing  at  the  sample  critical  current. 

However  the  circulating  current  does  not  begin  to  maximize  for  applied  fields  below  about  10 
mT,  whereas  at  the  higher  intensities,  apparent  saturation  in  the  photoresponse  occurs  for  fields 
above  3  mT.  It  is  tempting  to  suggest  that  the  flattening  of  the  photoresponse  is  due  to  an  in¬ 
crease  in  the  efficiency  of  the  process.  However,  at  the  present  time  such  a  claim  would  be  pre¬ 
mature.  This  question  will  addressed  in  much  greater  detail  during  Phase  11.  Similarly  the  lack 
of  any  photoresponse  with  a  trapping  field  of  about  0.3  mT  may  indicate  that  for  low  circulating 
currents  there  was  insufficient  Lorentz  force  to  pull  apart  a  vortex-antivortex  pair.  This  must 
also  be  addressed  in  Phase  U. 

4.2.3  Comparison  with  Recent  Research  by  Others 

Limited  efforts  toward  demonstrating  the  photofluxonic  effect  via  resistance  measurements  have 
been  made  by  groups  at  Rochester  [10]  and  TRW.  In  a  resistance  measurement,  the  movement 
of  the  vortices  under  strong  current  bias  results  in  an  induced  voltage  across  the  film  which  can 
be  measured.  However  it  must  be  borne  that  all  voltage  based  measurements  suffer  from  the 
problem  that  it  is  very  difficult  to  separate  the  voltage  due  to  movement  of  the  vortex  pair  from 
the  myriad  other  voltages  of  similar  magnitude  which  are  almost  certain  to  be  present. 

Spurious  voltage  signals  could  arise  from  changes  in  film  resistance  due  to  bolometric  heating, 
the  motion  of  partially  pinned  fluxons  in  nearby  regions  of  the  superconductor,  and  insufficient 


Quantum  Magnetics,  Inc. 


18 


Final  Report 
F49620-94-C-0058 


screening  of  environmental  magnetic  noise  which  changes  the  local  magnetic  field  at  the  junc¬ 
tion.  In  fact,  the  difficulty  of  an  unambiguous  subtraction  of  all  possible  spurious  junction  volt¬ 
ages  probably  accounts  for  lack  of  direct  observations  of  the  photofluxonic  effect  in  Josephson 
junctions  at  the  present  time. 

In  the  time  since  our  original  Phase  I  proposal,  the  only  notable  recent  publication  regarding 
photofluxonics  has  come  from  Walko  and  Van  Harlingen  [11]  who  proposed  an  alternative  de¬ 
tection  scheme.  While  this  scheme  is  similar  to  ours  in  that  a  SQUID  is  used  to  detect  the  vor¬ 
tex-antivortex  pair  magnetically,  the  vortices  are  measured  while  moving  and  no  attempt  is  made 
to  trap  the  flux  created. 

There  have  been  several  basic  experimental  investigations  into  the  nonbolometric  IR  response  of 
superconducting  films.  Taken  as  a  whole,  or  in  some  cases  individually,  these  results  are  incon¬ 
clusive  regarding  the  question  of  vortex-antivortex  pair  creation  by  photons.  Ghis  et  al.,  [12] 
have  observed  the  now  familiar  "ultra  fast  nonbolometric  response"  followed  by  a  slower  bolo- 
metric  decay  in  the  voltage  IR  response  of  YBCO  thin  films.  Chera  et  al.,  [13]  have  found  that 
the  microwave  response  of  T1  and  Bi  based  high  Tc  thin  films  displays  both  bolometric  and  non¬ 
bolometric  components  but  is  predominantly  the  latter,  while  the  effect  of  a  cw  HeNe  laser  is 
mainly  bolometric. 

Matsuda,  et  al.  have  shown  that  the  motion  of  vortex-antivortex  pairs  accounts  well  for  the  tem¬ 
perature  dependence  of  the  Hall  number,  the  diagonal  pxx  resistivity,  and  Hall  pxy  resistivities  in 
single  unit  cell  thick  films  of  YBCO.  It  is  postulated  that  a  Kosterlitz-Thouless  (KT)  transition 
[14]  occurs.  Below  the  KT  transition  temperature  vortex-anti  vortex  pairs  are  bounded  by  a 
quasi-long-range  order.  Above  the  KT  transition  the  pairs  dissociate  in  the  presence  of  a  finite 
transport  current.  It  is  central  to  the  argument  that  both  vortices  and  antivortices  are  present.  In 
the  model  the  vortex-antivortex  pairs  are  excited  by  thermal  fluctuations  and  are  extremely  nu¬ 
merous;  there  are  about  ten  times  more  pairs  than  vortices  induced  by  application  of  a  10  T  mag¬ 
netic  field.  These  most  recent  results  will  be  utilized  when  planning  the  research  and  develop¬ 
ment  to  be  performed  in  Phase  II. 

4.3  Pixel  Configuration 

The  heart  of  an  infrared  imaging  system  is  the  detector  focal  plane.  The  focal  plane  detector 
consists  of  individual  detector  elements  or  pixels.  The  detected  signal  from  each  pixel  is  com¬ 
bined  and  analyzed  to  obtain  an  image  of  the  source  using  standard  mathematical  techniques. 
For  actual  detectors,  it  is  important  that  the  detector  noise  dominates  the  noise  of  subsequent 
electronics  such  as  preamplifiers.  Also,  the  detector  should  be  sufficiently  fast  to  accommodate 
the  necessary  system  bandwidth. 

At  this  very  early  stage  it  is  difficult  to  be  very  specific  about  the  design  of  the  pixels.  Our 
overall  concept  for  the  array  is  shown  schematically  in  Figure  9.  The  photosensitive  detection 
regions  (pixels)  are  electrically  isolated  superconducting  rings.  Given  their  simplicity  there  are 
three  basic  questions; 

a)  What  is  the  optimum  high  temperature  superconductor  for  this  application. 
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b)  What  is  the  optimum  pixel  size:  what  diameter  and  track  width  should  the  loop  have  and 
what  thickness  of  film  to  use. 

c)  How  can  the  coupling  of  light  into  the  detection  region  be  optimized. 

Selection  of  the  High  Tc  Material.  Given  the  quantum  nature  of  the  detection  process,  selection 
of  the  optimum  high  Tc  film  really  requires  a  detailed  assessment  of  the  intrinsic  noise  mecha¬ 
nisms  within  the  detector  material.  For  example,  semiconductor  detectors  are  limited  by  the 
presence  of  a  background  current  due  to  minority  carriers.  For  photofluxonic  detection  the 
dominant  noise  source  is  most  likely  flux  creep. 

The  conventional  way  to  reduce  the  level  of  flux  creep  is  to  increase  the  number  of  flux  pinning 
centers  in  the  material  or  by  lowering  its  operating  temperature.  However,  making  the  trapped 
flux  less  mobile  also  increases  the  probability  that  a  vortex-antivortex  pair  will  become  trapped 
in  the  film  instead  of  being  pulled  apart.  Thus  decreasing  flux  creep  may  actually  prevent  one  of 
the  pair  being  swept  into  the  center  and  recorded  as  a  detection  event. 

One  strategy  would  be  to  use  a  high  Tc  film  with  very  little  flux  pinning.  In  previous  studies  we 
observed  a  greater  ratio  of  signal  to  background  drift  in  TBCCO  films  than  in  YBCO,  although 
this  cannot  be  seen  as  a  definitive  study.  However,  it  is  well  known  that  TBCCO  films  are  less 
effective  at  pinning  flux  than  YBCO  which  suggests  that  reducing  pinning  may  help.  Depending 
on  the  field  trapped,  the  level  of  flux  creep  could  initially  be  high  but  one  could  wait  or  devise 
flushing  procedures  to  reduce  its  effect.  However,  a  material  with  a  very  low  level  of  pinning 
would  also  have  a  rather  low  critical  current  density.  Thus  the  question  is  then:  what  is  the 
minimum  current  necessary  to  separate  a  vortex-antivortex  pair?  At  the  present  time  the  answer 
to  this  is  unknown.  One  of  first  tasks  of  Phase  n  would  be  to  make  measurements  with  several 
different  superconducting  materials  to  determine  in  a  systematic  way  the  current  necessary  to 
separate  vortex-antivortex  pairs.  One  clear  test  of  the  effect  of  trapped  vortices  would  be  to  use  a 
type  I  superconductor,  which  cannot  trap  vortices  in  the  same  manner  as  a  type  II  superconductor. 
Different  procedures  for  depositing  the  films  and  annealing  them  could  also  be  tried,  to  see  the 
effect  on  vortex  trapping. 
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Figure  9.  Schematic  of  Array  Concept. 
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The  Optimum  Pixel  Size.  One  way  to  reduce  the  chance  of  a  vortex  pair  becoming  trapped  in 
the  film  is  to  make  the  width  of  the  loop  and/or  the  thickness  of  the  film  as  small  as  possible. 
The  optimum  values  for  the  pixel  dimensions  will  depend  on  a  number  of  factors,  including  the 
required  resolution,  sensitivity  and  fabrication  constraints.  In  general  one  would  like  to  make  the 
pixel  dimensions  as  small  as  possible.  The  following  tradeoffs  are  apparent. 

Pixel  Diameter  Decreasing  individual  pixel  size  will  increase  the  pixel  density  allowing  for  more 
resolution  of  the  image  and  reducing  the  cooling  power  requirements.  The  trade-offs  are 
that  small  loops  are  more  difficult  to  fabricate.  In  addition,  a  smaller  loop  will  have  a 
smaller  inductance  reducing  the  number  of  flux  quanta  that  may  trapped  in  the  loop  before 
the  loop  current  goes  supercritical.  This  can,  in  effect,  reduce  the  available  imaging  time, 
before  the  array  has  to  be  recharged  by  applying  a  magnetic  field  to  reestablish  the  circulat¬ 
ing  current.  The  inductance  of  the  pixel  loop  is  approximately  proportional  to  its  diameter 
and  so  10  pm  loops  will  have  an  inductance  300  times  less  than  our  samples.  The  circulat¬ 
ing  current  depends  on  the  linewidth  and  thickness  of  the  loop  (see  below)  so  a  3  pm  x  1 
pm  cross  section  would  be  approximately  30  times  less  than  used  for  our  experiments. 
Thus  we  might  expect  to  trap  of  order  100  fluxons  in  such  a  pixel.  Whether  such  a  stored 
flux  is  adequate  depends  on  the  system  application.  However  it  seems  unlikely  that  pixels 
smaller  than  this  will  be  feasible. 

Pixel  Linewidth.  As  mentioned  above,  reducing  the  pixel  width  reduces  the  chance  that  the  vor¬ 
tex  pair  will  become  trapped  within  the  body  of  the  loop  and  not  reach  its  center.  Of 
course  the  loop  width  cannot  be  larger  than  its  radius  and  so  for  small  diameter  pixels  this 
is  a  further  reason  to  make  the  linewidth  small.  Advantages  of  increasing  the  width  include 
fabrication  issues,  getting  a  higher  circulating  current,  and  the  need  to  maximize  the  photon 
absorption  by  using  as  large  a  film  area  as  possible. 

Film  Thickness  In  this  case  one  must  trade  off  the  desire  for  a  large  circulating  current  with  the 
issue  of  unwanted  trapping  of  the  vortex  pair.  One  must  also  consider  the  absorption  depth 
of  the  photon  within  the  high  Tc  material.  This  is  discussed  further  below. 

Coupling  Radiation  into  the  Detector.  As  for  many  of  the  above  pixel  design  issues,  maximi¬ 
zation  of  the  coupling  of  radiation  into  the  sensor  requires  a  detailed  understanding  of  the  physics 
of  the  process.  At  this  very  early  stage  we  have  little  experimental  data  on  which  to  base  our  dis¬ 
cussion.  Since  we  believe  the  photofluxonic  process  to  be  basically  caused  by  local  heating  it 
should  be  possible  to  increase  the  efficiency  of  photon  absorption  by  depositing  absorption  pads 
of  antireflection  coatings  on  the  superconducting  loop.  Alternatively,  one  can  arrange  for  the 
photon  to  be  absorbed  in  the  substrate.  For  example,  detectors  have  been  constructed  from  su¬ 
perconducting  films  deposited  on  GaAs  in  which  the  photons  are  incident  on  the  GaAs  substrate. 
Noise  equivalent  powers  of  3  X  lO  '^  W/VHz  have  been  achieved  on  a  200  element  series  array 
10  X  10  )im“  junctions  [15].  The  choice  of  absorption  material  depend  strongly  on  the  wave¬ 
length  of  the  radiation  of  interest. 

Although  we  are  not  able  to  proceed  in  detail  with  the  pixel  design,  from  the  above  discussion 
and  our  experimental  results  we  are  able  to  deduce  the  following  overall  guidelines; 
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•  The  basic  pixel  configuration  will  be  a  simple  single  turn  loop. 

•  The  minimum  pixel  size  will  be  of  order  10  |im. 

•  Thallium  or  bismuth  based  high  Tc  materials  may  be  better  photofluxonic  detectors. 

•  These  general  observations  will  be  investigated  in  detail  during  the  research  in  Phase  II. 

5  Task  3:  Analyze  Pixel  Bias  Issues 
5.1  Refreshing  the  Pixels 

A  description  of  the  biasing  technique  used  in  Phase  I  work  was  presented  in  Section  3.1.3 
above.  In  biasing  the  pixels  one  can  choose  between  either  having  discrete  magnetic  field  coils 
behind  each  pixel  in  order  to  set  them  up  individually  or  whether  to  use  a  single  large  coil  around 
the  entire  array  to  perform  a  global  refresh.  It  is  not  necessary  for  the  coils  to  be  cooled  and  so 
there  would  be  no  additional  heat  load  on  the  refrigeration  system.  However,  biasing  the  pixels 
individually  requires  that  the  magnetic  field  coils  be  sufficiently  close  to  the  high  Tc  pixels  that 
one  coil  affects  only  the  desired  array  element.  Thus  using  discrete  coils  will  increase  the  system 
complexity  and  may  interfere  with  the  readout  scheme.  Accordingly  we  will  use  a  global  refresh 
method  in  Phase  11. 

6  Task  4:  Analyze  Array  Readout  Issues 
6.1  The  Array  Configuration 

Present  IR  focal  plane  array  technology  involves  the  assembly  of  elements  which  consist  of  a 
high  impedance  semiconducting  detector  followed  by  a  semiconducting  amplifier  [15].  The  out¬ 
puts  are  mixed  together  to  produce  an  analog  signal.  The  signal  corresponds  to  the  sequential 
interrogation  of  the  infrared  energy  falling  on  each  of  the  individual  detector  elements.  The  de¬ 
tectors  are  usually  fabricated  in  a  two-dimensional  array  on  one  substrate  and  the  amplifier  and 
mixer  circuitry  in  a  corresponding  array  on  a  second  substrate.  The  elements  are  connected  to¬ 
gether  by  means  such  as  indium  bump  bonds  sandwiched  between  the  two  substrates. 

It  is  a  nontrivial  task  to  fabricate  large  arrays  of  semiconducting  elements  on  a  single  substrate 
having  uniform  electrical  characteristics.  It  is  difficult  to  achieve  high  yields  because  of  the 
highly  precise  doping,  etc.  that  is  required.  This  limits  single  array  size;  larger  arrays  consist  of  a 
group  of  single  arrays  connected  together.  A  problem  with  using  large  semiconducting  arrays  is 
the  amount  of  heat  that  is  generated.  Noise  considerations  require  that  thermal  imaging  detectors 
be  operated  in  the  40  -  100  K  range  [15].  In  addition  to  this,  the  cooling  system  has  to  absorb  the 
heat  generated  by  the  operating  device. 

Superconducting  arrays  generate  little  or  no  heat;  current  low  Tc  systems  have  to  be  cooled  to 
lower  temperatures  than  semiconducting  arrays,  but  require  less  cooling  power  to  be  applied  be¬ 
cause  of  the  reduced  heat  load.  Superconducting  arrays  can  use  different  types  of  detection 
mechanisms.  Present  superconductor-based  device  technologies  include  devices  based  on  bolo- 
metric  detection  and  the  Josephson  effect. 
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Bolometric  detection  can  be  difficult  to  achieve  in  an  array,  because  of  the  need  for  temperature 
uniformity  (to  a  millidegree  or  so)  over  the  array.  The  sharp  resistance  transitions  of  conven¬ 
tional  superconductors  will  make  a  detector  difficult  to  achieve.  High  temperature  superconduc¬ 
tors  have  a  broader  transition,  and  have  been  used  to  detect  IR  radiation  [2],  but  conventional  su¬ 
perconducting  electronics  cannot  be  used  (because  of  temperature  considerations)  [15]. 

The  absence  of  interconnects  between  the  pixel  elements  of  our  detector  allows  us  great  flexibil¬ 
ity  in  the  design  of  the  array.  There  is  no  need  to  attempt  to  fabricate  very  large  arrays  on  a  sin¬ 
gle  substrate  which  can  lead  to  difficulties  in  fabrication.  Instead,  modest  size  arrays  (10  X  10  or 
even  5X5)  can  be  made  and  arranged  together  to  make  a  larger  array.  Large  arrays  are  needed 
to  maximize  the  detector  area  for  detecting  weak  signals.  Also,  in  order  to  properly  image  the 
source,  it  is  advantageous  to  collect  information  over  as  wide  an  area  as  possible. 

Construction  of  arrays,  per  se,  will  thus  be  very  straightforward.  At  the  present  time,  almost  all 
IR  detector  arrays  are  cooled  to  cryogenic  temperatures  and  so  the  need  to  cool  the  proposed  HTS 
materials  will  not  increase  the  system  complexity.  For  the  Phase  n  prototype,  we  will  mount  the 
array  directly  on  the  cold  finger  of  a  cryocooler  as  is  customary  for  IR  detectors.  We  will  use 
some  form  of  integrated  dewar  assembly  and  so  the  new  system  will  be  functionally  equivalent  to 
existing  IR  arrays.  The  other  issues  to  consider  are  how  to  bias  the  pixels  (see  Sections  3.1.3  and 

5. 1  above),  and  how  to  read  their  state.  This  issue  is  discussed  below. 

6.2  The  Pixel  Readout  Scheme 

One  of  the  innovations  of  the  proposed  new  photodetector  is  the  possibility  of  using  a  non- 
contact  inductively  coupled  readout  scheme  by  detecting  the  magnetic  field  present  in  the  high  Tc 
loops.  The  signal  is  simply  tpo  divided  by  the  loop  area,  which  is  approximately  2  x  10'®  T  and  2 
X  10'^  T  for  1  mm  and  10  pm  loops  respectively.  Hall  effect  sensors  and  magneto-optic  devices 
are  two  possible  options  to  use  as  readout  devices.  Each  of  these  are  magnetic  field  detectors  and 
should  be  able  to  measure  the  field  due  to  the  flux  trapped  in  the  superconducting  loops  at  least 
for  smaller  loop  diameters.  They  each  have  their  own  advantages  and  disadvantages  which  will 
be  discussed  below. 

Hall  effect  sensors  [16]  can  be  thin  film  and  can  be  integrated  into  the  IR  array.  A  convenient 
way  to  fabricate  them  is  to  start  with  a  silicon  wafer  that  is  suitably  doped  for  low  temperature 
Hall  operation.  If  the  doping  is  on  the  surface,  individual  Hall  sensors  can  be  made  by  etching 
patterns  in  the  silicon  wafer.  The  current  and  sensing  leads  can  be  patterned  on  the  wafer  using, 
say,  aluminum  thin  film  deposition  and  subsequent  photolithography.  The  superconduc;;ng  focal 
plane  array  can  then  be  deposited  on  the  other  side  of  the  wafer.  If  properly  aligned,  c.:  jh  indi¬ 
vidual  Hall  sensor  will  be  positioned  on  the  back  of  one  pixel  element. 

This  technique  of  using  Hall  sensors  can  be  inconvenient  because  of  the  large  number  of  electri¬ 
cal  leads  which  may  be  necessary  to  drive  and  sense  the  various  Hall  elements.  In  addition  to  the 
heat  leak  through  these  electrical  lines,  the  operation  of  the  Hall  sensors  will  generate  additional 
heat  to  be  overcome  by  the  refrigerator/cryogen.  The  main  problem  is  sensitivity.  Typical  de¬ 
vices  can  only  detect  fields  down  to  10'^  T  and  so  will  only  work  for  pixels  smaller  than  50  [im  in 
diameter. 
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An  alternate  approach  is  to  use  magneto-optic  sensors.  A  plane-polarized  beam  of  light  incident 
on  a  magneto-optic  material  has  its  plane  of  polarization  rotated  upon  transmission  through  it; 
the  amount  of  rotation  is  proportional  to  the  thickness  of  the  material  and  the  ambient  magnetic 
field.  The  constant  of  proportionality  is  material  dependent  and  is  known  as  the  Verdet  constant. 
This  is  the  Faraday  effect  [17].  There  is  a  related  effect  influencing  the  reflected  light  called  the 
Kerr  magneto-optic  effect.  For  example,  a  normally  incident  plane  polarized  light  becomes  el- 
liptically  polarized  on  reflection,  the  degree  of  the  shift  in  the  plane  of  polarization  being  de¬ 
pendent  on  the  ambient  magnetic  field.  In  one  possible  detection  scheme,  the  Kerr  effect  can  be 
exploited  to  detect  the  field  at  the  superconducting  loops  by  shining  light  on  a  magneto-optic 
material  placed  close  to  the  superconducting  loops.  Field  sensors  based  on  the  Faraday  and  Kerr 
effects  have  been  made  and  used  [18,19]. 

The  advantage  of  using  magneto-optics  is  that  the  heat  load  into  the  array  sensor  can  be  reduced 
significantly,  since  fiberoptics  can  be  substituted  for  the  electrical  leads  that  will  have  to  connect 
up  the  Hall  sensors.  There  is  a  strong  possibility  that  a  single  fiberoptic  line  can  scan  the  whole 
array,  further  reducing  the  heat  load.  The  pixel  readout  sensor  will  likely  be  simpler,  the  possi¬ 
bility  of  using  a  continuous  magnetooptic  film  on  the  backside  of  the  pixel  array  exists.  If  real¬ 
ized,  this  will  further  reduce  the  fabrication  complexity  of  the  array.  Lasers  to  scan  the  pixel  ar¬ 
ray  should  be  inexpensively  available;  technology  for  this  exists,  for  example,  in  magnetic  stor¬ 
age  devices  and  compact  disc  players. 

To  estimate  the  sensitivity  we  note  that  the  magnetic  field  of  the  loop  will  extend  into  the  mag¬ 
neto-optic  to  a  distance  approximately  equal  to  its  diameter.  For  a  material  with  a  very  high 
Faraday  rotation  such  as  YIG,  we  expect  rotations  of  5  x  10"*  and  0.05  degrees  for  1  mm  and  10 
pm  loops  respectively.  These  should  both  be  detectable  with  existing  apparatus. 
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INTERNET  LOCATIONS 

Note:  The  following  URLs  are  current  as  of  the  date  of publication. 

Optical  Detection  Laboratory  -  http://www.opt-scLarizona.edu/summaries/Eustace_Dereniak/homepage.html 
The  primary  focus  of  the  laboratory's  research  is  the  detection  of  light,  not  necessarily  in  the  visible  spectrum,  and 
conversion  of  detected  light  to  electrical  signals.  Under  the  umbrella  of  the  Optical  Detection  Laboratory  are  five 
individual  laboratories  working  on  different  detection  technologies. 

Infrared  Information  Analysis  Center-  IRIA  -  http://www.erim.org/IRIA/iria.html 

The  Infrared  Information  Analysis  Center  (IRIA)  was  established  in  1954  at  the  Willow  Run  Laboratories  to 
facilitate  the  exchange  of  information  within  the  Department  of  Defense  (DoD)  infrared  community.  The  center  is 
now  a  DoD  Information  Analysis  Center  (lAC)  sponsored  by  the  Defense  Technical  Information  Center  and 
monitored  by  the  Associate  Director  for  Science  and  Technology,  Army  Night  Vision  and  Electronic  Sensors 
Directorate. 

IRIA’s  mission  is  to  collect,  analyze,  and  disseminate  information  on  infrared  and  electro-optical  (IR/EO) 
technology  with  an  emphasis  on  military  applications. 

International  Society  for  Optical  Engineering  -  http://www.spie.org/ 

A  non-profit  professional  society  dedicated  to  advancing  research,  engineering,  and  applications  in  optics, 
photonics,  imaging,  and  electronics. 

National  Research  Council  of  Canada-Photonic  Systems  Group  -  http://alpha.ps.iitnrc.ca/ 

Photonics,  the  synthesis  of  optics  and  electronics  for  data  transmission  and  manipulation,  is  an  area  of  pivotal 
importance  in  information  processing  and  it  is  a  key  enabling  technology  for  the  information  highway. 

University  of  Maryland, 

Photonics  Switching  and  Integrated  Optoelectronics  Lab  -  http://www.ee.umd.edu/photonics 

Research  in  the  Photonics  Lab  spans  a  wide  range  of  topics  in  the  field  of  Integrated  Optoelectronics.  These  include 
semiconductor  lasers  and  amplifiers,  integrated  receivers,  optoelectronic  interconnects  and  packaging,  wavelength 
conversion  devices,  optical  bistability,  semiconductor  laser  reliability  studies  and  laser  theory  and  modeling. 

Select  “Related  Links”  to  get  all  the  following  Laboratories  and  Research  Centers: 

NIST 

Jet  Propulsion  Laboratory 

David  Samoff  Center 

Lawrence  Livermore  Laboratories 

AT&T  Bell  Laboratories 

Sandia  National  Laboratories 

Army  Research  Laboratories 

Rome  Labs  Surveillance  and  Photonics  Home  Page 

AT&T  Bell  Laboratories 

Sandia  National  Laboratories 

Optical  Sciences  Center,  University  of  Arizona 

Laboratory  for  Optics  and  Quantum  Electronics,  University  of  Virginia 

Center  for  Photonics  Research,  Boston  University 

Institute  of  Optics,  University  of  Rochester 

Center  for  Quantized  Electronic  Structures,  University  of  California,  Santa  Barbara 
Photonic  Systems  Group,  University  of  Illinois  at  Urbana  Champaign 
CREOL,  University  of  Florida 
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Hyperspectrum  News  Letter  -  http://www.techexpo.com/WWW/opto-knowledge/hyperspectrum/hypersp3.html 
Technology  Insertion  Engineering: 

the  Advanced  Technology  Dilemma  -  http://www.mcclellan.af.mil/TI/TIE/homepage.html 


Links  to  research  in  specialized  areas  including: 


Advanced  Electronics  Technology  Center  AETC 
Microelectronics  Technology  TIEFB 
Systems  and  Technology  Integration  TIEBD 
**Electro-Optics  Technology  TIEFC  -  http://www.mcclellan.af.mil/TI/TIEAaser.html 


Guide  to  NASA  Online  Resources: 

Alphabetical  List  of  Resources  -  http://naic.nasa.gov/nic/guide/alphabeticaLhtml 

Guide  to  NASA  Online  Resources  -  Alphabetical  List  of  Resources.  This  page  is  a  simple  listing  of  the  resources  in 
this  Guide  by  name.  If  you  are  not  sure  of  the  name,  but  know  you’d  like  to  find  a  resource  that  can  be  generally 
classified  as  scientific  try  this  URL  site. 


Sample  Papers  and  Documents  Available: 

Optical  and  Infrared  Interferometric  Imaging  -  http://dopey.cattech.edu/compastro/bispectrum_reporthtmi 
Scientific  Progress  in  Optical  and  Infrared  Interferometric  Imaging  Report  prepared  July  1990:  1.  Overview  of 
Techniques  2.  Scientific  Objectives  3.  Computational  Aspects  of  Bispectral  Imaging  4.  Results  5.  References 
6.  Publications. 


Near-Infrared  Spectroscopy  -  http://www.scimedia.com/chem-ed/spec/vib/nir.htm 

Near-Infrared  Absorption  Spectroscopy  (NIR)  Introduction  NIR  spectroscopy  is  the  measurement  of  the  wavelength 
and  intensity  of  the  absorption  of  near-infrared  light  by  a  sample.  Near-infrared  light  spans  the  800 
nm  -  2.5  um  (12,500  -  4000  cm-1)  range  and  is  energetic  enough  to  excite  overtones  and  combinations  of  molecular 
vibrations  to  higher  energy  levels. 


Infrared  Technology  -  http:/Aancelotgcs.redstone.army.milAthtml 

Infrared  Technology  Functions  Perform  the  common  functions  with  respect  to  techniques,  components,  devices, 
subsystems,  and  systems  which  function  in  the  infrared  (IR)  region  of  the  spectrum  such  as  focal  plane  arrays, 
detectors,  scanners,  radomes,  optics. 
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■^AD-A311078 

ARIZONA  STATE  UNIV  TEMPE 

(U)  Scale-Based  Techniques  in  Secure  and  Jam-Resistant 
Communication  Systems. 

DESCRIPTIVE  NOTE:  Final  rept.  Nov  92-Jan  96 
JUL  96  29P 

PERSONAL  AUTHORS:  Cochran,  Douglas 
CONTRACT  NO:  F49620-93- 1-0051 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  This  three-year  project  has  focused  on 
theoretical  development  and  preliminary  performance 
analysis  for  wavelet-based  multiple  access  communication 
systems.  Sets  of  orthonormal  wavelet  symbols  suitable  for 
asynchronous  multiple  access  spread  spectrum 
communications  were  developed.  Investigation  of  possible 
countermeasures  to  wavelet-based  communication  systems 
led  to  development  of  a  cyclostationary  signal  detector  for 
which  detection  thresholds  corresponding  to  desired 
false-alarm  probabilities  can  be  determined  analytically.  This 
detector  was  applied  against  both  traditionally  modulated  and 
low  probability  of  detection  signals  including  wavelet 
modulated  signals.  Some  basic  research  on  wavelet  analysis 
performed  in  connection  with  this  project  led  to  a  class  of 
multidimensional  generalizations  of  the  wavelet  transform. 

DESCRIPTORS:  (U)  *ANTIJA]V[MING,  *SECURE 
COMMUNICATIONS,  *SPREAD  SPECTRUM, 
♦MULTIPLE  ACCESS,  DETECTION,  DETECTORS, 
PERFORMANCE  TESTS,  THRESHOLD  EFFECTS, 
SIGNALS,  COMMUNICATION  AND  RADIO  SYSTEMS, 
COUNTERMEASURES. 

IDENTIFIERS:  (U)  *MULTIPLE  ACCESS 
COMMUNICATIONS,  *ASYNCHRONOUS 
COMMUNICATIONS. 


AD-A310  906 
ROME  LAB  ROME  NY 

(U)  Anti- Jamming  Optical  Beam  Nuller. 

DESCRIPTIVE  NOTE:  Rept.  for  Oct  92-Sep  95 
MAY  96  56P 

PERSONAL  AUTHORS:  Turbyfill,  Michael  E.;  Lutsko, 
Jeffrey  M. 

REPORT  NO:  RL-TR-96-65 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  The  development  of  an  Acousto-optic  (AO) 
based  anti-jamming  optical  beam  nuller  (A JOB),  including 
laboratory  and  live  radar  tests,  is  presented.  The  purpose  of 
the  AJOB  system  is  the  cancellation  of  multipath  jamming 
interference  in  advanced  surveillance  radars.  AJOB  is  a 
multichannel  adaptive  optical  system  which  performs 
cancellation  of  multiple  wideband  (10  MBz)  interference 
sources  in  the  presence  of  multipath.  The  live  radar  test 
consisted  of  using  a  downconverted  80  MHz  received  signal 
from  the  main  and  subarrays  of  a  C-band  radar  to  correlate 
jamming  signals  produced  by  stationary  jammers.  The 
correlation  parameters  fed  a  tapped  delay  line  filter  to  form  an 
estimate  of  the  noise,  which  was  subtracted  from  the  main 
antenna  signal.  For  the  scenarios  tested,  the  long  integration 
time  for  the  correlation  data  provided  accurate  estimates  of 
the  jammer  delays,  and  therefore  single-step  convergence  was 
achieved. 

DESCRIPTORS:  (U)  *RADAR  ANTIJAMMING, 
♦ACOUSTOOPTICS,  ALGORITHMS,  OPTICAL 
EQUIPMENT,  ACCURACY,  TRANSMISSION  LINES, 
OPTICAL  CORRELATORS,  INTEGRATION, 
INTERFERENCE,  ADAPTIVE  SYSTEMS, 
CONVERGENCE,  NULLS(AMPLITUDE),  RADAR 
ANTENNAS,  OPTICAL  PROCESSING,  BEAM 
STEERING,  MULTIPATH  TRANSMISSION,  PHASED 
ARRAYS,  C  BAND,  MULTICHANNEL,  BRAGG  ANGLE, 
ELECTRIC  FILTERS,  ACOUSTIC  DELAY  LINES. 

IDENTIFIERS:  (U)  PHASED  ARRAY  RADARS,  VCSEL 
(VERTICAL  CAVITY  SIDE  EMITTING  LASERS). 


^  Included  in  The  DTIC  Review,  October  1996 
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AD-A310  484 

NATIONAL  AIR  INTELLIGENCE  CENTER 
WRIGHT-PATTERSON  AFB  OH 

(U)  Developments  in  Military  Infrared  and  Laser  Technology. 
JUL  96  12P 

PERSONAL  AUTHORS:  Huifeng,  Xiong 
REPORT  NO:  NAIC-ID(RS)T-03 13-96 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  CAMA  China 
Astronautics  and  Missilery  Abstracts,  v3  nl  pl73-176  1996. 

ABSTRACT:  (U)  Infrared  and  laser  technology  is  widely 
used  in  aviation  and  space,  in  reconnaissance,  remote  sensing 
and  remote  metering,  missile  warning,  aircraft  night  vision 
and  night  navigation,  and  plays  an  important  role  in  target 
detection,  identification  and  tracking,  aiming  and  precision 
guidance  as  well  as  communications.  Their  application  can 
greatly  increase  the  precision  of  weapons  systems,  increase 
ECM  capabilities,  counter  stealth  capabilities,  counter  surface 
object  interference  capabilities,  greatly  reduce  costs,  and 
increase  the  effectiveness  of  weapons  systems.  Therefore, 
they  are  recognized  as  power  multipliers  for  the  weapons 
with  which  the  military  is  equipped.  The  military  of  all 
countries  are  playing  serious  attention  to  developments  in 
infrared  and  laser  technology,  especially  the  development  of 
elements  components  and  their  basic  technology.  For 
example,  of  the  1 1  different  key  technologies  of  the  United 
States  Defense  Science  and  Technology,  1.162  billion  dollars 
was  invested  in  optics  and  electronics  (primarily  composed  of 
infrared  and  laser  technologies)  during  the  years  1992,  1993 
and  1994  (not  including  SDI).  This  constituted  13  percent  of 
all  investments  in  key  technologies  during  these  three  years. 

DESCRIPTORS:  (U)  *INFRARED  DETECTION, 

*LASERS,  GUIDED  MISSILES,  OPTICS,  UNITED 
STATES,  AIRCRAFT,  DEFENSE  SYSTEMS, 
INVESTMENTS,  WEAPON  SYSTEMS,  TRACKING, 
COSTS,  NAVIGATION,  REMOTE  DETECTORS, 
PRECISION,  TRANSLATIONS  GUIDANCE,  NIGHT, 
NIGHT  VISION,  TARGET  DETECTION,  INFRARED 
RADIATION,  ELECTRONIC  COUNTERMEASURES, 
RECONNAISSANCE,  CHINA,  WARNING  SYSTEMS, 
CHINESE  LANGUAGE. 


AD-A310  003 

NATIONAL  AIR  INTELLIGENCE  CENTER 
WRIGHT-PATTERSON  AFB  OH 

(U)  Applications  of  GPS  in  Airborne  Electronic 
Countermeasure  Reconnaissance. 

APR  96  lOP 

PERSONAL  AUTHORS:  Zhigang,  Zhang 
REPORT  NO:  NAIC-ID(RS)T-0082-96 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  CAMA,  China 
Astronautics  and  Missilery  Abstracts  (China)  v2  n3  p39-43 
1995. 

ABSTRACT:  (U)  When  implementing  electronic  counter 
reconnaissance  or  other  electronic  countermeasure  missions 
on  moving  platforms,  operating  personnel  working  on  the 
platforms  must  grasp  in  real  time  the  exact  position  of  the 
platform  itself.  In  command  posts  or  control  centers,  there  is  a 
need  to  understand,  in  real  time,  the  direction  of  platform 
movements.  When  implementing  the  positioning  of  emitting 
sources,  precise  platform  locations  are  even  more 
indispensable.  In  the  past,  on  aircraft,  reliance  was  put  on 
inertial  navigation  systems  and  aviation  instruments  to 
provide  data  and,  after  processing,  precise  positions.  The 
limitations  associated  with  making  use  of  this  type  of  method 
are  relatively  large;  precisions  are  not  high.  Real  time 
characteristics  are  relatively  bad.  Opting  for  the  use  of  digital 
transmission  navigation  display  systems  based  on  global 
satellite  navigation  systems  avoids  the  shortcomings 
discussed  above.  Moreover,  it  is  possible  to  conveniently 
generalize  application  to  various  types  of  mobile  platforms. 

DESCRIPTORS:  (U)  *AERIAL  RECONNAISSANCE, 
♦GLOBAL  POSITIONING  SYSTEM,  *ELECTRONIC 
COUNTERMEASURES,  POSITION(LOCATION),  REAL 
TIME,  MOTION,  CONTROL  CENTERS,  AERONAUTICS, 
PLATFORMS,  MOBILE,  MISSIONS,  PRECISION, 
NAVIGATION  SATELLITES,  ARTIFICIAL  SATELLITES, 
INSTRUMENTATION,  TRANSLATIONS,  CHINA, 
CHINESE  LANGUAGE. 
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AD-A309  111 

NAVAL  AIR  WARFARE  CENTER  AIRCRAFT  DIV 
PATUXENT  RIVER  MD 

(U)  F-14  LANTIRN  Flight  Testing.  The  Cat  is  Back!! 

DESCRIPTIVE  NOTE:  Professional  paper 
APR  96  30P 

PERSONAL  AUTHORS:  Mnich,  William 
UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  The  system  is  currently  operational  on  both 
the  F-15  and  F-16.  Shown  here  on  the  Strike  Eagle  are  the 
key  elements  of  LANTIRN,  which  features  two  independent 
pod-mounted  sensors.  On  the  left  side  is  the  navigation  pod 
containing  both  the  Terrain-Following  Radar  and  Navigation 
FLIR  images  from  which  can  be  presented  to  the  pilot  on  the 
Wide  FOV  raster  HUD  for  lowlevel  navigation.  The  targeting 
pod  is  carried  on  the  other  side,  and  incorporates  the  targeting 
FLIR  and  laser  rangefinder  systems.  With  the  full  two-pod 
system,  the  F15E  can  perform  the  low-altitude  strike  mission 
day  or  night  in  almost  any  weather  conditions.  It’s  a 
magnificent  weapons  system  with  an  outstanding  combat 
record,  and  was  one  of  the  stars  of  Desert  Storm. 
Unfortunately,  the  F-15  can’t  land  on  an  aircraft  carrier.  So, 
why  did  the  Navy  decide  to  incorporate  the  LANTIRN 
targeting  pod  on  the  F-14,  and  how  did  we  do  it? 

DESCRIPTORS:  (U)  *FLIGHT  TESTING,  ^FORWARD 
LOOKING  INFRARED  SYSTEMS,  *NIGHT  FLIGHT, 
*LOW  ALTITUDE,  *TERRAIN  ANALYSIS  RADAR, 
♦NIGHT  RECONNAISSANCE,  DETECTORS,  WEAPON 
SYSTEMS,  LASERS,  NAVIGATION,  EXTERNAL 
STORES,  LOW  LEVEL,  DAY,  TARGETING,  MOUNTS, 
RANGE  FINDING,  STRIKE  WARFARE,  RASTERS, 
TERRAIN  FOLLOWING. 

IDENTIFIERS:  (U)  F-14  AIRCRAFT,  LANTIRN(LOW 
ALTITUDE  NAVIGATION  TARGETING  INFRARED  FOR 
NIGHT). 


AD-A309  765 

NAVAL  AIR  WARFARE  CENTER  AIRCRAFT  DIV 
PATUXENT  RIVER  MD 

(U)  Air  Combat  Environment  Test  and  Evaluation  Facility 
(ACETEF)  Support  for  Multi-Sensor,  Multi-Spectral  Sensor 
Fusion  Testing. 

MAY  96  22P 

PERSONAL  AUTHORS:  Macone,  Dan 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  This  paper  presents  an  outline  of  the 
following:  Description  of  ACETEF  Multi-Sensor/ 
MultiSpectral  Sensor  Fusion  Testing  Issues;  Example 
Applications;  Distributed  Training  and  Test  Applications; 
Summary. 

DESCRIPTORS:  (U)  *DATA  FUSION,  ♦MULTISENSORS, 
DETECTORS,  TEACHING  METHODS,  AERIAL 
WARFARE,  MULTISPECTRAL. 

IDENTIFIERS:  (U)  ACETEF(AIR  COMBAT 
ENVIRONMENT  TEST  AND  EVALUATION  FACILITY), 
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^AD-A309  670 

ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
ARNOLD  AFBTN 

(U)  Direct  Write  Scene  Generation  Development  Efforts  for 
Closed-Loop  Evaluation  of  Focal  Plane  Arrays. 

DESCRIPTIVE  NOTE:  Final  rept.  1  Oct  94-30  Sep  95 
MAY  96  48P 

PERSONAL  AUTHORS:  Steely,  S.  L.;  Fugerer,  R.  H.; 
Lowry,  H.  S.,  Ill;  Holt,  L.  L. 

REPORT  NO:  AEDC-TR-95-34 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  Laser  based  closed  loop  Direct  Write  Scene 
Generation  (DWSG)  methods  are  being  developed  at  AEDC 
to  simulate  dynamic  sensor  operations  and  complex  infrared 
scenes.  New  photonic  image  synthesis  methods  are  being 
developed  to  employ  image  to  object  Whittaker  Shannon 
sampling,  anisoplanatic  optical  convolution  by 
quasiisoplanatic  spatial  decomposition,  and  high  speed  digital 
electronics  for  acoustooptic  modulation.  Because  of  the  large 
terabyte  volume  of  data  to  be  processed,  the  increased 
bandwidth  requirements,  and  the  increased  simulation  fidelity 
required  for  DT&E  and  OT&E  of  focal  plane  array  sensors, 
the  laser  based  DWSG  methodology  is  being  extended  to 
accommodate  optical  and  computational  decomposition 
methods  to  better  exploit  highly  and  massively  parallel  real 
time  image  processing  schemes.  Optical  and  computational 
decomposition  will  not  only  provide  high  fidelity  optical 
simulation  for  anisoplanatic  optical  sensors  and  complex 
infrared  scenes,  but  will  also  facilitate  high  speed  parallel 
processing  schemes  for  real  time  closed  loop  DWSG  and 
sensor  operations.  A  computer  simulation  has  been 
constructed  and  is  being  used  to  develop  and  evaluate  the 
high  speed  computational  algorithms  required  for  scene 
extraction  and  convolutions.  A  Proof  of  Principle 
demonstration  using  existing  DWSG  hardware  also  illustrated 
the  feasibility  of  the  concept. 

DESCRIPTORS:  (U)  *FOCAL  PLANES,  *INFRARED 
HOMING  ALGORITHMS,  COMPUTERIZED 
SIMULATION,  REQUIREMENTS,  FLIGHT  TESTING, 
DIGITAL  SYSTEMS,  COMPUTATIONS,  PARAMETRIC 
ANALYSIS,  COST  EFFECTIVENESS,  REALTIME, 
ARRAYS,  LASER  APPLICATIONS,  INFRARED  IMAGES, 
STATISTICAL  PROCESSES,  BANDWIDTH, 
DECOMPOSITION,  CLOSED  LOOP  SYSTEMS, 
ACOUSTOOPTICS. 

IDENTIFIERS:  (U)  *FOCAL  PLANE  ARRAYS, 
DWSG(DIRECT  WRITE  SCENE  GENERATION), 
WHITTAKER  SHANNON  SAMPLING. 


AD-A309  229 

NORTHWESTERN  UNIV  EVANSTON  IL 
(U)  New  Infrared  Materials  and  Detectors. 

1995  6P 

PERSONAL  AUTHORS:  Razeghi,  M.;  Kim,  U.  D.;  Park,  S. 
U.;Choi,Y.H.;  Wu,D. 

CONTRACT  NO:  N00014-92-U-1951,  $N000 14-94- 1-0902 
UNCLASSIFIED  REPORT 

AVAILABILITY:  Pub.  in  Inst.  Physics  Conference  Sen,  nl45 
c8  pl085-1090  1995.  Available  only  to  DTIC  users.  No 
copies  furnished  by  NTIS. 

ABSTRACT:  (U)  In  this  paper,  we  present  an  overview  of  the 
growth  of  InAs(lx)Sbx,  InAs(l  -x)SB(x),  In(l-x)Tl(x)  Sb, 
and  (InP)lx(TlP3)x  as  potential  materials  for  long 
wavelength  (8-12  microns)  infrared  photodetector 
applications.  Incorporation  of  thallium  into  InSb  and  InP 
resulted  in  a  bandgap  decrease  of  the  alloys.  Thallium 
incorporation  has  been  evidenced  by  various  techniques  such 
as  X-ray  diffraction.  Auger  electron  spectrometry  and 
photoconductivity  measurements.  Preliminary  photodetectors 
fabricated  from  the  grown  materials  are  also  reported.  InAsSb 
photodetectors  showed  room  temperature  photoresponse  up  to 
13  microns. 

DESCRIPTORS:  (U)  *INFRARED  DETECTORS, 
*PHOTODETECTORS,  *INFRARED  OPTICAL 
MATERIALS,  REPRINTS,  HALL  EFFECT,  X  RAY 
DIFFRACTION,  INDIUM  ALLOYS,  FOCAL  PLANES, 
LONG  WAVELENGTHS,  AUGER  ELECTRON 
SPECTROSCOPY,  PHOTOVOLTAIC  EFFECT,  CRYSTAL 
GROWTH,  INFRARED  SPECTRA 
PHOTOCONDUCTIVITY,  THALLIUM,  INDIUM 
COMPOUNDS,  PHOTOSENSITIVITY,  INFRARED 
PHOTOCONDUCTORS. 


^  Included  in  The  DTIC  Review^  October  1996 
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AD-A308  808 

WAYNE  STATE  UNIV  DETROIT  MI  DEPT  OF 
ELECTRICAL  AND  COMPUTER  ENGINEERING 

(U)  Photoconductive  Gain  and  Generation-Recombination 
Noise  in  Quantum  Well  Infrared  Photodetectors. 

FEB  95  6P 

PERSONAL  AUTHORS:  Shadrin,  V,  D.;  Mitin,  V  V; 
Kochelap,  V.  A.;  Choi,  K.  K. 

CONTRACT  NO:  DAAH04-94-G-0097 

UNCLASSIFIED  REPORT 

AVAILABILITY:  Pub.  in  Jnl.  of  Applied  Physics,  v77  n4 
pl771-1776,  15  Feb  95.  Available  only  to  DTIC  users.  No 
copies  furnished  by  NTIS. 

ABSTRACT:  (U)  Photocurrent  and  excess  current  noise  in  a 
quantum  well  infrared  photodetector  are  considered  using  a 
drift-diffusion  model  of  charge  carrier  transport.  The  effect  of 
quantum  well  recharge  under  the  influence  of  the  nonuniform 
generated  charge  carriers  is  addressed.  The  recharging  effect 
drastically  changes  the  dependency  of  both  photoconductive 
gain  and  excess  current  noise  gain  upon  detector  parameters. 
We  have  found  that  for  uniform  generation,  both  gains 
coincide.  For  nonuniform  generation,  noise  gain  is  essentially 
different  from  photoconductive  gain.  This  distinction  is  of  the 
order  of  100%  for  the  real  device  parameters.  The  existing 
discrepancy  in  formulae  for  photoconductive  gain  and  excess 
current  noise  derived  in  different  models,  which  implicitly 
assumed  drift  transport  of  electrons,  is  cleared  up. 

DESCRIPTORS:  (U)  ^INFRARED  DETECTORS, 
♦RECOMBINATION  REACTIONS, 

♦PHOTODETECTORS,  REPRINTS,  PARAMETERS, 
QUANTUM  WELLS,  CHARGE  CARRIERS,  TRANSPORT 
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(U)  Experimental  Study  of  Unusual  IR  and  Far  IR 
Vibrational  Energy  Transfer  Processes  in  Solids. 

DESCRIPTIVE  NOTE:  Final  rept.  1  Aug  92-31  Jan  96 
JAN  96  60P 

PERSONAL  AUTHORS:  Sievers,  A.  J. 

CONTRACT  NO:  DAAL03-92-G-0369 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  Seven  different  experimental  studies  are 
described,  (a)  Hole  burning  in  chalcogenide  glasses:  The 
dephasing  time  is  a  function  of  the  average  coordination 
number  of  the  glass  and  is  independent  of  the  chemical 
composition,  (b)  Vibrational  relaxation  of  diatonic  molecules 
in  alkali  halide  crystals:  The  matrix  isolated  CN-  molecule 
obeys  an  energy  gap  law  while  the  hydride  molecule  does  not. 
(c)  Vibrational  relaxation  of  diatonic  molecules  in  glass:  Our 
IR  pump  -  probe  examination  of  SH  in  As2S3  glass  has 
shown  that  this  impurity  molecule  is  actually  weakly 
hydrogen  bonded  in  the  glassy  structure,  (d)  Pocket 
vibrational  modes  in  crystals:  A  systematic  study  of  a  point 
defect  in  a  crystal  lattice  has  demonstrated  that  there  are 
localized  modes  with  the  maximum  mode  amplitude  at  sites 
far  removed  from  the  defect  itself,  (e)  Intrinsic  localized 
vibrational  modes  in  perfect  crystals:  The  inclusion  of  both 
cubic  and  quartic  terms  in  the  potential  has  produced  not  only 
localized  vibrational  modes  but  also  a  localized  dc  expansion 
at  the  mode  site,  (f)  Extinction  sum  rules  for  particles  of 
arbitrary  size:  We  have  uncovered  a  general  extinction  sum 
rule  which  is  independent  of  the  size  of  the  particle  in  relation 
to  the  wavelength  and  independent  of  the  composite  nature  of 
the  particle,  (g)  Measurement  of  a  charged  particle  bunch 
from  the  coherent  transition  radiation  mm-wave  spectrum: 

The  coherent  far  infrared  radiation  induced  from  relativistic 
electron  bunches  submillimeter  length  provides  a  new  way  to 
characterize  the  electron  bunch  shape. 
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(U)  Spatial-Temporal  Detection  of  Electro-Optic  Moving 
Targets. 

DESCRIPTIVE  NOTE:  Journal  article 
APR  95  lOP 

PERSONAL  AUTHORS:  Pohlig,  Stephen  C. 

REPORT  NO:  JA-6920 
CONTRACT  NO:  F19628-90-C-0002 

UNCLASSIFIED  REPORT 

AVAILABILITY:  Pub.  in  IEEE  Transactions  on  Aerospace 
and  Electronic  Systems,  v3 1  n2  p608-616,  Apr  95.  Available 
only  to  DTIC  users.  No  copies  furnished  by  NTIS. 

ABSTRACT:  (U)  A  maximum  likelihood  algorithm  to  detect 
moving  targets  in  space  time  electrooptic  data  is  derived 
using  a  model  of  temporally  stationary  and  spatially 
nonstationary  clutter  statistics.  Performance  is  evaluated  in 
terms  of  the  probabilities  of  false  alarm  and  detection.  This 
algorithm  is  applied  to  a  variety  of  image  sequences:  visible 
band  and  infrared  (IR)  sensors,  with  terrestrial  and  celestial 
clutter  backgrounds.  Comparison  of  theoretically  predicted 
and  experimentally  derived  statistics  shows  excellent 
agreement,  validating  the  model  and  theoretical  predictions. 
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Infrared  Upgrades. 
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PERSONAL  AUTHORS:  Chase,  R.;  Singkofer,  K. 
CONTRACT  NO:  F30602-91-D-0042 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  This  report  covers  upgrades  made  to  the 
General  Research  Corporation  Advanced  Surveillance 
Testbed  (AST).  Under  this  task,  major  upgrades  were  made  to 
the  AST  phenomenological  models.  These  included  upgrades 
to  the  synthetic  background  models,  development  of  a 
statistical  persistence  model,  and  development  of  a  target 
intensity  model  based  on  CHARM  and  SIRRM.  In  addition, 
selected  AST  signal  processing  and  data  processing  modules 
were  upgraded  to  accommodate  the  newly  defined 
phenomenological  inputs. 

DESCRIPTORS:  (U)  *RECONNAISSANCE  SATELLITES, 
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REPORT  NO:  MS-10596A 
CONTRACT  NO:  F19628-95-C-0002 

UNCLASSIFIED  REPORT 

AVAILABILITY:  Pub.  in  Spaceflight  Mechanics  v89  pl681- 
1690  1995.  Available  only  to  DTIC  users.  No  copies 
furnished  by  NTIS. 

ABSTRACT:  (U)  Space  Surveillance  is  the  task  of  keeping  a 
current  catalogue  of  information  on  man  made  earth  bound 
resident  space  objects  (RSOs)  to  support  military  and  civilian 
needs.  The  mission  areas  include  new  launch  processing, 

RSO  catalogue  maintenance  (including  debris)  and  catalogue 
augmentation,  ASAT  support,  satellite  attack  warning,  RSO 
mission  and  payload  assessment,  and  treaty  monitoring.  Some 
of  the  necessary  functions  to  perform  these  missions  are 
search  and  detection  acquisition,  tracking,  data  processing, 
tasking  and  scheduling,  and  network  integration.  At  present,  a 
mature  network  of  sensors,  radar  and  optical,  provide  the  data 
to  accomplish  these  missions.  Evolving  requirements  suggest 
a  smart,  space  based  sensor  could  enhance  this  system.  The 
Midcourse  Space  Experiment  (MSX)  has  three  sensor 
systems  covering  four  Infrared  bands  (Spirit  3),  the  visible 
band  (SBV),  and  five  UV  and  Visible  bands  (UVISI)  that  will 
provide  space  surveillance  data  to  demonstrate  concepts  for 
space  based  space  surveillance.  For  MSX  there  are  18  space 
surveillance  experiments  in  four  categories.  They  are 
incremental.  First,  there  are  four  sensor  technology  tests  - 
sensor  characterization  and  calibration  -  seeking  to  learn 
about  the  performance  of  this  suite  of  space  based  sensors  in 
stressing  backgrounds  for  search,  target  detection,  clutter 
suppression,  tracking  and  metric,  photometric,  and 
radiometric  accuracy. 
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(U)  Detection  Range  Performance-Horizon  Infrared  Surveil¬ 
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PERSONAL  AUTHORS:  Dezeeuw,  Patrick  A. 

REPORT  NO:  NSWCDD/MP-94/363 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Original  contains  color  plates. 
All  DTIC/NTIS  reproductions  will  be  in  black  and  white. 

ABSTRACT:  (U)  The  Horizon  Infrared  Surveillance  Sensor 
(HISS)  Phase  2  system  was  involved  in  field  testing  at 
Wallops  Island,  Virginia  from  November  1993  through  April 
1994.  This  report  discusses  the  HISS  project  and  presents 
results  from  the  analysis  of  system  detection  range  perfor¬ 
mance.  The  HISS  Phase  2  detection  range  performance  has 
been  used  to  demonstrate  IR  contributions  to  an  integrated 
sensor  system  and  to  verify  detection  range  predictions  and 
improve  the  fidelity  of  current  detection  range  performance 
models. 
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(U)  Proceedings  of  the  Third  International  Symposium  on 
Long  Wavelength  Infrared  Detectors  and  Arrays:  Physics  and 
Applications  III  Held  in  Chicago,  Illinois  on  8-13  October 
1995,  Volume  95-28. 

APR  96  273P 

PERSONAL  AUTHORS:  Calvo,  Roque  J, 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  The  third  symposium  on  Long  Wavelength 
Infrared  Detectors  and  Arrays:  Physics  and  Applications  was 
held  at  the  Fall  Meeting  of  the  Electric  Chemical  Society, 
October  8-13,  1995  in  Chicago.  A  collection  of  papers  was 
published  in  the  beginning  of  the  Third  International  Sympo¬ 
sium. 
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DESCRIPTIVE  NOTE:  Journal  article 
OCT  95  lOP 

PERSONAL  AUTHORS:  Green,  Thomas  J,,  Jr.;  Marcus, 
Stephen;  Colella,  Barry  D. 

REPORT  NO:  MIT-JA-7186 
CONTRACT  NO:  F19628-95-C-0002 

UNCLASSIFIED  REPORT 

AVAILABILITY:  Pub.  in  Applied  Optics,  v34  n30 
p694 1-6949, 20  Oct  95.  Available  only  to  DTIC  users.  No 
copies  furnished  by  NTIS. 

ABSTRACT:  (U)  We  report  the  operation  of  an  imaging  Nd: 
YAG  microchip  laser  synthetic  aperture  radar,  with  which  we 
imaged  two  dimensional  (2-D)  models  of  military  targets.  The 
images  obtained  showed  spatial  resolution  significantly  better 
than  the  diffraction  limit  of  the  real  aperture  in  the 
along-track  dimension.  The  signal  processing  is  described, 
and  the  measurement  sensitivity  is  both  predicted  and 
verified.  In  addition,  2-D  images  with  high  resolution  in  both 
dimensions  were  generated  by  using  an  asymmetric  aperture 
to  match  the  along-track  synthetic  aperture  resolution  with  the 
across-track  diffraction  limited  resolution. 
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(U)  Developments  in  Infrared  Detection  and  Tracking 
Technology  and  Applications  to  Missiles. 

MAR  96  49P 

PERSONAL  AUTHORS:  He,  Liping 
REPORT  NO:  NAIC-ID(RS)T-0633-9 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  unidentified  Chinese 
language  periodical,  45p. 

ABSTRACT:  (U)  The  article  generally  describes  the  func¬ 
tions,  work  modes,  cooling  methods  and  major  performance 
indicators  of  infrared  detectors.  Detailed  presentations  are 
given  on  development  of  infrared  detection  tracking  technol¬ 
ogy,  in  addition  to  several  infrared  detectors  with  develop¬ 
ment  prospects  in  weapon  systems.  Moreover,  applications  of 
this  technology  are  enumerated  in  the  field  of  tactical 
missiles,  SDI  project,  and  in  the  Gulf  War.  Finally,  several 
proposals  are  made  on  future  developments. 
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(U)  Countermeasures  of  Military  Space  Systems. 

FEB96  18P 

PERSONAL  AUTHORS:  Bo,  Tian 
REPORT  NO:  NAIC-ID(RS)T-0620-95 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  unknown  Chinese 
source  np  nd, 

ABSTRACT:  (U)  Countermeasures  can  effectively  ensure  the 
safety  of  space  systems  in  coping  with  space  challenges  in 
military  space  systems.  Briefly,  the  article  presents  the 
features  and  difficulties  of  space  system  countermeasures, 
along  with  the  types  and  attack  principles  against  hostile 
threats  against  satellites,  protective  techniques  of 
satellite-borne  systems,  selection  of  countermeasure  plat¬ 
forms,  countermeasure  types,  countermeasure  techniques 
with  multiple  platforms,  space  C3  countermeasures,  space 
countermeasure  techniques,  and  integration  of  space  counter¬ 
measures. 
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(U)  The  Development  of  Background  Limited  Infrared 
Detection  Systems  Outside  of  China. 

FEB  96  28P 

PERSONAL  AUTHORS;  Hui,  Li;  Zibin,  Wang 
REPORT  NO:  NAIC-ID(RS)T-0631-95 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  Guo  Wai  Bei  Jing 
Xian  Hong  Wai  Tan  Ce  Xi  Tong  De  Fa  Zhan  (China) 
p460-470. 

ABSTRACT;  (U)  This  article  summarizes  the  development  of 
background  limited  infrared  detection  systems  outside  China, 
describing  a  number  of  infrared  space  detection  systems.  It 
primarily  introduces  long  wave  infrared  focal  plane  array 
technologies.  Finally,  it  puts  forward  our  proposals. 
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REPORT  NO:  NAIC-ID(RS)T-06 10-95 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  unknown  Chinese 
source  np  nd. 

ABSTRACT:  (U)  The  heat  model  of  laser  irradiated  PC  t5q)e 
HgCdTe  detectors  is  built  with  the  help  of  measuring  the 
resistance  temperature  relation  and  the  time  variation  when 
laser  irradiated.  The  laser  damage  thresholds  under  three 
damage  mechanisms  are  calculated. 
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Countermeasures  Technology:  Development  of  an  Active 
Aircraft  Defense  System  for  the  C/KC-135  Aircraft.  Volume 
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UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  Modem  Surface  to  Air  Missiles  (SAMs) 
present  a  significant  threat  to  today’s  military  and  civilian 
aircraft.  Current  countermeasure  systems  such  as  flares  and 
chaff  rely  on  decoying  the  missile  threat  and  do  not  provide 
adequate  protection  against  advanced  computerized  missiles 
(Schaffer,  1993:1).  An  aircraft  defense  system  that  actively 
seeks  out  and  defeats  an  incoming  missile  by  placing  a 
physical  barrier  in  the  missile’s  path  offers  a  promising 
alternative  to  current  countermeasures  technology.  This  thesis 
reports  the  preliminary  design  of  an  active  aircraft  defense 
system  for  the  protection  of  the  C/KC-135  aircraft  from 
SAMs.  The  developed  system  utilizes  a  kinetic  kill  mecha¬ 
nism  to  protect  the  aircraft  from  shoulder  launched  missiles 
while  the  aircraft  is  in  the  takeoff  and  climbout  configura¬ 
tions.  Both  smart  anti-missile  expendables  and  dumb  projec¬ 
tile  expendables  are  evaluated.  The  iterative  Systems  Engi¬ 
neering  approach  is  used  to  narrow  the  solution  set  to  the 
optimal  design.  The  final  outcome  is  the  refined  design  of  two 
candidate  aircraft  defense  systems  employing  a  kinetic  kill 
mechanism.  Both  systems  utilize  a  modified  ultra-violet 
tracker  and  employ  one  of  two  types  of  nets,  one  made  out  of 
Detonation  Cord  and  the  other  made  out  of  Spectra. 
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ABSTRACT:  (U)  Automatic  detection  of  time-critical  mobile 
targets  using  spectral-only  infrared  radiance  data  is  explored. 
A  quantification  of  the  probability  of  detection,  false  alarm 
rate,  and  total  error  rate  associated  with  this  detection  process 
is  provided.  A  set  of  classification  features  is  developed  for 
the  spectral  data,  and  these  features  are  utilized  in  a  Bayesian 
classifier  singly  and  in  combination  to  provide  target 
detection.  The  results  of  this  processing  are  presented  and 
sensitivity  of  the  class  separability  to  target  set,  target 
configuration,  diurnal  variations,  mean  contrast,  and  ambient 
temperature  estimation  errors  is  explored.  This  work 
introduces  the  concept  of  atmospheric  normalization  of 
classification  features,  in  which  feature  values  are  normalized 
using  an  estimate  of  the  ambient  temperature  surrounding  the 
target  being  observed  and  applying  the  Plauck  radiation  law 
with  those  estimates.  This  technique  is  demonstrated  to 
reduce  the  total  error  rate  associated  with  classification 
processing  to  less  than  one-fourth  of  that  observed  using 
non-normalized  features.  Classification  testing  of  spectral 
field  measurements  made  on  an  array  of  U.S.  and  foreign 
military  assets  reveals  a  total  error  rate  near  5%  with  a  95% 
probability  of  detection  and  a  concurrent  false  alarm  rate  of 
4%  when  a  single  classification  feature  is  employed.  Multiple 
feature  classification  on  the  same  data  yields  detection 
probabilities  near  97%  with  a  concurrent  false  alarm  rate  of 
2.5%,  Sensitivity  analysis  indicates  that  the  probability  of 
detection  is  reduced  to  70-75%  in  the  hours  preceding 
daylight,  and  that  for  the  total  error  rate  to  be  less  than  10%, 
the  target-to-background  mean  contrast  must  be  greater  than 
0. 
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FUNCTIONS,  TARGET  RECOGNITION,  PROBABILITY, 
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IDENTIFIERS:  (U)  C-135  AIRCRAFT,  KC-135  AIRCRAFT. 


The  DTIC  Review 


Defense  Technical  Information  Center 


AD-A304  832 

NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA 

(U)  Existing  and  Emerging  Maritime  Surveillance 
Technologies. 

DESCRIPTIVE  NOTE:  Master’s  thesis 
SEP  95  135P 

PERSONAL  AUTHORS:  Burkholder,  Gary  A. 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  This  paper  presents  an  assessment  of 
present  and  near-future  maritime  surveillance  technologies. 
Radar,  optical  and  infrared  sensor  systems  are  each  discussed. 
Thirty-two  major  domestic  and  foreign  manufacturers’  market 
position,  current  products  and  their  applications  are  reviewed. 
State-of-the-art  research  is  presented  and  analyzed  to  see  in 
what  direction  the  future  of  maritime  surveillance  is  headed 
In  examining  current  systems  and  those  being  researched, 
rapid  signal  processing,  data  fusion,  interpretation  and 
dissemination  constitute  the  major  challenges  and  provide  the 
most  opportunities  for  improvment. 

DESCRIPTORS:  (U)  *AERIAL  RECONNAISSANCE, 
*RADAR,  ^INFRARED  DETECTORS,  ^OPTICAL 
DETECTORS,  *OCEAN  SURVEILLANCE,  SIGNAL 
PROCESSING,  COASTAL  REGIONS, 

MANUFACTURING,  MARKETING,  REALTIME,  STATE 
OF  THE  ART,  ACOUSTIC  WAVES,  SIGNAL  TO  NOISE 
RATIO,  DOPPLER  RADAR,  THESES,  INTEGRATED 
CIRCUITS,  SYNTHETIC  APERTURE  RADAR, 
SEMICONDUCTORS,  AIR  TRAFFIC  CONTROL 
SYSTEMS,  FIBER  OPTICS  TRANSMISSION  LINES, 
CAMERAS,  PHOTOGRAPHY,  ELECTROMAGNETIC 
INTERFERENCE,  TRAVELING  WAVE  TUBES,  DATA 
FUSION,  RADAR  CLUTTER,  RADAR  ANTENNAS, 
ANALOG  TO  DIGITAL  CONVERTERS,  LAW 
ENFORCEMENT,  TRANSMITTER  RECEIVERS, 
INFRARED  TRACKING,  EARLY  WARNING  SYSTEMS, 
KLYSTRONS,  MAGNETRONS. 
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(U)  Electronic  Warfare  in  China's  Past,  Present,  and  Future. 

NOV  95  12P 

PERSONAL  AUTHORS:  Xu,  Shaoxing 
REPORT  NO:  NAIC-ID(RS)T-0273-95 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  unidentified  Chinese 
language  document,  n.d. 

ABSTRACT:  (U)  Chinese  electronic  warfare  has  a  recorded 
beginning  in  February  of  1 93 1 .  Before  that  war,  personnel  on 
our  side  maintained  intercept  of  enemy  stations  for  more  than 
twenty  days,  grasping  the  enemy  points  of  departure,  times, 
routes,  and  causing  our  forces  to  achieve  a  very  brilliant 
victory.  Praise  and  commendation  were  received  from 
Committee  Member  Mao  and  Commander  in  Chief  Zhu. 
During  the  more  than  thirty  years  of  history  from  1931  to 
1965,  this  is  the  breeding  period  of  Chinese  electronic 
warfare.  During  the  entire  phase,  it  is  possible  to  make  use  of 
the  six  characters’  analysis,  renovation,  system  (illegible)  in 
order  to  summarize  the  development  of  the  specialty  in 
question.  There  were  no  electronic  countermeasures 
personnel  in  the  early  period.  With  study  and  analysis  on  the 
one  hand  and  renovation  of  captured  Japanese  and  U.S. 
intercept,  reconnaissance,  and  jamming  equipment  on  the 
other— eager  to  meet  the  needs  of  war— this  equipment  was  put 
to  new  uses  in  battle.  In  application,  there  was  familiarization, 
mastery,  and  study.  In  particular,  during  the  ten  or  so  years 
after  liberation,  stress  was  placed  on  the  analysis  of  electronic 
countermeasures  equipment  on  the  U.S.  P2V7  (illegible)  low 
altitude  reconnaissance  aircraft  and  the  U-2  high  altitude 
reconnaissance  aircraft  as  well  as  the  Soviet  built  cp6-l  radar 
reconnaissance  machine  and  cp6-5  radar  jammer, 
communications  reconnaissance,  jamming  devices,  pulse 
analyzers,  and  so  on.  At  the  same  time  as  this,  a  number  of 
pieces  of  airborne,  ground,  and  submarine  borne 
reconnaissance  and  jamming  equipment  were  manufactured 
in  imitation. 
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(U)  Enhanced  Capabilities  of  Advanced  Airborne  Radar 
Simulation. 

DESCRIPTIVE  NOTE:  Final  technical  rept.  Oct  93-Sep  94 
JAN  96  117P 

PERSONAL  AUTHORS:  Babu,  V.  N.;  Tortes,  J.  A. 
CONTRACT  NO:  F19628-94-C-0001 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  MITRE  has  enhanced  the  capabilities  of  the 
Advanced  Airborne  Radar  Simulation,  a  flexible  end-to-end 
steady-state  simulation  that  can  be  used  to  evaluate  the 
critical  issues  affecting  the  performance  of  advanced  airborne 
radars  that  use  adaptive  antenna  techniques,  i.e.,  space-time 
processing  (STP),  for  both  jamming  and  clutter  suppression. 
For  a  large  phased-array  antenna,  the  limitations  of  a  practical 
implementation  of  optimal  STP  will  generally  require 
suboptimal  adaptive  array  architectures,  such  as  element 
space,  beam  space,  and  subarrays.  The  simulation  was 
modified  to  augment  these  suboptimal  architectures  with 
spatial  degrees  of  freedom  in  elevation.  The  received  signal 
models,  i.e.,  target  jamming  and  clutter,  were  modified  to 
include  the  effects  of  near-field  scattering  using  a  multiple 
point  scatterer  model. 

DESCRIPTORS:  (U)  ^SIMULATION,  *RADAR,  STEADY 
STATE,  SPATIAL  DISTRIBUTION,  SCATTERING, 
MODELS,  AIRBORNE,  ARRAYS,  DEGREES  OF 
FREEDOM,  SIGNALS,  ANTENNAS,  ADAPTIVE 
SYSTEMS,  JAMMING,  RADAR  CLUTTER,  ELEVATION, 
ARCHITECTURE,  NEAR  FIELD,  SUPPRESSION, 
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(U)  Utility  Analysis  of  High  Resolution. 

DESCRIPTIVE  NOTE:  Final  rept.  26  Jan-1  Oct  94 
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PERSONAL  AUTHORS:  Gerber,  Garth  L.;  Reiley,  Michael 
F;  Eismann,  Michael  T.;  Jacobs,  James  M. 

REPORT  NO:  ERIM-253885-2-F-VOL-1 
CONTRACT  NO:  DLA900-88-D-0392 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  ADA303475  ADB206544 

ABSTRACT:  (U)  This  report  summarizes  the  results  of  two 
studies  conducted  to  address  specific  issues  related  to  the 
utility  of  a  specific  high  resolution  multispectral  sensor, 
designated  the  F  sensor.  The  first  study  assessed  the  validity 
of  a  sensor  model,  entitled  the  Performance  Prediction  Model 
(PPM),  which  was  provided  by  the  sensor  contractor.  This 
was  performed  by  cross  checking  results  with  an  independent 
sensor  model.  In  most  cases,  good  agreement  was  achieved. 
The  second  study  assessed  the  expected  performance  of  a 
vital  component  of  the  F  sensor  a  magnetic  bearing  system 
(MBS)  for  line  of  sight  control  and  stabilization. 
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DESCRIPTIVE  NOTE:  Final  rept.  26  Jan  94-31  May  95 
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PERSONAL  AUTHORS:  Eismann,  Michael  T;  Ingle, 
Stephen  D. 

REPORT  NO:  ERIM-253885-2-F-VOL-3 
CONTRACT  NO:  DLA900-88-D-0392 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  ADB  206544,  ADA  303476 

ABSTRACT:  (U)  The  Image  Based  Sensor  Model  (IBSM)  is 
a  modular  set  of  numerical  tools  for  designing,  evaluating, 
and  modeling  electrooptical  and  infrared  (EOJIR)  imaging 
sensors.  The  primary  motivation  which  led  to  the 
development  of  IBSM  was  the  need  for  a  model  which  (a) 
could  produce  simulated  sensor  imagery  (based  on  high 
fidelity  input  imagery)  in  addition  to  sensor  performance 
metrics  to  better  characterize  the  imaging  performance  of  a 
sensor  system,  and  (b)  provides  the  flexibility  to  evaluate  and 
compare  sensors  and  imaging  configurations  with  differing 
characteristics  without  rewriting  computer  code.  The  model 
operates  within  the  Khoros  Cantata  environment,  and  can 
perform  realistic  simulation  of  image  degradations  and 
parametric  modeling  with  a  wide  range  of  atmospheric, 
sensor,  and  processing  effects.  This  report  provides  a 
comprehensive  overview  of  the  model. 
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PERSONAL  AUTHORS:  Kulp,  T.  U.;  McRae,  T,  G.; 
Kennedy,  R.;  Garvis,  D. 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  This  technical  report  describes  the 
methodology  to  build  and  operate  a  remotely  operated,  radio 
controlled,  battery  powered  infrared  detector  for  measuring 
hydrogen  chloride  (HCL)  emissions  from  Space  Launch 
operations. 
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CONTRACT  NO:  DAAL03-92-G-0355 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  Adaptive  FSK/PSK  is  a  spectrally  agile 
pulse  corpression  radar  signaling  scheme  that  learns 
target-specific  information  over  time  and  integrates  that 
information  into  the  design  of  its  radar  signal  processor  in 
real-time.  This  use  of  target-specific  information  into  the 
real-time  transmitter  signal  design  process  enhances  the 
return  from  the  target,  thus  improving  the  performance  of  the 
radar  system  in  both  additive  noise  environments  (jamming) 
and  ground  clutter  limited  environments  (low  altitude 
targets).  Since  FSK/PSK  signaling  is  a  pseudorandom  pulse 
compression  scheme  with  a  large  time-bandwidth  product,  it 
has  inherent  low  probability  of  intercept  characteristics. 
Additionally,  adaptive  FSK/PSK  signaling  is  able  to  earn 
target  specific  information  which  can  be  used  to  identify  the 
target  using  existing  target  identification  techniques.  In  this 
study,  the  cross-ambiguity  function  of  a  general  FSK/PSK 
waveform  is  derived  and  analyzed  in  a  probabilistic  manner. 
The  results  of  this  probabilistic  analysis  are  used  to  develop  a 
design  methodology  that  allows  the  signal  to  be  matched  to 
an  arbitrary  target.  It  is  shown  that  both  the  signal  to  clutter 
ratio  and  signal  to  noise  ratio  of  a  matched  FSK/PSK  radar 
are  superior  to  those  of  a  conventional  signaling  radar  that 
utilizes  no  target-specific  information. 
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REPORT  NO:  NAIC-ID(RS)T-03  80-95 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  Yu  Jing  Ji  Dui  Kang 
Chu  Jan  (China)  p33-35,  n.d. 

ABSTRACT:  (U)  This  article  makes  a  simple  introduction 
and  explanation  with  regard  to  the  roles  of  early  warning 
aircraft  in  modem  warfare  as  well  as  their  electronic 
equipment.  It  points  out  that  countermeasures  concerning 
early  warning  aircraft  will  be  combined  countermeasures. 
They  include  destruction  of  the  fuselage  of  early  warning 
aircraft  as  well  as  the  use  of  electronic  countermeasure 
equipment  to  carry  out  jamming.  In  order  to  achieve  effective 
countermeasures  against  early  warning  aircraft,  the  author 
puts  forward  research  work  associated  with  8  areas  which 
should  be  developed,  such  as,  active/passive  combined 
detection  systems. 
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PERSONAL  AUTHORS:  Baichuan,  Liang 
REPORT  NO:  NAIC-ID(RS)T-0384-95 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  unknown  Chinese 
source. 

ABSTRACT:  (U)  Based  on  his  research  of  the  characteristics 
of  different  electronic  systems  of  early  warning  aircraft,  the 
author  proposes  methods  of  countering  early  warning  aircraft, 
and  studies  techniques  of  jamming  their  radar, 
communications  systems,  and  Identification  Friend  or  Foe 
IFF  systems. 
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SUPPLEMENTARY  NOTE:  Trans,  of  Yi  Zhong  Zi  Shi  De  Ji 
Zai  Wu  Yuan  Guang  Dlan  Dui  Kang  Xi  Tong  (China)  PI 5-26. 

ABSTRACT:  (U)  Passive  jamming  and  infrared  jamming  are 
among  important  component  parts  making  up  combined 
electronic  countermeasures  systems  at  the  present  time.  In 
modem  warfare,  following  along  with  the  widespread 
application  of  precision  guided  weapons— in  particular, 
photoelectric-infrared  guidance  weapons— they  constitute  a 
threat  to  the  survival  of  combat  aircraft  which  grows  more 
severe  by  the  day.  Modem  limited  wars— in  particular,  the 
practical  realization  of  the  Gulf  war— clearly  show  that  most 
damaged  or  lost  aircraft  were  all  due  to  being  attacked  by 
infrared  homing  missiles.  At  the  present  time,  approximately 
80%  of  missiles  opt  for  the  use  of  infrared  guidance. 
Moreover,  the  U.S.  has  already  taken  passive  photoelectric 
detection  equipment  to  act  as  the  principal  military  sensors  of 
the  future.  It  has  also  been  used  to  supply  target  acquisition, 
tracking,  and  positioning  capabilities  for  stealth  aircraft  in 
order  to  strengthen  covering  capabilities.  As  a  result,  the 
militaries  of  various  nations  pay  serious  attention  to 
passive-photoelectric  countermeasure  system  technology 
equipment,  which  has  already  become  one  type  of  effective 
means  for  modem  combat  aircraft  to  carry  out  defense 
penetrations  and  to  increase  effective  guarantees  of  airborne 
security. 
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UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  During  November  1994,  the  USCO 
conducted  a  field  evaluation  of  several  commercially 
available  portable  infrared  (IR)  sensors  for  their  ability  to 
detect  oil  on  water  at  night.  The  evaluation  was  conducted 
over  the  natural  oil  seeps  off  the  coast  of  Santa  Barbara,  CA. 
The  portable  sensors  were  operated  from  the  open  door  of  a 
Coast  Guard  HH-60J  helicopter.  Sensors  evaluated  included 
the  Agema  Thermovision  1000  and  Texas  Instruments 
LOCUSP  operating  in  the  long  wave  IR  (LWIR)  and  the 
Cincinnati  Electronics  IRC-160ST  and  IRRIS-256ST 
operating  in  the  medium  wave  IR  (MWIR).  The  installed 
FLIR  2000  LWIR  system  provided  a  baseline  reference  of 
current  CG  IR  capabilities.  The  LOCUSP  was  the  only 
uncooled  sensor  evaluated. 
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UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  The  Infrared  Astronomy  Satellite  (IRAS) 
was  launched  and  operated  during  a  10-month  period  in  1983. 
The  orbit  was  close  to  the  MSX  orbit  design,  but  the  science 
data  were  collected  in  a  mode  where  the  focal  plane  was 
pointing  directly  away  from  the  Earth.  The  Space  Research 
Institute  at  Groningen,  Netherlands,  collected  approximately 
139,000  tracks  of  data  that  had  focal  plane  motion  different 
than  astronomical  sources,  and  the  IRAS  Processing  and 
Analysis  Center  (IPAC)  determined  the  boresite  pointing  of 
IRAS  to  within  20  arcsec.  This  report  will  focus  on  the 
nonastronomical  detections  from  IRAS,  many  of  which  are 
Resident  Space  Objects  (RSOs).  In  particular,  the  focus  of  the 
study  is  on  how  many  are  correlated  to  the  known  RSO 
catalogue  for  1983  and  the  calibration  and  characterization  of 
the  metric  accuracy  for  the  correlated  data.  This  study  was 
undertaken  to  prepare  for  analysis  of  RSO  detections  from  the 
MSX  satellite,  and  in  particular,  so  that  automatic  analysis 
tools  designed  for  analysis  of  surveillance  experiment  data 
could  be  tested.  The  supporting  analysis  tools,  required 
corollary  data,  and  metric  calibration  procedure  will  be 
described,  and  the  results  of  the  accuracy  of  the  IRAS 
ephemeris  and  metric  RSO  detections  will  be  presented. 
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SUPPLEMENTARY  NOTE:  Prepared  in  cooperation  with 
Syracuse  Univ,,  NY  ADA300898  ADA300900  ADA300901 
ADA300902 

ABSTRACT:  (U)  This  thesis  addresses  the  performance  of 
Locally  Optimum  Detectors  in  radar  weak  signal  detection  for 
finite  sample  sizes  where  the  radar  disturbance  is  modeled  as 
a  correlated  non-Gaussian  random  process.  The  theory  of 
Spherically  Invariant  Random  Process  is  used  for  statistical 
characterization  of  non-Gaussian  radar  clutter.  In  particular, 
the  K-distribution  and  the  Student-T  distributions  have  been 
considered  as  models  for  radar  clutter.  A  canonical  form  is 
established  for  the  Locally  Optimum  Detector  that  is  a 
product  of  the  Gaussian  linear  receiver  and  a  zero  memory 
nonlinearity.  The  functional  form  of  the  zero  memory 
nonlinearity  depends  on  the  approximation  used  for  the 
underlying  radar  clutter  probability  density  function.  Since 
the  weak  signal  detector  is  nonlinear,  thresholds  for  specified 
false  alarm  probability  cannot  be  established  in  closed  form. 
Given  a  specified  false  alarm  probability  a  new  method  for 
threshold  estimation  based  on  extreme  value  theory  is  derived 
that  reduces  by  orders  of  magnitude  the  computation  and 
sample  size  required  to  set  the  threshold.  Once  the  threshold 
is  set  the  performance  of  the  Locally  Optimum  Detector  is 
carried  out  for  finite  sample  sizes  through  computer 
simulations. 
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SUPPLEMENTARY  NOTE:  Trans,  of  Zhong  Bo  Duan  Hong 
Wai  Cheng  Xiang  Xi  Tong  Ji  Qi  (China)  pi -9,  11  Sep  95. 

ABSTRACT:  (U)  This  article  makes  only  a  simple 
introduction  of  the  status  of  developments  and  applications 
associated  with  3-5  micron  wave  band  infrared  imaging 
technology  outside  China.  Concise  analysis  and  comparisons 
are  carried  out  on  its  performance,  discussing  the  situations 
and  causes  giving  precedence  to  an  option  for  the  use  of  3  -  5 
micron  infrared  imagery. 
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UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  Hong  Wai  Tu  Xiang 
Shi  Bie  Chu  Li,  (China)  nl  p50-54  1994. 

ABSTRACT:  (U)  This  document  addresses  in  simple  terms 
the  flow  of  information  in  infrared  imagery  recognition  and 
processing,  and  proposes  using  knowledge-based  expert 
systems  and  advanced  neural  networks  to  bring  about  a 
development  direction  for  automated,  intelligent  imagery 
recognition  and  processing.  This  document  also  presents  a 
kind  of  effective  missile-borne  infrared  imagery  detection, 
tracking,  and  extrapolation  method. 
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DESCRIPTIVE  NOTE:  Professional  paper 
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UNCLASSIFIED  REPORT 

AVAILABILITY:  Pub.  in  Sensor  and  Propagation  Panel 
Symposium  on  Propagation  Assessment  in  Coastal 
Environments,  pi 8-1  - 18-10, 19-22  Sep  94.  Available  only  to 
DTIC  users.  No  copies  furnished  by  NTIS. 

ABSTRACT:  (U)  This  paper  describes  a  1993  field 
experiment  entitled  Marine  Aerosol  Properties  and  Thermal 
Imager  Performance  Trial  (MAPTIP)  conducted  by  NATO 
AC/  243  Panel  04/RSG.8  and  04/RSG.5  in  the  Dutch  coastal 
waters.  Objectives  were:  to  improve  and  validate  vertical 
marine  aerosol  models  by  providing  an  extensive  set  of 
aerosol  and  meteorological  measurements,  within  a  coastal 
environment  at  different  altitudes  and  for  a  range  of 
meteorological  conditions;  make  aerosol  and  meteorological 
observations  in  the  first  10  m  of  the  ocean  surface  with  a  view 
to  extending  existing  aerosol  models  to  incorporate 
near-surface  effects;  and  to  assess  marine  boundary  layer 
effects  on  thermal  Imaging  systems. 
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UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Trans,  of  unknown  source 
(China),  p45-47,  n.d. 

ABSTRACT:  (U)  Up  to  the  present  time,  the  majority  of  laser 
radars  which  have  been  reported  were  very  expensive,  had 
large  volumes,  and  were  very  complicated  in  terms  of 
utilization  and  maintenance.  As  far  as  this  type  of  low  cost 
laser  range  finding  and  detection  system  reported  in  this 
paper  is  concerned,  use  is  made  of  a  type  of  semiconductor 
laser  diode  emission  device  and  a  small  model  computer  to  do 
control.  With  regard  to  making  use  of  semiconductor  laser 
diodes  to  replace  Nd:  YAG  lasers,  the  results  were  to  improve 
optical  emission  efficiencies,  electronic  control  compatibility, 
and  reduce  system  dimensions.  One  type  of  high  sensitivity 
PIN  diode  receiver  is  capable  of  receiving  radar  pulses 
reflected  back  from  targets,  and,  in  conjunction  with  that, 
directly  supplying  digital  output.  Electronic  methods  are  used 
to  measure  laser  pulse  travel  time  of  flight.  In  conjunction 
with  that,  an  IBM- AT  compatible  minicomputer  is  used,  and 
it  is  possible  to  realize  the  needed  calculations  as  well  as 
displaying  target  distance,  speed,  and  position.  Systems  are 
capable  of  tracking  horizon  range  targets  more  than  50m 
away  from  emitters. 
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SUPPLEMENTARY  NOTE:  Trans,  of  Bulletin  of  Nanjing 
University  of  Science  and  Engineering  (China)  n7(77),  1994. 

ABSTRACT:  (U)  The  effectiveness  of  infrared  guidance  and 
radar  guidance  jamming  with  infrared  bait  and  chaff  has  been 
proven  in  modem  high-tech  warfare.  Their  jamming  regime 
relies  on  infrared  radiation  of  the  infrared  bait  and  radar 
reflective  waves  from  chaff  to  simulate  the  target  (such  as 
aircraft,  warships,  tanks,  and  point  targets)  performance  to 
confuse,  disturb,  and  jam  the  infrared  guidance  and  radar 
guidance,  as  well  as  their  probing  and  aiming  systems.  In 
tactical  applications,  these  are  divided  into  centroid  mode, 
dump  mode,  dilution  mode,  and  confusion  mode.  Study  of  the 
jamming  performance  of  infrared  bait  and  chaff  is  required  in 
developing  the  smoke,  flame,  light,  and  electric/sourceless 
jamming  technologies.  It  is  very  important  and  practical  to 
enhance  the  jamming  performance  of  infrared  bait  and  chaff, 
as  well  as  to  extend  the  jamming  frequency  spectrum. 
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ABSTRACT:  (U)  The  study  was  to  investigate  in-cylinder 
events  of  a  direct  injection-type  diesel  engine  by  using  a  new 
high-speed  infrared  (IR)  digital  imaging  system  for  obtaining 
information  that  was  difficult  to  achieve  by  the  conventional 
devices.  For  this,  a  new  highspeed  dual-spectra  infrared 
digital  imaging  system  was  developed  to  simultaneously 
capture  two  geometrically  identical  (in  respective  spectral) 
sets  of  IR  images  having  discrete  digital  information  in  a 
(64x64)  matrix  at  rates  as  high  as  over  1 ,800  frames/sec  each 
with  exposure  period  as  short  as  20  micron  sec.  At  the  same 
time,  a  new  advanced  four-color  W  imaging  system  was 
constructed.  The  first  two  sets  of  spectral  data  were  the 
radiation  from  water  vapor  emission  bands  to  compute  the 
distributions  of  temperature  and  specie  in  the  gaseous  mixture 
and  the  remaining  two  sets  of  data  were  to  find  the 
instantaneous  temperature  distribution  over  the  cylinder 
surface.  More  than  eight  reviewed  publications  have  been 
produced  to  report  many  new  findings  including: 

Distributions  of  Water  Vapor  and  Temperature  in  a  Flame; 

End  Gas  Images  Prior  to  Onset  of  Knock;  Effect  of  MTBE  on 
Diesel  Combustion;  Impact  of  Oxygen  Enrichment  on 
In-cylinder  Reactions;  Spectral  IR  Images  of  Spray  Plume; 
Residual  Gas  Distribution;  Preflame  Reactions  in  Diesel 
Combustion;  Preflame  Reactions  in  the  End  Gas  of  an  SI 
Engine;  Postflame  Oxidation;  and  Liquid  Fuel  Layers  during 
Combustion  in  an  SI  Engine.  In  addition,  some  computational 
analysis  of  diesel  combustion  was  performed  using  KJVA-II 
program  in  order  to  compare  results  from  the  prediction  and 
the  measurements  made  using  the  new  IR  imaging  diagnostic 
tool, 
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ABSTRACT:  (U)  This  report  compares  measured  forward- 
looking  infrared  (FLIR)  system  detection  ranges  of  a  target 
ship,  the  Research  Vessel  Point  Sur,  with  predictions  from  the 
Electro-Optical  Tactical  Decision  Aid  (EOTDA)  version  3.0. 
The  EOTDA  was  primarily  developed  by  the  Air  Force,  with 
only  minimal  attention  applied  to  the  complexities  of  the 
marine  environment.  The  objective  of  this  case  study  was  to 
evaluate  the  EOTDA  performance  for  a  target  in  an  open 
ocean  background. 
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ABSTRACT:  (U)  Rapid  advances  have  been  made  recently  in 
the  development  and  applications  of  infrared  thermal  imaging 
technology  used  primarily  for  military  purposes.  In  this 
article,  advances  in  infrared  thermal  imaging  are  introduced 
in  brief,  and  discussion  is  focused  on  technical  progress  in  the 
field  of  infrared  focal  plane  arrays  (IRFPA).  Also,  the  trend 
for  infrared  thermal  imaging  technique  to  expand  rapidly 
from  purely  military  use  to  civil  applications  is  demonstrated. 
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(U)  Focal  Plane  Characterization  Capabilities  at  NCCOSC 
RDT&E  Division. 

DESCRIPTIVE  NOTE:  Professional  paper 
JUN95  9P 

PERSONAL  AUTHORS:  Frisbie,  S.  A. 

UNCLASSIFIED  REPORT 

AVAILABILITY:  Pub.  in  Smart  Focal  Plane  Arrays  and  Focal 
Plane  Array  Testing,  v2474,  pl87-192,  Apr  95.  Available  only 
to  DTIC  users.  No  copies  furnished  by  NTIS. 

ABSTRACT:  (U)  The  infrared  measurement  facility  at  the 
NRaD  (formerly  NOSC)  division  of  NCCOSC  has  been 
involved  in  infrared  radiometric  testing  for  over  forty  years. 
This  versatile  facility  provides  complete  radiometric  testing 
of  single  element  detectors,  detector  arrays  and  hybrid  focal 
plane  arrays  (FPAs).  Measurement  capabilities  include 
spectral  response  flood  illumination,  spot  illumination  and 
ionizing  radiation  testing  in  low  to  medium  background 
photon  flux  environments. 
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PERSONAL  AUTHORS:  Dityatyev,  A.  A.;  Gurey,  A.  E.; 
Mazaev,  A.  A.;  Semenov,  O.  G.;  Tikhomirov,  A.  A. 

REPORT  NO:  SLPI-1/95 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  On  the  base  of  analysis  of  experimental 
results  and  theoretic  researches  the  phenomenological  model 
of  a  charge  carrier  generation  under  the  action  of  irradiation  is 
presented,  and  the  carrier  transport  in  thin  film  photoadvices 
is  discussed.  A  possibility  of  IR  sensible  devices  realization 
on  the  base  of  high  fullerenes  or  doped  fullerenes  is 
considered.  It  is  proposed  to  study  Shottky  transition  and 
hetero  Junction  with  polymers  for  realization  of  all 
photosensitive  properties  of  fullerenes.  On  the  base  of 
extended  Huckel  model  the  numeric  calculations  of  electron 
structure  of  high  fullerenes  and  the  search  of  isomers  with 
defined  band  gap  value  are  made.  Experimental  researches  of 
absorption  spectra  for  UV-,  optic  and  IR  bands  are  made, 
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(U)  Giant  Enhanced  Backscattering  in  Laser  Radar 
Signatures. 

1995  13P 

PERSONAL  AUTHORS:  Gu,  Zu-Han;  Dummer,  Richard  S.; 
McKenna,  Paul;  Maradudin,  Alexei  A.;  Estep,  Jeff  A. 
CONTRACT  NO:  DAAH04-93-G-0014 

UNCLASSIFIED  REPORT 

AVAILABILITY:  Pub.  in  SPIE,  v2271,  p97-108,  1995. 
Available  only  to  DTIC  users.  No  copies  furnished  by  NTIS. 

ABSTRACT:  (U)  Due  to  underestimation  of  rough  surface 
scattering  at  retro-reflection  direction,  consideration  of 
enhanced  backscattering  (i.e.,  opposition  effect)  can  be 
critically  important  in  developing  laser  radar  signatures.  This 
is  illustrated  by  the  comparison  of  two  signature  models  of 
the  same  airplane  -  one  developed  using  bidirectional 
reflectance  measurement  data  at  one  degree  from  the  retro- 
direction,  and  the  other  developed  from  monostatic  measure¬ 
ment  data  at  true  retro-direction.  The  unique  instrument  used 
to  obtain  the  laser  reflectance  data,  the  reflectance  data  itself, 
and  the  software  ensemble  used  to  develop  the  signature  are 
described  and  discussed. 
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B,;  Silverman,  Harvey  F.;  Blane,  Michael  M. 

REPORT  NO:  LEMS-TR-126 
CONTRACT  NO:  F49620-93- 1-0501 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  Our  original  research  on  the 
multicovariance  matched  filter  deals  with  optimum  low 
resolution  target  detection  in  a  single-frame  multicolor  image, 
such  as  a  multispectral  infrared  or  polarimetric  synthetic 
aperture  radar  picture.  The  multicovariance  method  com¬ 
pletely  uses  all  the  joint  variability  of  the  problem,  in  both 
space  and  frequency,  in  a  way  that  generalizes  both  the 
traditional  spatial  matched  filter  and  also  techniques  involv¬ 
ing  scalar  ratios  between  frequency  bands.  The  main  new 
focus  of  our  work,  directed  toward  achieving  the  best  target 
detection  performance  that  is  possible,  is  to  develop  a 
preprocessing  step  involving  optimal  adaptive  estimation  of 
the  local  clutter  background.  This  involves  segmenting  the 
image  into  regions,  which  correspond  to  different  back¬ 
ground/clutter  statistical  models.  Statistics  of  real  data  are 
being  studied  and  used  in  new  state-of-the-art  hierarchical 
segmentation  algorithms  based  on  Markov  Random  Field, 
polynomial  and  autoregressive  models  for  vector- valued 
random  processes.  The  major  algorithmic  challenges  here  are 
in  estimating  the  best  possible  background/clutter  models  and 
in  accurately  estimating  the  boundaries  between  different 
model  regions.  We  are  in  the  process  of  developing  extremely 
efficient  and  robust  algorithms  to  estimate  these  clutter 
models.  These  are  similar  to  familiar  algorithms  from 
mainstream  signal  processing,  but  solve  the  interpolation 
problem  for  Markov  Random  Fields,  which  is  different  than 
the  usual  linear  prediction  problem. 
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REPORT  NO:  NRL/MR/7220-95-7652 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Original  contains  color  plates. 
All  DTIC/NTIS  reproductions  will  be  in  black  and  white. 

ABSTRACT:  (U)  Many  military  and  environmental 
applications  require  wide  area  surveillance  on  a  timely  and 
global  basis.  Passive  microwave  and  millimeter  wave 
imaging  has  the  advantage  over  visual  and  infrared  systems 
of  being  a  near  all-weather  surveillance  system  with  the 
ability  to  see  through  clouds  and  most  adverse  weather 
conditions  with  equal  day  and  night  sensitivity.  The  major 
disadvantage  of  microwave  and  millimeter  wave  radiometric 
systems  is  that  the  relatively  long  wavelengths  require 
aperture  several  orders  of  magnitude  larger  than  at  infrared  to 
obtain  comparable  spatial  resolution.  However,  even  low 
resolution  microwave  and  millimeter  wave  radiometers  can 
contribute  significantly  to  intelligence  gathering  and 
environmental  remote  sensing.  This  document  presents  a 
study  of  the  capabilities  of  microwave  and  millimeter  wave 
imaging  from  space  with  systems  which  provide  resolutions 
of  100-400  m  at  observation  frequencies  of  1.4-140  GHz. 
Possible  radiometer  architectures  and  their  space 
development  are  examined.  The  surveillance  and 
environmental  remote  sensing  benefits  of  multifrequency 
high  resolution  microwave  and  millimeter  wave  radiometry 
are  discussed  in  depth. 
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(U)  Physics-Based  Infrared  Terrain  Radiance  Texture  Model. 

DESCRIPTIVE  NOTE:  Final  rept. 

FEB  95  123P 

PERSONAL  AUTHORS:  Weiss,  Richard  A.;  Sabol,  Bruce 
M.;  Smith,  James  A.;  Bales,  John  W. 

REPORT  NO:  WES/TR/EL-95-5 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Prepared  in  collaboration  with 
Tuskegee  Univ.,  AL. 

ABSTRACT:  (U)  This  report  develops  a  physics-based 
semiempirical  model  for  predicting  the  texture  of  an  infrared 
(IR)  image  of  terrain  areas  including  soil,  grass,  and  forest 
canopies.  The  model  is  based  on  an  empirical  determination 
of  the  relationship  between  the  standard  deviation  and 
correlation  length  of  the  terrain  radiance  variations  in  an  IR 
image  and  the  average  value  of  the  radiance  in  the  image.  The 
average  value  of  the  terrain  radiance  is  then  theoretically 
predicted  by  the  WESTHERM  terrain  radiance  prediction 
model  in  terms  of  local  terrain  and  weather  conditions  and  in 
terms  of  the  time  of  day.  In  this  way,  the  diurnal  variation  of 
the  standard  deviation  and  correlation  length  of  the  radiance 
or  brightness  temperature  can  be  predicted  for  specified 
terrain  and  weather  conditions. 
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PERSONAL  AUTHORS:  Shaffer,  William  A,;  Schaum,  Alan 
R;  Hinnrichs,  Michelle;  Massie,  Mark  A.;  Fields,  David 
REPORT  NO:  NRL/MR/5620--95-7648 

UNCLASSIFIED  REPORT 

SUPPLEMENTARY  NOTE:  Prepared  in  collaboration  with 
Pacific  Advanced  Technology,  Solvang,  CA  and  Lawrence 
Livermore  National  Lab.,  CA. 

ABSTRACT:  (U)  From  September  12-15, 1994,  two 
midwave  hyperspectral  imagers  (the  Infrared  Multi-spectral 
Sensor  (IMSS)  of  Pacific  Advanced  Technology  and  the 
Livermore  Imaging  Fourier  Transform  Infrared  Spectrometer 
(LIFTIRS)  of  Lawrence  Livermore  National  Laboratory) 
made  measurements  of  particular  patches  of  forest 
background  and  a  Lockheed  Santa  Cruz  Facility  test  stand  at 
various  times  of  day  and  observed  the  firing  of  the  rocket 
engine  from  a  sensitive  piece  of  hardware  mounted  in  the 
Lockheed  test  stand.  The  hyperspectral  image  cubes  of  the 
same  patch  of  background  observed  at  various  times  will  be 
used  to  investigate  the  temporal  nature  of  the 
constituent-temperature  interaction.  The  data  of  the  rocket 
firing  was  collected  as  part  of  a  classified  program  of  the  U.S. 
Army  Missile  Command.  Two  sets  of  data  cubes  taken  of  the 
same  background  region  by  IMSS  and  LIFTIRS 
simultaneously  will  be  used  both  to  compare  the  two 
instruments  and  to  aid  in  instrument  characterization.  This 
paper  describes  the  instruments  used  in  this  data  collection, 
the  types  of  backgrounds  observed  and  the  experiments 
performed,  and  the  data  collected. 
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(U)  Inductively  Coupled  IR  Staring  Array  Based  on 
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DESCRIPTIVE  NOTE:  Final  rept.  15  Jul  94-14  Jan  95 
MAR  95  29P 

PERSONAL  AUTHORS:  Hibbs,  Andrew  D.;  Singsaas,  Alan 
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REPORT  NO:  CM1198F 
CONTRACT  NO:  F49620-94-C-0058 

UNCLASSIFIED  REPORT 

ABSTRACT:  (U)  The  goal  of  this  study  was  to  investigate 
the  feasibility  of  fabricating  a  new  type  of  monolithic  focal 
plane  multipixel  staring  array  for  infrared  (IR)  radiation 
detection.  This  sensor  works  by  the  recently  postulated 
process  of  photofluxonic  conversion  in  which  an  infrared 
photon  creates  a  vortex-antivortex  pair  in  a  superconducting 
material.  The  novelty  of  the  new  sensor  lies  in  the  detection 
mechanism  itself,  and  in  the  way  the  vortex  (or  antivortex)  is 
trapped  in  the  device  allowing  individual  sensing  pixels  to  be 
read  out  without  direct  electrical  connections.  The  overall 
objective  of  Phase  I  was  to  determine  the  feasibility  of 
building  an  inductively  coupled,  non-contact  staring  array 
based  on  the  new  photofluxonic  detection  elements.  To  this 
end,  we  have  performed  a  series  of  experiments  on  a  single 
array  element,  a  2  mm  diameter  loop  of  YBa2Cu307,  and 
considered  pixel  requirements  for  constructing  an  array.  The 
new  system  offers  tremendous  potential  for  improvement 
over  current  infrared  detection  devices  via  both  the  enhanced 
detection  sensitivity  and  the  greatly  reduced  system 
complexity. 
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(U)  LNTlSb  as  a  Long-Wave  Infrared  (LWIR)  Material: 
Defects  and  Transport  Properties. 
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FEB  95  67P 
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CONTRACT  NO:  N00014-93-C-0091 
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ABSTRACT:  (U)  We  have  evaluated  three  III-V 
semiconductor  alloys  as  possible  candidates  for  future 
long-wave  infrared  (LWIR)  detector  materials.  The  cohesive 
energies,  elastic  constants,  band  structures,  electron  nobilities, 
and  phase  diagrams  are  calculated  and  are  compared  to  those 
of  Hg  sub  1-xCd  sub  xTe  (MCT)  alloys.  All  three  of  these 
III-V  alloys  have  their  band  gap  change  from  negative  to 
positive  values  as  the  alloy  composition  x  decreases  from  1  to 
0.  The  X  values  for  the  0.1-eV  gap  are  estimated  to  be 
0.67,0.15,  and  0.08,  respectively,  for  ITP,  ITA,  and  ITS. 

While  ITP  and  ITA  form  stable  zincblende  solid  solutions  for 
all  alloy  compositions,  zincblende  ITS  is  stable  only  for  a 
range  of  x  less  than  0.15.  ITP  and  ITA  have  considerably 
larger  cohesive  energies  and  elastic  constants  than  does  MCT, 
indicating  that  they  are  structurally  robust.  At  a  0.1  eV  gap, 
the  band  structures  near  the  gap  and  the  electron  nobilities  in 
ITP,  ITA,  and  ITS  are  also  found  to  be  comparable  to  those  of 
MCT.  Because  the  lattice  constants  of  ITP  and  ITA  are  less 
than  2X  larger  than  the  respective  values  in  InP  and  InAs,  the 
latter  should  provide  natural  substrates  for  the  growth  of 
active  LWIR  alloys. 
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